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SELENIUM PHOTOCELLS 


for use with direct reading instruments or valve amplifiers. 


TEST TUBE COLORIMETERS 


for Haemoglobin and other determinations. 


- REFLECTOMETERS 


for comparing reflectance from all surfaces. 


EVANS ELECTROSELENIUM 


LIMITED 
BESHOP’S STORTFORD, HERTS. 
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ADVERTISEMENTS 


in the 


PROCEEDINGS 


Rates for one insertion 
Page - £10 0s. Od. t-page - £2 12s. 6d 
$-page- £5 2s. 6d. $4-page - £1 7s. 6d. 
10°Jo and 20°/o reductions for series of three and six, 
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Applications for space should be made to 
THE PHYSICAL SOCIETY 
1 Lowther Gardens, Prince Consort Road, 
London S.W.7. 


OPTICAL BENCH 
PHOTOGRAPHIC ENLARGER 


PHOTOGRAPHIC COLOUR 
FILTERS 


OPTICAL FLATS 


THE 
LEECH OPTICAL Co. Ltd. 
WHETSTONE LEICS. 


Telephone: NARBOROUGH 2213 


TRANSFORMERS 
AND CHOKES 


to meet any practical specified requirement for 


SMALL 


Research and Experimental 
Purposes 


We have considerable experience in the design 
and manufacture of special transformers and 
chokes for a wide variety of purposes and in- 
clude amongst our clients many well-known 
Research Laboratories, Government Institu- 
tions, Technical Colleges and Industrial Re- 
search Departments. 


On receipt of data giving the desired electrical 
performance and circuit conditions we will 
quote for prototypes or quantities without 
obligation and are in a position to give very 
quick delivery. Our prices are reasonable, our 
goods of high quality and our technical service 
unequalled. We also stock most of the con- 
ventional types of transformers and chokes. 


SOWTER TRANSFORMERS, 


27 Homefield Road, Wimbledon, 
S.W. 19. 
’Phone : WIMBLEDON 0244. 


PROJECTION EQUIPMENT 
MIcROSCOPY 


APLANATIC MAGNIFIERS 
PRISMS TO ANY SPECIFICATION MIRRORS OF ALL TYPES 


LEECH 


BiuRNsDSUNGeoGaA SES 
for the 


Paco: Color DoloN2Gs5 


Binding Cases for volume 57 (1945) and 
previous volumes may be obtained for 
6s., inclusive of postage. For 9s. the six 
parts of a volume will be bound in the 
publisher’s binding cases and returned. 


THE PHYSICAL SOCIETY 


1 Lowther Gardens, Prince Consort Road, 
London S.W.7 


BLOOMED SURFACES 
ALUMINISED SURFACES 


RHODIUMISED SURFACES 


London Office : 
94, HATTON GARDEN, E.C.1 
Telephone: HOLBORN 1752 


VITREOSIL 
APPARATUS 


This photograph—by the way, it won 
International Exhibition in 1939—shows 
some pre-war apparatus made to speci- 
fication in VITREOSIL pure fused silica, 
but you need hardly be reminded that 
this work can always be done for you. 


THE THERMAL SYNDICATE LTD. 


Head Office: Wallsend, Northumberland. 
London Depot: 12-14, Old Pye Street, 
Westminster, S.W.|1. 
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SIMPLIFIED 


BALDWIN MEASUREMENT 


INSTRUMENT COMPANY LIMITED 


Laboratory methods of determination requiring time AND 
and skill can often be replaced by direct-reading 
instruments that can be successfully used by the non- 
technical worker. Our development service for industry 


covers the field of Idea, Prototype and Manufacture. CONTROL 
Ask us about it. 


——— 


BALDWIN INSTRUMENT COMPANY LIMITED 
London Office: GRAND BUILDINGS, TRAFALGAR SQUARE, W.C.2 


Originators and makers of scientific instruments for measurement and control. 


for HEAVY DUTY _ 
ELECTRICAL AND INSTRUMENT 


Gresrer tensile, elastic and fatigue strengths F 
than ‘any other non-ferrous alloy, a higher conductivity than BERYLLIUM COPPER 
any of the bronzes and excellent resistance to corrosion and \ 
wear—these characteristics of Mallory 73 Beryllium Copper 
have made it first choice for instrument springs, diaphragms : 
and bellows, current-carrying springs, snap action switch Properties of 


blades, contact blades and clips. MALLORY 73 BERYLLIUM COPPER 


Supplied annealed or lightly cold worked, it has good | after heat treatment 

forming properties and is readity fabricated into springs : P 

and parts of complicated shape. A simple heat treatment Se a ee EEE OG 

then develops its remarkable properties. Ror “i ; aa 
Limit of proportionality 

Available as sheet, strip and wire, in a range of tempers to tons per sq. inch 47.50 


suit users’ requirements, and as rod, tube, precision rolled 


; ; . : ; , Fatigue limit 
hair-spring strip, and silver-faced contact bi-metal strip. 8 


tons per sq. inch + ]9-20 
Vickers pyramid hardness 
350-420 


Electrical Conductivity 
per cent 1.A.C.S. 23-25 


JOHNSON, MATTHEY & CO. LIMITED, 
Controlling MALLORY METALLURGICAL PRODUCTS LTD. 
73-83 Hatton Garden, London, E.C.| 
Telephone: HOLborn 9277 


Full details are given in our booklet, 
which will be sent en request. 


G.D.40 - 
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SCI 
ENTIFIC BOOKS 


Messrs H. K. LEWIS can supply from stock or to order any book on 
the Physical and Chemical Sciences. 


CONTINENTAL AND AMERICAN works unobtainable in this 
country can be secured under Board of Trade licence in the 
shortest possible time. 

SECOND-HAND SCIENTIFIC BOOKS. An extensive stock of books In 
all branches of Pure and Applied Science may be seen in this department. 
Large and small collections bought. Back volumes of Scientific Journals. 


Old and rare Scientific Books. Mention Interests when writing. 
140 GOWER STREET. 


SCIENTIFIC LENDING LIBRARY 


Annual subscription from One Guinea. _Detai's of terms and prospectus free on request. 


THE LIBRARY CATALOGUE revised to December 1943, containing a classified index 
of authors and subjects: to subscribers 12s. 6d. net., to non-subscribers 25s. net., postage 8d. 


Bi-monthly List of Additions, free on application 


Taiephones eusion a2 1 bY K LEWIS & Co. Ltd. 


Telegrams: ‘‘ Publicavit, 


Westcent, London” 136 GOWER STREET, LONDON, W.C.1 


Progress hy Quality 


AVIMO 
35 mm. CAMERA 


Designed and developed by 
Avimo for Standard Cathode Ray 
Tube photography, the camera 
illustrated is rapidly becoming the 
standard Continuous Film Camera 
in the leading research labora- 
tories. Avimo specialise in in- 
dustrial cameras for specific uses. 


DESIGN - DEVELOPMENT - MANUFACTURE OF OPTICAL MECHANICAL:ELECTRICAL INSTRUMENTS 


AVIMO LTD TAUNTON SOMERSET ENGLAND 
8440 D 


a2 
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VALVES FOR RESEARCH 


AND 
ELECTRONIC CONTROL 
Photo-cells, Ignitrons, 
Thyratrons,9 Rectittetgse 
Safety Gaps 


@ S 


VALVES FOR THE 
CENTIMETRE WAVE BAND 
Magnetrons, Klystrons, 
Crystal Rectifiers 


Please send details 


of your requirements to :— 


RUGBY | 


PANY LIMITED, RUGBY, ENGLAND. 
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Pee hem Model fire 300 
HIGH SPEED CAMERA 


Available on hire for high speed research 


work requiring accurate observation of 


< 


fast moving machinery, projectiles, etc. 


COST OF HIRE:- 
£15.15.0 PER DAY (including operator) 


(See eS 
WoeoV NSE EN: EME 


forth Circular Road, Cricklewood, N.W.2. 
Telephone : GlAdstone 6373 


Beimoes i MPLY..BE ~-ACCURATE: = 
in pH determination 


For speed and accuracy in pH determina- 
tion Marconi pH Meters, using the glass 
electrode, satisfy all requirements in 
industry and laboratory. Electronics 
provides a means of making spot tests at 
any time, independent of colour compari- 
sons—yet in so simple a fashion that 
neither special skill nor training is neces- 
sary. The pH Meters are available either 
for mains or battery operation, each pro- 
viding rapid readings accurate to ‘02 pH. 
A demonstration can be arranged to suit 
your convenience, 


For process control there are Marconi record on a twenty-four hour chart, Full particulars 
instruments which provide continuous are available on request; when writing, please 
visible indication of pH or a permanent describe the process involved. 


| MARCONI INSTRUMENTS LTD. 
MEASUREMENT FOR INDUSTRY 


ST. ALBANS, HERTS. ’Phone: St. Aibans 4323-6 Northern Office: 30 ALBION STREET, HULL. 'Phone: Hull 16144 
Western Office : 10 PORTVIEW ROAD, AVONMOUTH, BRISTOL. ’Phone : Avonmouth 432. 
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| | Glee ee = See oe eae 
| | | | | | ape seeeeee re SA PERFORMANCE 
This record ofanactual | The SUNVIC Propor- 
experiment shows how tioning Head is nor- 
the SUNVIC Proportion- mally used in conjunc- 
ing Head, ‘used with any tion witha SUNVIC Type 
Standard Toluene Regu- EA2/T Electronic Relay. 
lator, improves the This apparatus is fully THE THERMOSTATIC 
fineness of regulation, described in Technical CONTROM ESLER TS 
giving a _ stability of Publication EA11/48. SUNVIC CONTROLS LIMITED 
the order of 5 milli- Please request a free  gtanhope House, Kean Street 
degrees. copy today. London, W.C.2 
yori, I NT yy wl IN UT m1 
tit il ‘ial ‘iy ol ae Hi ual 
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| THE G-L-C CONSTANT VOLUME. 
GAS ANALYSIS APPARATUS 


A new apparatus, simple to operate, for use in the | 


accurate analysis of practically any gaseous mixture | 


@ A modification of the ‘Bone and Wheeler’’ apparatus, | 
developed in the Fulham Laboratory of the Gas Light & Coke’ | 
Co. as a result of many years’ experience. 


® Using mercury as the confining liquid, it is a precision 
apparatus essential in gasworks and where producer and 
other gas plants are in operation. 


PRINCIPLE: A water-jacket burette, an absorption pipette and two gas- 
heated reaction tubes containing specially prepared copper oxide and a 
copper oxide-iron oxide catalyst. 


ADVANTAGES: Accuracy, rapidity, simplicity, freedom from solubility, 
absorption and calibration errors, complete absorption assured by the 
use for each analysis of fresh reagents, etc. 


Grrr INand [ATLOCK La 


1938, 57, 388 


Ref. J.S.C.1., 
PI ly for Booklet GT1302/31 ; err: ; 
CASS APR ESOS | Established as Scientific Instrument Makers in 1826 
LONDON MANCHESTER GLASGOW EDINBURGH 
Kemble Street, W.C.2 19, Cheetham Hill Rd., 4 45, Renfrew Street, C.2 7, Teviot Place, ! 
BIRMINGHAM : 


STANLEY BELCHER & MASON LTD., CHURCH STREET, 3 
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THE “SHIRLEY” 
MOISTURE METER 


Originally developed for the rapid determination 
of the moisture regain of cotton, the instrument is now 
availabe for a wide range of materials, including wool, 
rayon, flax, jute, hemp, etc. 


Information concerning Meters suitable for 
testing other materials will be given on request. 


The instrument is both quick and accurate, giving 
instantaneous direct readings of Moisture Regain. 
Worked from A.C. Supply (consumption is only 25 
watts), the meter is so simple that an unskilled person 


can obtain correct readings. 


THE RECORD ELECTRICAL CO. LTD. 
BROADHEATH °: ALTRINCHAM : CHESHIRE 
’Phone : Altrincham 3221/2 "Grams : “ Infusion ’’ Altrincham 


London Office : 28 Victoria Street, Westminster, S.W.!. Phone: Abbey 5148. 


\V \V FOR NEARLY A CENTURY WRAY LENSES 
) MA\ HAVE ENJOYED THE CONFIDENCE OF SCIENTISTS 


SCIENTIFIC 
LENSES 


Vv 


THE LENS HERE ILLUSTRATED 
HAS BEEN SPECIALLY DESIGNED 
FOR USE ON CONSTANT SPEED 
CAMERAS TO RECORD ON 35 mm. 
FILM FROM A CATHODE RAY TUBE. 
IT IS AN ANASTIGMAT OF 2” FOCAL 
LENGTH, HAS AN APERTURE OF F.!, 
AND IS CORRECTED FOR AN OBJECT 
TO IMAGE RATIO OF 4:1. 


WRAY (OPTICAL WORKS) LTD. BROMLEY, KENT 


vili 


SESS 


VISUAL TWO-DIMENSIONAL 
delineation of any recurrent law. 


RELATIVE TIMING OF EVENTS 
and other comparative measure- 
ments with extreme accuracy. 


PHOTOGRAPHIC RECORDING 
of transient phenomena. 


SIMULTANEOUS INDICATION 
of two variables on a common 
time axis. 


Completely embracing all the above func- 
tions, the unique Cossor DOUBLE BEAM 
Cathode Ray Tube, as embodied in Model 
339 Oscillograph, is inherently applicable 
to all problems arising in the 


RECORDING, 
INDICATING and 
MONITORING 


of effects and events in Electrical engineer- 
|} ing ; and in Mechanical engineering when 
the effects can be made available as a vol- 
tage. Recurrent traces are studied visually 
and transients may be photographically 
recorded with Model 427 Camera. 


A. C. COSSOR LTD., 


Instruments Dept. 
Cossor House, London, N.5 


*Phone : CANonbury 1234 (30 lines). 
’Grams : Amplifiers Phone London. 
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CARPENTER HIGH SPEED) 
POLARIZED RELAY 


Main features of Standard model : High 
speed. Short transit time. Negligible dis- 
tortion up to 300 c/s. Contact gap a | 
function of input power, hence small dis- | 
tortion almost down to failure point. High 
contact pressures. No contact chatter. 
High sensitivity—robust operation at | 
5 mv.a. at 100 c/s. or 0-2 mw. D.G 
Great ease of adjustment. Magnetic bias | 
adjustment giving absolutely smooth con- | 
trol. Balanced armature—hence immunity 
to considerable vibration and no positional | 
error. 


Full details and Test data available on request. 


TELEPHON E 
MANUFACTURING 
co. LTD. 


Contractors to British - Dominion - Colonial and * 
Foreign Governments. 


HOLLINGSWORTH WORKS, 
DULWICH, LONDON, S.E.2I 


Telephone 
GIPsy Hill 2211 (10 lines). 
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DRAYTON 


“Hydroflex’ 
METAL BELLOWS 


These bellows are formed from the 
initial tube in one gradua! con- 
tinuous operation, resulting in a 
uniformity of wall thickness un- 
attainable by any other method. 
The tough resilient product is 
tested to many times the maxi- 
mum rated working pressure dur- 
ing formation. Customers’ end 
plates can be fitted prior to form- 
ing so that the soldered joint is 
also pressure tested. 


For Gland Seals ; Refrigeration 
Control; Thermostatic and 
Pressure Operated Devices, etc. 


A PACKLESS GLAND WITH 
OUTSTANDING ADVANTAGES 


@ Every Gland uniform @ Absolutely reliable in 
in life and performance. operation. 


@ Pretested during @ Provision of closed 
forming. end eliminates one joint. 


Root Diam. 3/8’ to3” Outside Diam. 9/16” to 43” 


SPR ES PA EY ST ET EE ST 


RAYTON REGULATOR & INSTRUMENT CO. LTD, 
EST DRAYTON West Drayton 2611.  MIDDX. 


EEE 


WESTINGHOUSE] @) 


CONSTANT 
VOLTAGE 


A GcaDEG: 


In the following systems, automatic 
compensation is provided for the effects 
of mains fluctuations and/or load 
variations. 


The “Stabilistor”’ 


This device provides an A.C. voltage of 
good waveform, maintained within limits 
of +2% for all loads between no load 
and full load and any input between 190 
and 260 volts 50 c.p.s. single-phase A.C. 
Closer limits are obtainable if either the 
input or the load are constant. Response 
is almost instantaneous and efficiency is 
about 85%. 


Where the load is constant, a_ special 
frequency-compensated model is available, 
which will hold the output voltage constant 
to within a few parts in 1,000 in spite of 
variations of + 6% in mains voltage and 
simultaneous or independent variations of 
+ 2% in mains frequency. 


““Westat”’ 


A special fully automatic battery charger, 
which maintains the battery voltage within 
fine limits despite simultaneous variations 
of mains and load. Action is instantaneous 
and, after a power supply interruption, 
the Westat automatically provides a heavy 
charge until the battery voltage reaches its 
normal value. No moving parts, saturated 
chokes, gas discharge tubes or barretters. 


“Noregg 39 


Similar to the Westat but for direct operation 
of D.C. appiratus from the A.C. mains. 
High efficiency and power factor; low 
percentage ripple. 


Write for publications E.E.2 and D.P.\1L to Dept P.S. 


Westinghouse Brake & Signal Co. Ltd. 
Pew Hill House, Chippenham, Wilts. 
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OPEN{CIRCUIT VOLTAGE\e 
OF A¥67mm. MEGATRON## 
PHOTO CELL ij 


~“\\\ 


INSTRUMENTS 


FOR INDICATING: RECORDING: CONTROLLING 


es 
TEMPERATURE = 8, MEGATRON LTD. 


We shall be glad to supply full particulars 115A FONTHILL ROAD, 


BOWEN INSTRUMENT Co. Ltd. AG 
9/13, NEWTON ROAD - LEEDS, 7 


Telephones: 41036/7 


Archway 3739 


RESISTANCES 


Although present circumstances 
render it difficult for us to give our 
pre-war service to all customers, we 
are still working in their interests. 


New materials and manufacturing 
processes, which we are now using 
to increase output, also contribute 
in large measure to improved 
performance and reliability of our 
products. Thus, when normal times 
return, all users“of Berce 
Resistances will benefit by our 
work to-day. 


fHE BRITISH ELECTRIC. RESISTANCE GO. LTD. 
QUEENSWAY, PONDERS END, MIDDLESEX 


: Telephone: HOWARD 1492. 
Telegrams VITROHM, ENFIELD.” ~ 


(R.1) 
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.F. OSCILLATORS ; U.H.F. GENERATORS ; Q-METERS ; 
ALVE VOLTMETERS ; NOISE-METERS ; MULTI-RANGE 
ETERS ; AUTOMATIC OHMETERS ; STRAIN GAUGE 
MPLIFIERS ; “CATHODE RAY. OSCILLOSCOPES, etc. 


LELAND 


INSTRUMENTS LTD. 


Electronic Engineers & Consultants 


tA NEW RANGE OF 


‘BEAT 
FREQUENCY 
OSCILLATORS 


MODEL 27 Mk. II-B 


is a compact portable mains- 
driven instrument covering the 
range 20-21,000 c.p.s. with 125 
milliwatts output (25 v. across 
5,000 ohms and 8v. across 
500 ohms). Distortion is less 
than 2:5 per cent. The scale is 
directly calibrated and direct 
reading. A ‘‘ magic eye ’’ indi- 
cator is provided for zero setting. 


® 
MODEL 27 Mk. Ill 


is similar, but has an A.C. output 
voltmeter and pilot lamp, instead 
of the ‘‘ magic eye ’”’ 


ce 


® 
MODEL 27 Mk. IV-A, 


illustrated, is a two-range instru- 
ment which provides a very open 
scale at the lower frequencies. 
The upper half of the tuning dial 
covers 10-21,000 c.p.s., while 
the lower half covers 10-1,000c.p.s. 
In all other respects this model is 
similar to the Mk. II-B. 


®@ 
OTHERINSTRUMENTS 


covering frequencies between 
5 c.p.s. and 5 m.c.p.s., are also 
available. Engineers are invited 
to apply for full data, stating the 
frequency range and output they 
require. 


e 
PARTICULARS 


of this range and other instru- 
ments are available in 1946 
catalogue. 


| JOHN STREET, BEDFORD ROW, LONDON, W.C.!. Te! : CHANCERY 8765 


Cough -Spangle  Gerris . Boonton » fallantine  Mipiyn 


XU 
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-LABGEAR ELECTRONIC VIEWING UNIT| 


A compact 
Cathode Ray 
Tube Unit 
complete with 
power supply 


Weight 
1234 Ib. 
Screen size | 
Size 


23 in, | 
3 ; 
64X 553 X 12 in. | 


Designed for application as a Modulation Indicator or Frequency Comparitor in the Laboratory. 
Specification and Details supplied on request 


Works : i A B ( E A R Head Office : 
WILLOW PLACE, 


FITZROY LANE, 
CAMBRIDGE INSTRUMENT AND ELECTRONIC CAMBRIDGE 
ENGINEERS 


& MUMETA| 


= MAGNETIC SCREENS 


MUMETAL is a 
renown 
S— Passed for the shieldin ed TELCON alloy unsur- 


& of delicate equipment 


: or alternati : 

ee ve During Wartime MUMETAS Base te 
ensively employed for a2 : as been 

other eens in Radar and 


The experience 


inquiries are invited. 


AEEREEES THE TELEGRAPH CONSTRUCTION & MAINTENANCE CO. L | 
Founded 1864 i 


Head Office: 22 OLD BROAD ST., LONDON, E.C.2. Tel: LONdon Wall 3! 
- Enquiries to: TELCON WORKS, GREENWICH, S.E.10. Tel: Greenwich I 
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HILIPS GEIGER-COUNTER X-RAY SPECTROMETER 


. The new Philips Geiger-Counter X-Ray 
Spectrometer is a revolutionary development 
based on newly applied principles of diffraction 
technique. Using a direct-reading method, no 
film is required. The speed of analysis is increased 
and greater accuracy obtained. The Spectrometer 
is simple, rapid and easy to operate. It provides a 
fast medium for the qualitative and quantitative 
analysis of crystalline and certain amorphous 
materials, quickly identifying chemical substances 
present and their state of chcmical combination. 
Under optimum conditions of resolution, diff- 
raction angles are obtained with an accuracy 
of 0:03 degree. The reflected beams are located 
by means of a highly sensitive Geiger-Counter 
tube mounted on a movable arm that traverses 
the scale. Radiation intensity is integrated and 
totalized by an electro-mechanical counter, while 
average intensity is indicated by a micro-ammeter, 
both connected to the Geiger-Counter through 
electronic circuits. 


HILIPS X-RAY DIFFRACTION APPARATUS 


. This apparatus constitutes the ideal flexible 
diffraction machine. With a power unit rated at 
60 kvp. 50 ma. any of the modern techniques can 
be readily employed. The full wave rectification 
means a longer tube life because the tube filament 
runs cooler, and for equal target loading the 
filament emissicn is only about half of that 
compared with self-rectified or half-wave recti- 
fication. | Furthermore, contamination of the 
tube target, due to deposits of tungsten, is thereby 
greatly reduced. This results in a clean tube 
target, no contaminated radiation to complicate 
the pattern on the film, and longer, more efficient 
tube life. Among the many applications for this 
machine are :—Chemical Indentification . Phase 
Analysis . Solid Solubility Study . Single Crystal 
Studies . Particle Size Measurements . Metallo- 
graphic Studies . Molecular Weight Determination 
Structure Identification . Stress Analysis . Fibre 
Studies . Process Control . Orientation Studies. 


PHILIPS = METALIX 


ILIPS LAMPS LTD, CENTURY HOUSE, SHAFTESBURY AVENUE, LONDON, W.C.2, 
(149) 


Xiv 
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Research in Applied Photography 


For many years Ilford Limited have played a leading 
part in rendering photographic service to scientists in 
the Universities and other research centres, both here 
and overseas, and in industry. They have taken pride 
in these services as part of their contribution to Science 
itself—a contribution which has often been out of all 
comparison to the actual quantities of sensitive materials 
involved. 


The demands for these services have grown apace. They 
have been intensified rather than reduced by the war, as 
photographic modes of investigation have found unex- 
pected uses. Ilford Limited have therefore expanded 
their organisation to cope with such demands by extend- 
ing the scope of their Technical Service Department, and 
by instituting a new Technical Research Department. 
While the former is prepared to deal with requests for 
technical data or assistance in problems arising in the use 
of standard products, the latter is concerned with new 
problems or novel applications which may involve new 
sensitive materials or photographic techniques. 


Among the fields to which Ilford Limited have made 
outstanding contributions are: 


SPECTROGRAPHY 
PHOTOMICROGRAPHY 
ELECTRON MICROGRAPHY 
CATHODE-RAY OSCILLOGRAPHY 
PLATES FOR NUCLEAR RESEARCH 
SELECTIVE AND NEUTRAL FILTERS AND WEDGES FOR 
RESEARCH PURPOSES GENERALLY 


Special ILFORD products for these purposes are backed 
by an unsurpassed range of plates, papers and films for 
general work. 


Consult ILFORD LIMITED about your photographie 
problems and the experience gained by us and by you 
will add to the common pool of scientific knowledge and 
technique. E 


ILFORD LIMITED - ILFORD - LONDON 
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[eal Lies Uscillator 


Ste 


BEAT FREQUENCY | 
OSCILLATOR 


TYPE 


AF 200 
FREQUENCY RANGE HARMONIC CONTENT 

50 to 20,000 c/s. Less than 3% over main portion of the 
ACCURACY frequency range at full rated output. 

+ 1% + 2 ¢/s. R.F, CONTENT 

OUTPUT Less than 025% at full output. 
2 watts into 600 or 10 ohms. DRIFT 

HUM Less than 10 c/s a day after first half- 

Less than 0°5% at full output. hour of operation. 


PRICE £50 


FURZERILL <i> pannATOAIES 


;OREHAM WOOOG-H EATS 


TELEPHONE : ELSTREE 1137 
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MULLARD-BT: 


ELECTROMETRIC 
TITRATION APPARATUS 
TYPE-E920 


600r T T , + CO 
H 
750 
re | 760 
| 
Fr PAS on On DUIS reed 650 
WITH N/20 TER ROUG AMMONIUM SULPHATE 
ELEC TRODLS, BRIGHT PLATINUM AND) 
coo! SATURATED CALOMEL 600 
2 Seer ee mV 
t 500 
ic anata 1 450 
Pe tes Res ae eS a | oa 3 400 
“CS FeWs 
CCS Fe SO4 
Titration of Vanadium in Steel with N/20 Ferrous 
Ammonium Sulphate Electrodes: Bright Platinum and 
Saturated Calomel. 


Tuis apparatus has been designed to supply the need 
for a robust, self-contained potentiometric titration 
apparatus which is sufficiently simple to be used in routine testing by unskilled 
operators and is yet capable of meeting the requirements of the industrial research 
chemist. 

The apparatus is operated from the 50-cycle supply mains and the end point is 
detected by a ““ Magic Eye” indicator. Thus there are no batteries to be replaced 
and no delicate galvanometer to be damaged by mechanical shocks or electrical 
overloads. A special circuit eliminates all possibility of drift during a titration, 
and changes of mains supply voltage do not give rise to any inaccuracies. 


OTHER TYPES OF MULLARD APPARATUS | 
HIGH AND LOW FREQUENCY A.C. BRIDGE EQUIPMENT. GRAIN MOISTURE METERS.) 


CATHODE RAY OSCILLOGRAPHS. CONDUCTRIMETIC APPARATUS. BEAT-FREQUENCY OSCILLATORS | 


i 


Bes 


MULLARD 


RGN OE eS: SESS SS — ue ooh Pil 
WIRELESS SERVICE €O. LTD., CENTURY HOUSE, SHAFTESBURY AVENUE, W.C i = 
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MULLARD 


“Ticonal” Permanent Magnets _ 


FOR INDUSTRY 


MULLARD “TICONAL" magnets are the most efficient 
permanent magnets available. A revolutionary process of 
manufacture, including a new alloy and heat treatment in 
a magnetic field, was developed in our own laboratories, 
resulting in magnets of thrce times the performance of 
any other commercial magnet produced at that time. 


‘’Ticonal”’ permanent magnets have set a new high standard 
and enable magnetic circuits to be designed with high 
efficiency, resulting in small physical dimensions and weight 
and in many cases an appreciable saving in cost. 


That is why more and more manufacturers in the electrical 
and radio industries are specifying “ ‘Ticonal ”? permanent 
magnets in their products. 


Below are a few typicai applications :— 


LOUDSPEAKERS e MOTORS & GENERATORS 
MICROPHONES e TELEPHONE HAND GENERATORS 
PICKUPS e COMPASSES e MAGNETIC FILTERS 
CYCLE DYNAMOS e MAGNETRON FIELDS 
TEMPERATURE CONTROL e SPEEDOMETERS 


ELECTRICAL INDICATING & RECORDING INSTRUMENTS 


The experience of our engineers is at your service on any problems connected 
with the design of magnetic circuits. 


THE MULLARD WIRELESS SERVICE CO. LTD., 


|(MAGNET DIVISION) 


CENTURY HOUSE, SHAFTESBURY AVENUE, LONDON, W.C.2. 
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THE KENT 


MULTELEC 


is a strong and sturdy instrument which, upon activation at two- 
second intervals by any minute change of potential, sets up a 
powerful servo-motor force for controlling any desired form of 
correction, yet without imposing any restraint upon the sensitivity 
of the instrument. Furthermore, it indicates and records. 


VARIOUS USES: 

The Kent Multelec is extensively used in various industries for 
recording, indicating and controlling temperatures, its potentiometric 
basis ensures a fine degree of accuracy under industrial conditions. 


It is also used in the control of pH value in various processes, and 
for the correction of boiler feed water by measurement of 
electrical conductivity. 


It will be recognised that the Multelec can be adapted to many 
forms of controlling, recording or indicating besides the above. 


George Kent, Ltd., also supply meters for accurate measurement 
of the flow of steam, water, gas, air, oil or any other fluid, 
whatever the quantity or pressure. 


The Multelec booklet (No. 897) is available on request. 


GEORGE KENT, LTD. LUTON & LONDON 
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Portable test set Type -32-A. 


This set constitutes a complete self-contained equipment for measurement of conductor resistance, dielectric 
resistance, capacity and for fault location on telegraph and telephone lines and cables. 

The set contains the following comprehensive equipment: BRIDGE NETWORK comprising variable ratios 
and 4-decade measuring arm—TWO GALVANOMETERS (Pointer type and sensitive mirror galvanometer) 
with 10-position- UNIVERSAL SHUNT, SHORT CIRCUIT KEY and REVERSING SWITCH — RESISTANCE and 
CAPACITY STANDARDS for determining galvanemeter constants—66-volt. TESTING BATTERY, BATTERY 
RESISTANCE, REVERSING KEY, MILLIAMMETER-PLUG and COMMUTATOR for selecting test circuits. 


Mounted in leather bound teak case with all accessories, spare parts and repair kit — Illustrations show the 
instrument closed and sct up for use. 


DIMENSIONS Closed, 24 in. x 124 in. x 112 in. deep. 
WEIGHT 56 Ibs. 


A> DIETS GoRSISP Sis ILVeE BUELL Eaten Neale aL 
FOR OVER 60 YEARS DESIGNERS AND MAKERS OF 
PRECISION INSTRUMENTS. 


MUIRHEAD & CO. LTD., ELMERS END, BECKENHAM, 
KENT. Telephone: Beckenham 0041-0042 


BEDS (UUPSP ESISE Da Orn REQUH 
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Testing Time 
Ahead... 


RECONSTRUCTION, as we now know too well, means something other than 


enjoyment of the non-existent fruits of victory. The amenities of a world at 


well-earned peace are not for us until we have replaced the ravages of war with 
the necessities of life and the realities of universal peaceful intent for victors, 
victims and vanquished alike. 
That is speaking collectively . 
far in the six years of ceaseless toil, urged on by dire necessity and poril. 
are not resting now. We are still pressing on, pressing into the service of those 
engaged on rebuilding the body and soul of a whole world the knowledge gained, 
the advancements perfected, the skill and craftmanship that outmatched the 


.. For ourselves, we learned much and progressed 
We 


efforts of our enemies. 
\\\ 


PRECISION ELECTRICAL MEASURING 


SS 
gy 


Regd. Trade Mark 
Sole Proprietors and Manufacturers : 


XxX1 


“Avo ’’ Electrical 
Testing Instruments set a 
standard by which other 
instruments are judged. 


THE UNIVERSAL 
AVOMETER 


THE" AVO. 7 
ALL-WAVE OSCILLATOR 
THE ‘‘ AVO”’ 
VALVE TESTER 
THE ‘‘ AVO”’ 
TEST BRIDGE 
THE UNIVERSAL 
AVOMINOR 
THE D.C. AVOMINOR 
THE ““‘ AVO” 
LIGHT METER 
THE ‘“‘ AVO”’ 
EXPOSURE METER 
BTC., ETC. 


Catalogues on application. 


INSTRUMENTS 


THE AUTOMATIC COIL WINDER & ELECTRICAL EQUIPMENT CO. LTD. 
’Phone: ViCtoria 3404/8. 


Winder House, Douglas Street, London, S.W.|. 
i IAF SE STA TESS DT AE SI Ns DEP CAI STIRS ARO 8 se 8 SNA SA GLE FS TN a TS 


Tie 
NEW MULTITONE UNIVERSAL 


POCKET HEARING AID 


The New Pocket Aid makes use of the new cascade circuit (Patent applied | 
for) which we have developed. Two valves only are used, a heptode as a two- | 
stage amplifier and a H.F. pentode as an output triode. The new circuit has 
enabled us to produce a hearing aid of outstanding performance, at the very | 
moderate price of 12 guineas. 


The new instrument will give excellent hearing to all cases of deafness | 
responsive to instrumental help. It can give, if required, higher amplification | 
than any other pocket aid we have produced in the past, whilst its controls can | 
be adjusted to make it equally helpful to those needing clarity of reproduction | 
rather than great volume. It can be used without alteration with any type | 
of earpiece, including crystal or bone-conductor, and is constructed throughout | 
with proved components and valves available the world over. | 

The instrument incorporates many new features, one of which regulates | 
the high tension current to the particular requirements of the individual user, | 
thus ensuring longest possible life for high tension batteries. In this direction 
it is particularly economical—-the cost of replacements for a severely deaf person | 
being estimated at Is. 6d. per week. 


The batteries used are of the new layer built type—more compact and of | 
greater capacity than any obtainable hitherto. They are all fitted with inter- | 
national connections and sockets, and are obtainable all over the world. 


Robust in construction, reliable and economical in use, the Universal 
Pocket Aid is small and light enough to be easily worn on the person. It is 
housed in an attractive plastic case, and connected with the earpiece and 
batteries by a single lead. 

SPECIFICATION. 
Size : 52 in. by 2% in. by @.in. Weight : 94 oz. 


Three-stage amplifier with crystal microphone. Button-base international valves. 
Low tension drain 90 m.a. High tension drain from + to 24 m.a. Maximum 


voltage amplification 4,000. High tension battery 45 volts layer-built. Low 
tension |4 volts. 


All our aids may be had on a week’s home trial, free of charge. 


There are Multitone Agents in most cities in Great Britain and in many 
countries abroad. 


MULTITONE ELECTRIC Co. Ltd. 
92, New Cavendish Street - - London, W.| 
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ANODE VOLTAGE. 


ANODE & SCREEN CURRENTS IN: 


THE KT 8 
The OSRAM KT8 is a transmitting beam tetrode particularly 
suitable for use as an R.F. amplifier, oscillator or frequency 
multiplier. [ts outstanding features include :— 


A useful output down to wavelengths of the / 
order of six metres. ja Ye 


“Aligned grids ’’ resulting in low screen current, 
largely removing the limit usually set by screen 
dissipation. 


Low interelectrode capacities normally making 
neutralising unnecessary. 


&£G.C. Osram 


CATHODE RAY TUBES VA EVES 


Osram 
PHOTO CELLS 


Advt. of The General Electric Co., Ltd., Magnet House, Kingsway, London, W.C.2, 
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APPARATUS FOR 


CHEMICAL ANALYSI 


Mieke OWA S Gi 
special interest to the 
industrial chemist, both for 
research work and for routine 
analyses. They are described 
in the publications mentioned 
below, which contain much 
valuable information. 


POLAROGRAPH 


For the rapid electro-chemical analy- 

sis of solutions, enabling traces of 

metals and other substances to be 

determined by the Heyrovsky drop- 

i ping mercury electrode method. 

tl Full notes on technique, with an 
extensive bibliography, in 


LIST 109-L 


An_ indicating polarograph (the 
‘“Valtamoscope’’) is also available. 


| TITRATION APPARATUS 
(Illustrated) 

f For applications where a titration 
method is desirable. Especially useful 


for determining end-points in co- 
loured solutions. Described in 


SHEET 247-L 


pH METERS 


The full range of instruments for 

measuring, recording and controlling 

pH in all branches of research and 
industry is described in 


LIST 108-L 


Any or all of these useful 
publications gladly sent on 
request 


13, GROSVENOR PLACE, LONDON. SWI... 
WORKS: LONDON & CAMBRIDGE. 
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GEOMAGNETIC SECULAR VARIATIONS AND 
SURVEYS 


By. J._A. FLEMING, 


Director, Department of Terrestrial Magnetism, 
Carnegie Institution of Washington 


Third Charles Chree Address, delivered 6 December 1945 


ABSTRACT. The secular variations of the Earth’s field call for frequent determinations 
of the geomagnetic elements at many selected stations on land and at sea. So far, general 
world magnetic surveys have of necessity been restricted to the surface of the Earth. 
Previous surveys, including those of the Department of Terrestrial Magnetism of the 
Carnegie Institution of Washington, have led to great improvement in our knowledge 
of the distribution of the field, especially during the past few decades. During the war, 
much attention has been devoted to the study and analysis of data well distributed over 
the Earth’s entire surface at more than 10,000 stations. These have resulted, for the first 
time, in the preparation of accurate isoporic charts, that is, charts of equal annual rate of 
change for magnetic declination, inclination, and the horizontal, vertical, eastward, north- 
ward and total components of the field for the four epochs of 1912-5, 1922-5, 1932:5, 
and 1942-5. The motions of the maximum and minimum isoporic foci during these four 
epochs indicate the complexity of the secular changes and interpretations. 

Isomagnetic charts based upon surface observations must always be limited, so far as 
faithful depiction of the field is concerned, because of the impracticability of obtaining 
observations at an infinite number of stations and of deductions for both the regular and 
irregular changes in the field. Progress in instrumentation during the war on the applica- 
tions of the geomagnetic field of the Earth have produced improvements which make 
feasible the early realization of magnetic surveys by airplanes at several different levels 
proposed by Professor Sydney Chapman in the first Charles Chree Address. Some of the 
potentialities, possibilities, and needs for intense national and international coordination in 
magnetic surveys by plane are reviewed. 


Sa UN ROD Cay LON 


AM humbly grateful to your Society for making me the third recipient of the 

Charles Chree Medal and Prize. ‘The award is, I take it, largely a recognition 

of the geomagnetic researches which my colleagues and I have made at the 
Department of Terrestrial Magnetism during more than four decades with the 
generous support of the Carnegie Institution of Washington. ‘T’he two previous 
awards were to scientific colleagues and friends (see Chapman, 1941, and 
Schonland. 1943), whose work has constantly inspired investigators everywhere, 
and especially those of us in the United States who have been associated with the 
Department of Terrestrial Magnetism. The award affords further evidence of 
that scientific co-operation of British and American physicists and geophysicists 
which has been so vital in forwarding our common war effort, particularly in 
geomagnetic applications. 


PHYS. SOC, LVIII, 3 15 
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Progress in geomagnetism owes much to Charles Chree as Superintendent of 
the Magnetic and Meteorological Observatory at Kew from 1893 to 1925, and 
thereafter until his death in 1928. He and his contemporaries—Schmidt, 
Bauer, Van Bemmelen, and Moos—built well on foundations laid by their eminent 
predecessors in the task of portraying the complex phenomena of the magnetic and 
electric fields of the Earth through surveys on land and sea, investigations in the 
laboratory, and theoretical hypotheses and interpretation. As a colleague and 
friend, Chree influenced greatly the work of Bauer in 1899 through the expansion 
of geomagnetic research in, and survey of, the United States Coast and Geodetic 
Survey, and again in 1904 through the realization of Bauer’s plan for a Department 
of Research in Terrestrial Magnetism submitted to the Trustees of the Carnegie 
Institution of Washington. His wise counsel carried great weight and played a 
large part in the subsequent development and growth of the Department. Chree’s 
share in international and national progress in our field is well known to all of you. 
His commonsense treatment of factual information gathered in all parts of the 
world—especially in polar regions—was outstanding; it has contributed much to 
applications of geomagnetic research touching so many lines of human peace-time 
and war-time endeavor. 

The foundation and provision by Miss Jessie Chree of the Charles Chree 
Medal and Prize through the Physical Society not only serve tohonourthememory 
of Chree’s accomplishments but also to stimulate future workers and progress in 
our fields. 

The first Charles Chree Address, in 1941, by Professor Sydney Chapman gave 
a scholarly general account of the magnetic field of the Earth. I have chosen for 
my subject Geomagnetic secular variations and surveys, and will attempt to review 
at greater length world magnetic surveys and potentialities so prophetically 
outlined in Part III of the first Address. 


§2. GEOMAGNETIC SECULAR VARIATIONS 


Several centuries of speculation and research on the old and challenging 
problem of the origin of the Earth’s magnetism and its changes have yielded as yet 
no adequate explanation. Much has been learned from past attempts at explan- 
ations, and we have also added, in the course of time, much new and rich material 
relating to the description of the phenomenon. This material has illuminated 
greatly the picture of what we are trying to explain, but understanding of the 
cause of geomagnetism is emerging relatively slowly. 

In spite of vast improvements in our knowledge of the geomagnetic field, that 
knowledge remains insufficient for various theoretical and practical problems. 
For example, Schrodinger (1943) has recently attempted to use the data from 
geomagnetism to test some theories of physical fields which he has developed. 
For these tests he required information on the magnitude of the external field 
and of the non-potential field. It was important to know the magnitudes of these 
quantities to a moderate degree of precision, but those who are working in geo- 
magnetic research recognize that we cannot even be sure, on the basis of existing 
observations, that there is an external field or a non-potential field. 

The magnetic field varies with time (Fleming, 1938a; Barker and McNish, 
1938). Its delineation at any instant is hence difficult because it is not practical 
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to undertake its simultaneous measurement at all points of the Earth’s surface. 
There are some 70 magnetic observatories (figures 1 and 2) where continuous 
simultaneous measurements of the field are made, but the complexity of the field 
and its changes with time do not permit satisfactory interpolation of values in 
intervening regions. The records at observatories must be supplemented by 
observations at many isolated points. The density of observations, however, per 
given area of survey is uneven, and great oceanic areas are sparsely represented. 
Measurements at these stations on land or sea are usually made on one day only, 
and the resulting data form the basis for isomagnetic charts. 
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Figure 1. Geomagnetic observations, Eastern Hemisphere, functioning continuously or 
at intervals, 1905-42, and recommended additional locations. 


The world magnetic survey of the Carnegie Institution of Washington from 
1905 to the beginning of the war improved knowledge of the general field and its 
changes, which are slow when compared with the life-span of a human being but 
are very rapid when compared with astronomical or geological scales of time. Nor 
are these secular changes small. During the past half-century a gradual change 
amounting to 10°% of the maximum value of the Earth’s field has occurred off the 
Guinea coast. By way of further illustration, consider a few indications based on 
long series of observations. The actually observed changes in magnetic direction 
at London and at several other places as determined from earliest known measure- 
ments, including those for 1945, are shown in figure 3. ‘The figure shows also 
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notable changes in the position of the magnetic equator since 1700. The pro- 
gression in direction of the compass by centuries is non-uniform and notably large. 
The general aperiodic character of secular change in direction is confirmed by the 
measured directions of magnetization of particles locked in annual varved deposits 
in lakes and age-old deposits in oceans. 

Radical changes in the estimated position of the lines of zero magnetic declin- 
ation (agonics) are noted since the year 1500. The line along which the compass 
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Figure 2. Geomagnetic observations, Western Hemisphere, functioning continuously or 
at intervals, 1905-42, and recommended additional locations. 


points true north (figure 4) has shifted westwards from Suez to Chicago in 445 
years, only a short interval in the Earth’s history. 

Clearly, secular change involves the general magnetic field of the Earth, and 
any attempt to explain either phenomenon must also account for the other. 

There are two hopeful present approaches. One is the derivation of more — 
adequate and complete data on the observed large-scale changes in the Earth’s 
magnetism. The other is the improved description of the permanent field itself. 
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Figure 3. Secular variation of magnetic equator in Western Hemisphere from 1700 to 1945 
and of magnetic declination at London, Eastport, San Francisco, Callao, and Rio de Janeiro, 


from earliest observed values. 
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Figure 4. Progressive change (secular variation) of lines of zero magnetic declination (agonics). 
(The dots indicate the side of the agonic lines on which the declinations were easterly.) 
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The special type of the research required calls in a unique degree for world-wide 
coordination of data and experiments. No single, well-planned experiment or 
observation can supply the solution. Observations must be made in all parts of 
the world and continued for a long period of time. Techniques for the organiza- 
tion and interpretation of the data must be developed and experimental research 
must be conducted along lines which will supply information on basic physical 
problems related to this subject. 

In the reduction of magnetic observations to epoch, the aim is to determine at all 
stations (at epochs not too far removed from those at previous times of observation) 
the intensity and direction of the Earth’s field freed from the effects of small 
superposed fluctuations in intensity. For practical reasons the problem of 
describing the Earth’s surface field is simplified by mapping only the main or 
permanent field. The strength of the main field at any time is called its normal 
value at that time. 

Although the rough and general description of the geomagnetic field in terms 
of its normal value as shown on magnetic charts is comparatively easy, the deri- 
vation of the normal values themselves is complex and difficult. ‘To obtain these 
normal values it is first necessary to remove from the individual magnetic obser- 
vations the contributions of extraneous fluctuations not a part of the secular 
variation. Often the extraneous fluctuations are large enough to affect seriously 
the estimates of secular change obtained from two or more observations made in 
different years at the same place. ‘These fluctuations may thus make difficult or 
impossible a reliable estimate of the normal values (since the secular change may 
not be determined accurately either with respect to magnitude or sign) at epochs 
other than those of observation. When the secular change has been corrected for 
fluctuations in field, so that it is known with accuracy, the magnetic observations 
are then readily reduced to the epoch desired for a magnetic chart. 

Of fundamental importance in charting the Earth’s magnetic field for a given 
epoch are charts showing rates of annual secular change, that is, isoporic charts. 
Too little attention has been given by various organizations responsible for the 
preparation of world isomagnetic maps to the importance of constructing isoporic 
charts. This importance arises because it provides the only feasible means of 
enhancing and extending the value measured at a station, say in 1920, to determine 
the value for the station at a later time, say 1945. 

The first comprehensive world isoporic charts were prepared by Fisk (1929 
and 1932) and by Fisk and Fleming (1931) at the Department of Terrestrial 
Magnetism for the epoch 1922-5 (specimen given in figure 5). These give the 
average annual secular change in the magnetic elements and their components, 
namely, declination, inclination, and the horizontal, vertical, total, northward, and 
eastward components of the field. 

Fisk’s isoporic charts, which have been widely used, were naturally tentative in 
a number of regions. ‘They did not take into account corrections of the observed 
data for geomagnetic fluctuations. They are also in several respects mutually 
inconsistent with the known nature of the geomagnetic field; and the singularities 
in field near the poles could not be taken into account because of the sparsity of 
data in high latitudes. 


The danger of using these charts today in mapping is that the pattern of the 
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isopors changes rapidly. For instance, in declination there is scarcely the slightest 
resemblance between the isopors of 1922 and those of 1942 in many parts of the 
world. Uncertainties of this kind account for relatively large discrepancies on 
some charts to actually observed values, In other words, it is necessary to know 
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Figure 5. Isoporic chart for magnetic declination in minutes per year, approximate epoch 1922. 
Prepared at the Department of Terrestrial Magnetism, Carnegie Institution of Washington, 
by H. W. Fisk, 1930 (slightly revised, October 1942). 


not only the pattern of the isopors at a given time but also the trend of change in 
pattern. ‘This trend can only be derived by mapping the isopors for several 
different epochs. This in turn requires that there be maintained adequate 
surveys at regular intervals and continuous observations at magnetic observatories. 


§3. GEOMAGNETIC SURVEYS 


The accuracy of magnetic observations depends mainly on the type and 
character of equipment, plan and technique of operation, and skill and care of the 
observer. ‘Thus measurements of a half-century ago are less accurate than those 
made in more recent years because of improved and superior instruments and 
techniques. An important illustration of this trend was the introduction of the 
earth-inductor for measuring magnetic dip to supplant the dip-circles of earlier 
generations, with substantial increase in accuracy. Loss of accuracy may result 
from deterioration of an instrument after years of use. The observational pro- 
gramme and the number of measurements required at a station vary with each 
data-gathering organization. Although procedures of measurements are well 
standardized on the whole, the results obtained by different observers show small 
differences when using the same instrument. 

There are many sources of error: for example, magnetic impurities in the 
metal parts of an instrument; changes with time in adopted standards, such as 
change in the moment of inertia of the oscillating magnet-system in measuring 
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horizontal intensity ; changes with time in deflection distances; wear in bearings 
and brushes; unspecified changes in instrumental constants not readily ascertained 
but detected by inter-comparisons with carefully maintainedstandardinstruments; 
poor levelling and influences of temperature-changes during measurements; 
infrequent inter-comparisons of magnetic instruments with reliable standards; 
‘and faulty computations of results—the last often very serious. Fortunately, in 
most surveys, errors from these sources are of minor character and-less important 
than the influences of extraneous fluctuations on data used in mapping the main 
field and its secular vaiation. 

The large portion of the Earth’s surface covered by the oceans makes the 
determination of accurate values at sea a major objective of the world-wide 
magnetic survey. It was not until 1905 that full realization of this objective had a 
definite beginning through systematic oceanic surveys then sponsored by the 
Carnegie Institution of Washington. 

As you are well aware, the first attempt to accomplish a magnetic survey at 
sea was the expedition of Halley between 1698 and 1700 (Chapman and Bartels, 
1940; Fleming, 1940). He made several voyages in the North and South 
Atlantic oceans, determining magnetic declination only. The results were 
embodied in his chart ‘‘ Lines of equal magnetic variation” of the Atlantic for the 
year 1700—the first isomagnetic chart. In the past 100 years extensive magnetic 
surveys have been made by many countries. Important undertakings following 
were in the Erebus, the Terror, and the Pagoda between 1840 and 1845, chiefly in 
southern waters; on these, all elements were observed utilizing the Fox dip-circle 
for measuring magnetic inclination and intensity. The frigate Novara obtained 
magnetic declinations while circumnavigating the globe in 1857-60. Obser- 
vations of the three magnetic elements were made over various oceans during the 
notable cruises of the Challenger in 1872—76 and of the Gazelle in 1874-76. More 
recently, observations were also made by the naval services of the various countries 
and by antarctic expeditions, notably the Discovery and the Gauss. The accom- 
panying charts (figures 6a, 66, 6c) show the tracks of the chief vessels on which 
magnetic observations were made from 1839 to 1916. 

All these observations were of varying degrees of accuracy set by available 
instruments and by the disturbing factors originating in the magnetic character of 
the vessels. ‘The distribution, both as regards position and epoch, was not such 
as to yield coordinated charts applying to definiteepochs. Therefore, in planning 
in 1904 to fill in the gaps in data for the remaining regions on land and sea for the 
world magnetic survey, the Carnegie Institution of Washington gave careful 
consideration to the oceanic survey. ‘The Institution’s earliest work was with the 
chartered vessel Galilee in the Pacific Ocean during 1905 to 1908. The Carnegie 
was designed in 1908 primarily for magnetic surveys and investigations, and the 
vessel and equipment were completed in 1909. Before the loss of the Carnegie at 
Apia in November, 1929, the ten cruises of that vessel and the Galilee had traversed 
some 362,000 nautical miles and from tens of thousands of individual observations 
had obtained values of declination at over 3800 stations and of both inclination and 
horizontal intensity at over 2300 stations. Figure 7a shows in a general way the 
scope of these activities from 1905, whilst figure 7b shows a chart of horizontal 
intensity obtained from the existing sources of information in 1940. 


Geomagnetic secular variations and surveys 22K 
¥ 

On the side of practical application, the increasing use of the oceans in the 

commerce of nations by sea and air makes the contin 


uation of the survey a matter of 
international concern and benefit. 
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Figure 6a. Tracks of chief vessels on which magnetic observations were made in the 
Indian Ocean, 1839-1916. 


continuation of the surface magnetic survey is the determination of the acceleration 
of secular variations; accumulated data show that they cannot be extrapolated 
reliably over periods as long as five years. A definite control is necessary for a 
number of epochs to facilitate the investigation of causes producing and governing 
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Atlantic Ocean, 1839-1916. 
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these progressive changes, which, it appears, would be favoured by accurate 
knowledge of their accelerations and distribution. The observations at repeat- 
stations on land have helped greatly, in conjunction with more detailed results of 
magnetic observatories, in delineating the complex patterns of secular change on’ 
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Figure 6c. Tracks of chief vessels on which magnetic observations were made in the 


Pacific Ocean, 1839-1916. 


land. By continuity these have served in some uncertain measure to bridge the 
gap across the oceans as well, when supported by the secular changes given by the 
coarser and less accurate ocean observations. It should be noted here that, though 
the cruises of the Galilee and Carnegie afforded usually only rough indications of 
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secular change in field, there were discovered great discrepancies in existing charts 
of the field itself. Great gaps in information remain in polar regions; for example, 


v2" 130° _140"_150°_ 160° 170° 180° 170° 160° 150° 140° 130° 120° 110" 100° 90° 80° 70° 
—— == i. a < —= = Be —— x x _ 


GEE 


CHELYUSKIN 2- 
1918-1919 7H 


..” POINT BARROW 
\932-1933 
a 


“A D Eyed aie Wy; a 


~ 60 
LEGEND~ SURVEY ON LAND < LEGEND - SURVEY ON SEA 
FIELD SIATIONS, MANY RECCCUPIED SEVERAL TOMES @ Z CRUSLS, CARNEG/E INSTITUTION OF WASHINGT OW 
= i GALEE (1905-1908) —— — — —— 
¢ MASNETIC AND ELECTRIC Ce Gren A beng ite eras x Y ry al 

a CARNEGIE INSTITUTION OF PASH (CLUETT (194) + Fy 

TUTTLE AMERICAR FeAPORARY COOPERATION, 6 70 18 ONTHS ie q dy i A 

(78°3475, 163° 56" Pemaotirwrannccomeetne = @ 2 GY CRUSE OF THE MAUD (918-1825) -- ~~~ =~ ~~ Ve G Ys Mp Grit Be 

1926-1930 AND 
1933-1935 21 LANO STATIONS COOPERATION WITH ayRo ANTAR 


CTC 
EXPEDITION I BETWEEN 76° T08I"S. AND I4S* TO 164" W. DLL, 2 
— o ig - . J o = = . == a . a5 z : as =a a —— . a € oF . . - 
130° 140° 150160" 170° 180 170° 160° 150° 140° 130° 120" 110° 100° 90° 80° 70° 60° so” 40° 30° 20° 10° 0° 10° 20° 3D% 40° 50° 60° 70° 80° 30° 100° 110° (20 


Figure 7a. Magnetic survey operations of Department of Terrestrial Magnetism, 
Carnegie Institution of Washington, from 1905. 
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Figure 7b. Isomagnetic chart of horizontal intensity of the earth’s magnetic force 
for the year 1940 (expressed in centimetre-gram-second units). 


the south polar regions have only about four badly distributed stations yielding 
inferior indication of secular change since 1900. 


It is gratifying that the British Admiralty designed and constructed the non- 
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magnetic vessel Research to continue the surface survey at sea and to fill in those 
serious gaps in the present data caused by the rapid change in the secular variation 
in certain regions of the Earth since the earlier surveys. But the task in the 
geophysical survey of the oceans is so great that other hydrographic services of 
maritime nations should also provide suitable surface or air vessels with equipment 
and personnel to take appropriate and coordinated share in the world survey. 


§4. ISOPORIC CHARTS AND GEOLOGIC FEATURES 


That irregularities of the distribution of the Earth’s magnetic field over conti- 
nental areas may bear definite relation to geologic features of the Earth’s crust has 
been recognized for along time. It is only in comparatively recent years that such 
magnetic data have been intensively considered as an aid in the determination of 
the character of geologic substructure. From accumulated observations it is 
possible to account for a considerable part of the Earth’s main field on the basis of a 
uniformly magnetized sphere. ‘The difference between such a uniform or normal 
magnetization and that actually existing serves to bring out the general irregular- 
ities as a residual field. The Earth’s outer crust or shell, perhaps 25 or more miles 
thick, is not homogeneous and, therefore, is not uniform in its magnetic behaviour. 
Thus there are many regions of local magnetic disturbance, such as are caused by 
magnetic ore-deposits, some so great as to give to local poles and other irregularities 
of intensity twofold or even threefold the normal value in such regions. Anomalies 
in Russia, in Germany, and in New South Wales are a few striking examples of 
such magnetic irregularities. ‘The great anomaly in South Africa has been 
apparently retained for over 100 million years by the huge Pilansberg system of 
Paleozoic volcanic dykes (figure 8); fourteen vertical dykes are found to be 
strongly magnetized in a direction opposite to the present-day magnetic field of 
the Earth (Piggot, 1938). In this general region of South Africa, the magnetic 
horizontal component of the Earth’s field has decreased 16°, of its present value 
in a 30-year period ! 

The anomalous magnetic features thus observed over continental areas must 
have their counterparts in the great oceanic basins (Fleming, 1938b). We 
may therefore look confidently to such features as an aid in the study of oceanic 
structure. Indeed, observations made on isolated islands and island groups 
almost invariably indicate abnormal magnetic conditions, and the magnetic 
survey of the oceans has revealed such conditions at sea. 

However, the application of magnetic-survey results to the investigation of 
crustal features of the ocean bottom depends upon the existence of a difference in 
substructural physical properties as well as upon the magnitude of this difference. 
Here enters the magnetic susceptibility of rocks, some of which are strongly 
magnetic, as iron-ores and igneous and metamorphic rocks, and some weakly 
_ magnetic, as most sedimentary rocks. Here, then, one must think of the observed 
magnetic anomalies as a kind of general reconnaissance method to indicate areas 
for other more intensive geophysical studies. However, the magnetic method of 
approach lacks depth-control, and, therefore, the size of geologic features must 
increase in proportion to depth if we are to detect such features. ‘Thus, unless the 
magnetic anomaly observed at sea is over shallow water or, if over deep water, is of 
enormous proportions, there may not be great hope of immediately deriving exact 
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information regarding particular areas of bottom substructure from distribution 
data as shown only by isomagnetic maps. 

The secular-variation change of the Earth’s magnetic field is a phenomenon 
which has large potential value in the study of the Earth’s crust. Herein is the 
great need of continuance of the accurate world surface magnetic survey, that we 
may determine the character of such changes with epoch and their accelerations at 
various places. 

Ever since the discovery of the change in the Earth’s magnetic field through the 
centuries, investigators have noted the apparent symmetry and regularity of 
secular variation when attention is confined to a single station or to a small part of 
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Figure 8. Magnetic anomaly of Pilansberg Dyke System, southern Transvaal, showing 
anomaly of vertical intensity expressed in gammas (after Gelletich). 


the Earth’s surface such as western Europe—only one-hundredth of the Earth’s 
surface. Considering, however, the actual observational data when the whole 
Earth is taken into account, secular variation is not so simple a phenomenon that 
it may be explained by periodic movements of the Earth’s magnetic poles and the 
like, but is a regional phenomenon connected with a large geological structure— 
oceans and continents—of the deeper layers of the Earth’s crust. 

The foci of rapid isoporic changes raise questions of broad geophysical 
significance. A large part of these changes may have origin in influences impressed 
by forces ceaselessly at work within the Earth. The rise and fall of the rate of 
secular change and the slow expansion and then the gradual retraction of the areas 
within which there have been excessive alterations in the magnetic elements are 
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significant of such changes. The distribution of isoporic foci (figures 9a, 95) is 
practically all in the hemisphere containing the great land-masses with the 
intervening Atlantic Ocean. Such foci as are found in the Pacific Ocean are of 
moderate intensity and not well defined. These relations to the surface structure 
of the Earth can be scarcely accidental; it is natural to suspect that there is a causal 


relation between crustal or subcrustal movements and conditions and these 
magnetic manifestations. 
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Figure 9a. Longitudinal distribution of proportion of annual change (4H/H) of horizontal 
intensity (after Fisk). (Lower curve represents average positive values of 4H/H in each 
lune between meridians, while upper curve represents numerical magnitude of average 
negative values ; thus shaded areas between curves are measures of excess of negative over 
positive annual change.) 


9° 60 


| 3 9c°w Yep 90° 180° ; eo°w 
q SCALE OF LONGITUDE 


eo oe OF LAND AREA 
a a 
c N 
PACIF(C_OCEAM Z 
ACTIVITY \ ———— <a 


NORTH AMERICA 
SOUTH AMERICA 


SCALE PER CENT OF LAND 


x 
9° 
w% 
~ 
x 
sey 
o 


EUROPE 
AFRICA 


AUSTRALIA 


Figure 9b. Variation with longitude of 4H/H (annual change averaged without regard to sign), 
of distribution of proportion of land and water areas, and of secular-variation activity 
approximately determined by density of distribution of isoporic lines (after Fisk). 


As an example, attention may be called to the apparent diminution of the 
intensity of the Earth’s magnetic field which is marked over oceanic areas, especially 
in the Southern Hemisphere (Fisk, 1932). ‘The interpretation of such data must 
be important in geophysical and geological research to advance understanding of 
Earth phenomena. 
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Thus the observed longitudinal distribution of magnetic secular variation 
agrees with that of land areas—as witness the moderate rates of annual change over 
the Pacific as compared with those over the Atlantic and adjoining continental 
areas. A definite control is necessary for a number of epochs to facilitate investi- 
gation of causes producing and governing these progressive changes. 

Further data bearing on correlation between the surface-distribution of 
secular-change activity promise conclusions concerning isoporic processes 
localized in the crustal layer. ‘Thus continued secular-variation surveys at sea 
should bring together seismic and magnetic methods of approach to crustal 
adjustments and possibly gravimetric work. 

Secular-variation information may be extended to remote times for which 
there are no observations. When lava cools and solidifies, following a volcanic 
outburst, it takes up a permanent magnetization dependent upon the orientation 
of the Earth’s magnetic field at the time. This, because of small capacity for 
magnetization in the Earth’s magnetic field after solidifying, may remain practically 


Figure 10. Apparatus for measuring magnetic moments and orientations of samples of 
varves and of ocean-bottom cores. 


constant. ‘Thus the direction of the originally acquired permanent magnetization 
can be determined by laboratory tests provided every detail of the orientation of 
the mass tested is carefully noted and marked when it is removed. 

Just as lavas cooling through the Curie point in the presence of a magnetic field 
become magnetized in the direction of the field, so sediments containing micro- 
scopic magnetic particles also assume the direction of the prevailing magnetic field. 
This is truly “ fossil” magnetization, for the magnetic field of the epoch of deposi- 
tion is embedded in the sediment as surely as the organic remains of that epoch. 
Whether this fossil magnetism remains as a connate property of the rock or whether 
it is destroyed by subsequent events, even as organic fossils are sometimes 
destroyed, depends upon the post-natal experience of the sediment. In so far as 
it remains, it may supply an important record of the Earth’s magnetism and an 
important clue to the history of the rock. 
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‘fy Lhe measurement of the Earth’s magnetism at great depths in the sea may now 
be consummated by magnetic examination of the core-samples taken with 
Dr. Piggot’s apparatus described in 1938. Inorder that the fossil magnetization of 
the specimens may not be disturbed, a method of measurement has been developed 
for tests of high sensitivity and accuracy without submitting the specimens to 
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Figure 11. Direction of magnetization of sediments in core No. 3 
(curve smoothed by formula (a+ 2b-+3¢-+2d-- e)/9). 


artificial magnetic fields. This equipment (figure 10) permits detection of a 
magnetic moment as minute as that of a good steel magnet one millimetre long 
and one-hundredth of a millimetre in diameter. Further refinements will 
increase many-fold the sensitivity of this method. 

Tests on one core show distinct differences in the direction of magnetization at 
various depths (figure 11). If it is found that these differences are the result of 
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differences in the Earth’s field at the time the sediments were deposited, we 
have in the core-samples a means of extending the past history of geomagnetism 
and its secular changes. Interrelating individual cores by similarities in their 
magnetism would determine more accurately than surface observations the 
extent and character of the regional magnetic anomalies of the ocean bottom. 
With improved knowledge of magnetic anomalies, the question of interrelations 
with gravimetric anomalies and deep-focus earthquakes may be answered. It 
is frequently noted that gravimetric and magnetic anomalies occur in the same 
regions. Apparently areas of deep-focus earthquakes are confined to the great 
regions of positive anomalies of magnetic vertical intensity around the Pacific 
Ocean. Vening Meinesz has found an excess of gravity in the oceans and a 
sudden decrease at the continental borders, and has demonstrated that these 
deviations of isostasy cannot be situated in the Earth’s thin crust but must extend 
far into the plastic layers below. There must result mass flow in the Earth’s 
interior, downward below the oceans and upward below the continents. The 
charts of isoporic foci and their motions with time certainly lend strength to the idea, 
indicating the interior of the Earth to be more mobile than the external layers, not 
only as a whole but regionally. May not, then, the anomalies of gravity and of geo- 
magnetism and the deep earthquakes have a common origin in mass flows in the 
inner Earth? ‘Thus the facts of the geophysical observer of geomagnetism in the 
field and in the laboratory may contribute much to a complete picture of the 
Earth’s crust and particularly of that great portion which is the ocean bottom. 


$5; RECENT DERIVATIONS OF SECULARS CHANGE 


The Department of ‘Terrestrial Magnetism has compiled isoporic charts of the 
geomagnetic field for the four epochs of 1912-5, 1922-5, 1932-5, and 1942-5 for 
declination, inclination, and the northward, eastward and total components of 
intensity. ‘hese were compiled preparatory to the construction of new world 
isomagnetic charts of these elements in seventeen sections for the epoch 1945-0. 

All available repeat-station observations have been used in compiling secular 
change for each of the seven elements and components. All have been corrected 
for the influence of extraneous geomagnetic variations, including those related to 
the sunspot cycle (not properly considered as secular change), post-perturbation, 
annual variation, and the quiet-day solar daily variation. Graphs of the observed 
values so reduced to “normal” were made for each station. Specimen graphs for 
the station Suez and the Alibag Observatory are reproduced in figures 12 (a, b, c) 
and 13(aand 6). ‘The values as actually observed are plotted in order to show the 
order of correction required to obtain reduced ‘“‘ normal” values. The corrections 
at Suez are largest for the horizontal and northward components, due to the 
considerable latitude-effect of the geomagnetic variation with sunspot-cycle and 
of post-perturbation. The corrections are of greatest percentage benefit when 
the secular changes are small. ‘The same considerations apply also in the case 
of observatories such as Alibag; here the corrections (figures 14.@ and 144), as for 
many series of data, have been of material assistence in improving estimates of 
secular change in intensity near 1910. 

Graphs of the foregoing kind were scaled at epochs 1912-5, 1922-5, 1932-5, 
and 1942-5, and plotted on world maps. A sample of these is given in figures 
15 (a, 6, c) and 16 (a, 6, c), showing the isopors in declination and horizontal 
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Figure 13 a. Specimen graphs for values of D, H and I from 1905, Suez. 
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Figure 13 6. Specimen graphs for values of X, Y, Z and F from 1905, Suez. 
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Figure 146. Specimen graphs for values of X, Y, Zand F from 1905 for Alibag. 
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a. Isoporic chart for declination (lines of equal annua 
per year), epoch 1912°5. 


PAeSer. 
OO 


YS 


Ny 
Vi 
NAY 
we 
TA 


NTH: 
[ | f Y VE eae. 


Bl) 
Ww) 


eg 


aw 
f PX 
z MG \ 


Figure 155. Isoporic chart for declination (lines of equal annual change in minutes of arc 
per year), North Polar region, epoch 1912°5. 
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Figure 15 c. Isoporic chart for declination (lines of equal annual change in minutes of arc 
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intensity. for epoch 1912-5. “These were-first drawn independently of one another 
and then adjusted to mutual consistency such that the vertical component of the 
curl of the field was nearly zero at the corner of each of the 10°-tessera. The 
adjustments were further continued to ensure vanishing of line-integrals of the 
northward component along meridians and of the eastward component along 
parailels of latitude. | | 
The isopors near the magnetic poles are based on a theoretical fit of available 
data, based on a power-series expansion for the potential due to Chapman (1942), 
The theoretical fit of secular changes near the dip-poles was also used to estimate 


Hy 
Saal we atias my e 


AME 
ved aa 


SEEPS KOK) 
SEEEPERE SS 


Isoporic chart for horizontal intensity (lines of equal annual change in 
gammas per year), North Polar region, epoch 1912:5. 


the positions of the north and south magnetic poles given in table 1, These 
tentative estimates assume t 


he revised location of the north magnetic pole by 
Wasserfall (1939) based on Amundsen’s observations of 1904-5 (Graarud and 
Russeltvedt, 1925) and of the recent flight of the Aries. The spacing between 
successive positions of the P | 


he pole at the various intervals of time: are adjusted in 
Proportion to the velocities of the pole estimated from secular change. The 
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Figure 16 c. Isoporic chart for horizontal intensity (lines of equal annual change in 
gammas per year), South Polar region, epoch 1912°5. 


Table 1. Coordinates of north and south magnetic poles. 


Based on secular-change charts for four 
epochs and positions observed in 1904:5 for 
north pole and 1912-5 for south pole 


Based on  isomag- 
netic charts for four 
epochs and positions 
observed in 1904:5 


pat and 1945 
North East South East North East 
lati- longi- lati- longi- lati- longi- 
tude tude tude tude tude tude 
1904°-5 70°5* 263 -°5* ary se 70:5* 263-5* 
1912-5 70:9 263-2 71-2+ 150-5 MAES 262:7 
1922°5 71:4 262°8 70:2 149-2 72:6 261°8 
1932°5 71-9 262°4 69-0 148-1 74-0 260-0 
1942°5 72°6 262°1 68°3 146-2 75:4 258:2 
1945:0 72°8 262:0 68-2 145-4 76:05 258-0t 


* After Wasserfall (1939), based on observations by Amundsen. 
+ After Farr (1944), based on observations by Mawson. 
{ Preliminary values communicated verbally as observed on Aries. 
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The isoporic charts show agreement with Chapman’s rule that the number of 
nodes equals the number of foci minus 2 (Chapman, 1942). For declination there 
are eight nodes and ten foci with no ray-poles. For horizontal intensity there are 
ten nodes, ten foci, and two ray-poles, the latter counting as foci in the application 
of the rule. New features also depicted for the first time are the distribution of 
lines near the magnetic poles. 

The ten-year successive positions of positive and negative foci in declination 
and horizontal intensity for the four epochs ten years apart are plotted in figure Nii 
The most notable feature is the large and systematic departures in position per 
ten-year period. ables 2 and 3 list the locations of foci, their magnitudes, and 


1912.5 1222. 
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Figure 17. Successive approximate ten-yearly positions and magnitudes of isoporic foci; 
declination in minutes of arc and horizontal intensity in gammas ; epochs 1912°5, 1922°5; 
1932-5 and 19425. 


positions. In the Eastern Hemisphere the foci of the secular change in horizontal 
intensity show in general circulation in position to the north and to the west. This 
finding is in agreement with the general results for the position of the magnetic 
equator indicated in figure 3. In the Western Hemisphere there seems some 
indication of motion to the south. The foci for secular change in declination, on 
the other hand, seem to show motion in earlier years which is in certain Be: 
towards the east, there being less regularity than that noted for the secular change 
in horizontal intensity. 

‘These new isoporic charts offer several new and interesting possibilities. 
Although constructed independently of one another, the differences of the chart 
epoch by epoch, reveal surprisingly simple consistent patterns. These on 
patterns represent the time-rate of change of secular variation per ten-year 
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interval. It should prove interesting to try and relate these new values to other 
geophysical material. The consistency of the differences found are encouraging 
in another way. ‘They may indicate that the isoporic charts of the northward and 
eastward components may be successfully subjected, for the first time, to spherical 
harmonic analysis. It is planned to undertake such a trial analysis, a check being 
provided by comparing the separate results for the two components. ifstie 
analysis is satisfactory, necessary and desirable improvements in the isoporic 
charts for vertical component can be made, for instance, in polar regions and over 
the oceans. he extrapolation of the secular-change field or its equivalent in 
terms of a thin current-distribution at various depths within the Earth is also 
contemplated. 

An interesting hint of the results of this extrapolation can perhaps be derived 
by noting the considerable complexity of the isoporic charts in horizontal intensity 
and declination (no attention should be paid to complexities near the magnetic 
poles). For this purpose we should preferably have the vector combination of the 
values of secular change in declination and horizontal intensity. ‘The complexity 
in pattern must increase rapidly downwards in the Earth, and it is highly probable 
therefore that this extrapolation would attain unreasonable complexity at modest 
depths, only a moderately small fraction of the distance to the Earth’s centre. 
We thus restrict the probable site of the origin of secular change to smaller spacial 
dimensions than hitherto, perhaps within the mantle of the crust of the Earth. 


§6. FUTURE GEOMAGNETIC SURVEYS 


The present slow and costly methods of geomagnetic surveys seem likely to be 
supplemented by new techniques and methods, the possibilities of which are now 
becoming gradually apparent. « Clearly, the old procedure of measuring the 
Earth’s field at different times in different places in one region during one year or 
years and in another at a different epoch is inadequate to all the needs. What is 
required is a description of the geomagnetic field during a given year based on 
measurements made during that year at an orderly set of points sufficiently close 
together and spaced in a manner nearly independent of topography, areas of land 
and sea, and climate. The prophecy of Chapman in the first Charles Chree 
Address (1941) that such measurements might well be made from planes promises 
early realization. A project of this kind is now feasible instrumentally and 
affords attractive post-war possibilities in application. 

The available data are already of necessity squeezed rather hard to obtain 
useful isomagnetic charts for the epoch 1945 in some regions. 

As stated above, charts of first accuracy could not be prepared without extensive 
and uncertain extrapolations for secular variations; it was impossible to meet 
military specifications for accuracy on charts on which our magnetic-mine warfare 
and degaussing were to be based. Mathematical adjustments were resorted to in 
making up for the paucity of observations. Believing that a non-potential field 
does not exist, maps were mutually adjusted so that the curl of the field vanished. 
It is probable that these adjustments give a more accurate representation of the 
field than could otherwise be obtained, but no amount of mathematical adjustments 
can ever compensate for absence of observations. 

We have sought to improve our knowledge by increasing the accuracy of field 
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observations, and this step has been important in accurately determining secular 
changes. However accurately the magnetic elements are determined at a given 
point, absence of observations at intervening points leaves our information still 
inadequate for surface isomagnetic maps unless we can establish almost an infinite 
number of stations. 

This condition may now be remedied by magnetic observations withinstruments 
airborne above the surface of the Earth, and thus farther away from magnetic 
anomalies which depreciate the accuracy of surface observations. Travelling 
at speeds which a few decades ago seemed fantastic, continuous observations 
may now be obtained over a large area in a single day, and ina relatively short time 
a survey of the normal geomagnetic field can be made. The accuracy of such 
techniques, thanks to the achievements during the past war, can be kept more than 
adequate for both practical and theoretical purposes. Electronic circuits, 
monitored by standard cells, can maintain standard electric currents in flight, so 
that magnetic profiles, accurate to a few tenths of gammas or better, in absolute 
value in any component can be obtained. The following limitations in the present 
airborne gear must be overcome:—(a) Stable operation must be provided over 
extended periods of time of measurements of intensity: this is synonymous with 
negligible base-line and scale drifts in the measuring circuits. (6) Modifications 
are necessary to permit measurement of the angle of declination, that is, the angle 
between the horizontal projection of the detector coil at any instant and true north 
or the true heading of the aircraft. (c) Observations of the angle of dip will involve 
the construction of a suitable gyro mechanism to establish an accurate horizontal 
plane; among problems in this connection are the elimination of centrifugal 
pseudo-gravity effects due to plane manoeuvres or bumpy air. 

The advantages of an effective airborne surveying instrument are great, and 
include:—(a) Ability to cover large areas quickly and continuously. (6) Elimination 
of purely surface anomalies, which are of no significance in the preparation of charts. 
of the normal field; thus, by flying at altitudes of, say, 20 000 feet, smoothed values 
of the geomagnetic elements are obtained directly, without the necessity of mathe- 
matical treatment or judgment onthe part of the person processing the data of many 
agencies and of various degrees of accuracy. (c) Ability to define uniquely the 
magnetic field by flights in two or more horizontal planes. 

Perfection of a suitable airborne instrument for magnetic surveying will require 
both considerable time and funds to overcome these limitations. Nevertheless, 
the end results possible would justify such expenditure of both time and funds. 

Azimuths for declination can be obtained by automatic or manual observations. 
of the Sun, smoothed by gyroscopic stabilization. An example of what can be 
done with airborne methods is supplied by the recent flights (August 1945) of the 
Aries under sponsorship of the Imperial School of Navigation on the initiative of the 
Astronomer Royal. In two flights of less than 24 hours over entirely different 
regions, a new position of the magnetic pole was determined, probably more 
accurately than ever before. This achievement may be compared with the 
tedious and dangerous expeditions of Ross and Amundsen to accomplish the same 
purpose many years ago. The instruments used were magnetic compasses of 
various types, including gyrostabilized compasses, and also a three-element, 
saturated-core type of magnetometer. The last, of a relatively low order of 
accuracy, was not designed for the flights; it had to be held level manually during 
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an observation. Furthermore, only one element could be observed at a time, and 
there was no means for continuously recording—thus truly average values of the 
observations obtained seem consistent. This pioneer work demonstrates 
what may be expected when magneticians take full advantage of the technical 
progress of the last five years. 


§7. CONCLUSION 


During the war, as by-products of military demands, many improved techniques, 
applicable not only for surveys and experiments but also for mass reductions, 
analyses, and correlations by machine methods, have been developed. Peace-time 
applications of these must serve both continuation of surface and initiation of 
aerial magnetic surveys by our several Governments, as well as interpretative 
investigations of the vast existing stock of data and that potentially possible in 
much greater degree in the near future. It is to be hoped that a long period of 
amity and cordial relations between all nations is now on the horizon. May it give 
firm foundation in forwarding and coordinating large activities in Earth physics 
based on that truly scientific principle of entire lack of desire for personal or 
national aggrandizement. 

Finally, such progress requires the combined efforts of the physicist, the 
mathematician, the geologist, the astrophysicist, and the engineer to enable the 
geophysicist to solve the riddles of the Earth’s main magnetic field and of its 
secular variations. 
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Figure 3. The original cavity magnetron detached from pumping system. 


Figure 4. The CV 64, a 9-10-cm. echelon-strapped magnetron much used by the Services. 
Peak power output ~35 kw. 
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Figure 5. Megawatt magnetron with end-plate removed to show strap construction. (Scale in cm.) 


Figure 6. Megawatt megnetron; second view, showing wave-guide output system. (Scale in cm.) 
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THE, CAVITY MAGNETRON* 


DY oe ho UND AL keno. 
University of St. Andrews 


An Address delivered 12 December 1945, with demonstrations, on the occasion of the 
presentation of the 22nd Duddell Medal to Prof. Randall 


MS. received 17 February 1946 


on me I wish to offer my most grateful thanks. I have very little to guide 

me by way of precedent on this occasion. Either by accident or design the 
remarks of previous recipients of the Medal have but rarely been reproduced 
in extenso. I must therefore chart my own course. Before I refer to the work 
you have singled out for special mention I would like to say this: 

Not only do I owe a great deal to the laboratories in which I have worked 
and to their Directors, but I have been specially fortunate in the young people 
who have collaborated with me. If I had not been trained in Manchester under 
W. L. Bragg, I might not have developed the interest in crystals which was so 
valuable as an approach to the problems of luminescence. If I had not spent 
some years in the Wembley laboratories my approach to valve problems would 
have been inadequate when I came to think of short-wave generators. And, 
indeed, if I had not been in the University of Birmingham when Oliphant started 
his great drive on centimetre-wave work, I should probably have had nothing 
at all to do with such instruments. I say these things because I believe that at 
any given time one’s past experience affects one’s approach to a problem, and 
determines more than many people think the way in which a problem ts solved. 

The luminescence work, which started a good deal earlier than that on high- 
frequency generators, was of two kinds: technical and fundamental. I might 
mention on the technical side the discovery, by McKeag and myself (1936), 
of zinc beryllium silicate activated by manganese as a phosphor suitable for 
excitation by a low-pressure mercury discharge. This phosphor is now a com- 
ponent in many of the mixtures used in the low-voltage fluorescent lamps (see 
Jenkins, 1945). I think, also, we were among the first to elucidate the technical 
processes required to obtain a number of high-quality silicate and phosphate 
phosphors. 

On the fundamental side I should like to mention the work of Randall and 
Wilkins (1946) on the mechanism of afterglow or phosphorescence in different 


[' accepting the great honour which you and your Council have bestowed 


* This short title has been chosen since it covers the greater part of the subject-matter of the 
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types of material. In a series of recent papers it is shown how the characteristic 
thermoluminescence is related to electron traps in the solid; further, that the 
temperature at which light is emitted is a measure of the depth of the trap below 
the conduction band in the solid. From a knowledge of the relative numbers 
of electrons in traps of different depths it was possible to calculate the type of 
decay law obeyed over a narrow temperature range. Wilkins and Garlick (1945) 
have followed this work with a paper in which many earlier misconceptions 
concerning bimolecular decay laws have been cleared away. here is one 
important problem in this field which someone may wish to investigate: the 
nature of electron traps in particular cases. No direct attack is possible for 
obvious reasons. The detailed study of phosphors with more than one activator— 
rare-earth ions would probably be suitable—will probably give information, 
not only with regard to the luminescence centre, but about the trap, and optical 
and thermal excitation energies. 

The work in which I suspect you are more interested, however, is the discovery 
by Dr. Boot and myself of the resonator magnetron in the winter of 1939-40. 
Later, in 1941, Dr. Sayers made an invaluable fundamental contribution. As the 
magnetron was immediately successful as a short-wave generator it was taken 
up for the Government by the large electrical firms, and I should like to pay 
tribute to the work which was carried out by the staffs of the B. T.-H. andthe G.E.C. 
in this connexion. 

As the story has not really been told of how this work arose I should like to 
spend a little time on the origins of the magnetron as we now know it, and to follow 
this with a brief reference to later experimental developments. At the close 
of my talk I hope to show a powerful magnetron in operation, and for generous 
help in this I wish to thank the Admiralty Signal Establishment. 

In the early autumn of 1939, with others whom Professor Oliphant had 
collected together for this purpose, I began to spend some time on problems 
connected with centimetre waves, as many other physicists were doing elsewhere. 
The chief problem was to produce a high-power centimetre wave oscillator. 
The centimetric nature of the wave was required for direction of beam, and the 
power for long range. At first Boot and myself spent a few weeks studying the 
Barkhausen-Kurz tube as a detector device; the greater number of workers inthe 
laboratory, however, were concerned with the klystron (Varian and Varian, 1939), 
both as an oscillator and as an amplifier, and we were naturally interested in the 
whole field. ‘The klystron was outstandingly important in that it used for the 
first time closed, or essentially closed, resonators, fitted with grids through which 
an electron beam passed, for the production of high-frequency power. Such 
resonators had been considered by Rayleigh in 1897, and the papers of 
Hansen (1938, 1939) showed how enclosures of this kind lend themselves to certain 
modes of electrical vibrations. On any given mode the frequency is largely, 
if not entirely, dependent on the dimensions of the resonator. When such 
resonators are constructed of copper, three important features are evident :— 
(i) low h.f. losses; (ii) wave-length stability; (iii) potential capability of large 
heat dissipation. 

To be set against this, however, is the serious difficulty of getting sufficient 
power into the electron beam, which is necessarily of small cross-sectional area. 
We began, therefore, to think of other possibilities. Was there any means 


of using the desirable features of cavity resonators, and at the same time avoiding 
the limitations of the electron beam ? 
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Conventional designs of magnetrons were equally notable in failing to meet 
the desirable features, and by November 1939 we had decided to try to devise 
a type of magnetron which did so. The Hansen papers already referred to 
were by this time available in the laboratory, but it was obvious that the types 
of resonator considered—hollow spheres, cubes, etc.—could not be associated 
with the cylindrical anode and cathode of the magnetron. The only other types 
of short-wave resonator circuit with which we were familiar were Hertz’s original 
loop-wire resonator and a short-circuited quarter-wave line. A three-dimensional 
version of the first of these gives a cylinder with a slot down one side, and the 
second becomes a parallel-sided slot A/4 deep (see figure 1). Either design was 
well suited to the arrangement of a number of segments with cylindrical symmetry 
round a central cathode, and it is also clear that a combined anode and resonator 
system could be machined from solid copper to allow of good heat transfer and, 
it was hoped, high-power dissipation. 

No theoretical calculations of the wave-length to be expected from a resonator 
of the cylinder-slot type had been made, but Hertz had shown, and Macdonald 
(1902) had confirmed theoretically, that the wave-length to be expected from 
a simple Hertzian dipole was A=7-94d, where d is the diameter of the ring. 

Six resonators, 1-2 cm. in diameter, with slots 0-1 x 0-1 cm., were employed, 
and the length of the resonator was 4 cm.; it was hoped that the wave-length 
of this system symmetrically surrounding an anode-hole of the same diameter 
would be approximately 10 cm. Preliminary calculations had shown that this 
design might be expected to operate reasonably well with an applied potential 
of 16,000 volts and a magnetic field of approximately 1000 oersteds, so that 
provision of these essentials was a necessary preliminary. In order to save time 
a D.C. power-supply system was built up, and not a modulated supply, and a large 
electromagnet was used to give the required field parallel to the cylindrical axis 
of the valve. Figure 2 shows a view of the first anode block, figure 3 the first 
resonator magnetron, continuously operated on the pump with glass-metal 
junctions closed with sealing wax. For various reasons the first trials were delayed, 
and it was not until 21 February 1940 that the valve was tested. It worked 
immediately. By various rather crude means it was shown that the output 
was about 400 w. ‘The wave-length was measured by means of Lecher wires 
and shown to be 9-8 cm. After this, things moved rapidly, and I will not 
enumeiate the various technical stages of development, which included the 
making of sealed-off valves by ourselves and by G.E.C., Wembley, incorporating 
the famous gold-seal device as a means of fixing the end-plates to the valve. 
We also received great help on the technical side from S. M. Duke, formerly 
of the Wembley laboratories. Later, Professor Oliphant organized a small 
production unit for experimental magnetrons, and during its life over a thousand 
experimental valves were produced. 

There have been many stages in the development of the magnetron since 
1940, and it is obviously impossible to go into these in detail. It was, of course, 
important to see to what extent power could be increased. ‘The Wembley and 
Birmingham laboratories simultaneously introduced large oxide cathodes, 
although both laboratories were aware that the French had been experimenting 
on these lines with the older types of split-anode magnetron. During the 
summer of 1940 we were able to show that the cavity-resonator principle was of 
general application; any fears that might have existed concerning the applicability 
to larger numbers of resonators and shorter wave-lengths proved groundless. 
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During July 1940 experiments on a 5-cm. valve were carried out, and a small 
amount of power was also obtained at 1:9 cm. using 30 slot resonators. A little 
later, September 1940, a scaled-down version of the original 10-cm. valve 
(6 resonators) was operated at 3 cm. 

In the early design work we were aided by some elementary considerations. 
For electrons passing close to the anode without grazing it the operating voltage 
turned out to be 


V=2X107 x 72D2/N2242 volts, 9 nn ewe (1) 


where D is the anode-hole diameter, N the number of resonators, and ¢ the phase 
difference between adjacent resonator tongues. 
The corresponding expression for the magnetic field is 


H = 7-24 x 104/NX oersteds.| 5 ee eee (2) 


These expressions were undoubtedly based on quite inadequate theory; 
but as they were found extremely useful in the design of new valves they deserve 
a small place in the record of developments. It was later shown experimentally 
at Birmingham that the resonator magnetron will operate with higher efficiency 
at higher magnetic fields and voltages than those derivable from equations (1) 
and (2). The first magnetron had an efficiency in the neighbourhood of 10%, 
and the pulsed power output of some of the early valves was approximately 
50kw. Under the new operating conditions efficiencies of more than twice the 
above figure were found. Later (1942) it became clear that the high operating 
conditions were simply an approach to the 7-mode. 

The mention of power and operating conditions brings up again a feature 
with which we were concerned inthe earliest stages—that of high current or peak- 
current. From oxide-coated cathodes it is possible to get peak-currents of the 
order of 5 amp./cm? of cathode surface, and from this fact it is seen that the limit 
of moderate beam current which obtains in the klystron is overcome. 

Any discussion of cathodes would be incomplete without mention of secondary 
emission. From various pieces of evidence it became clear that electrons striking 
the cathode with finite velocity on completion of an orbit produced secondary 
electrons. ‘This led to an interesting experiment in which a magnetron was made 
with a simple aluminium cathode entirely unheated. A small quantity of hydrogen 
was admitted to the valve from a palladium tube and was sufficient to start the 
build-up of current, and such a magnetron worked normally apart from eventual 
clean-up of the gas. ‘This experiment showed for the first time the possibility 
of cold or even water-cooled cathode magnetrons, and foretold the ultimate, 
but still unrealized, development of such valves for c.w. purposes. 

The next stage of development in the magnetron has come to be known 
as “strapping”. In this development I have no personal share, but the story 
of the magnetron would be grossly inadequate without it. Essentially this is the 
work of Dr. Sayers at Birmingham. 

So far we have considered the magnetron as radiating only one frequency, 
and, indeed, this can be so. The resonator system has an infinite number of 
resonant modes with wave-lengths extending from some maximum value to lower 
and still lower values. The symmetry of the system causes the modes to be 
bunched together into groups of different frequencies, and each group by analogy 
with ordinary spectroscopy can be called a multiplet. In general each multiplet 
will contain N—1 modes if there are N resonators. Most of the modes are 
degenerate, ie. more than one mode is found for a given wave-length. All 
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resonant wave-lengths but one (if N is even), or all (if N is odd) are degenerate. 
Thus in a 10-hole magnetron there are five possible wave-lengths in a multiplet 
[1+(N/2—1)]; one of these corresponds to a single mode, and each of the others 
to two modes, giving a total of 1+2(N/2—1)=N-—1 modes or 1+2x4=9. 
The fundamental fact of the existence of different modes in a magnetron resonator 
system showed itself in the phenomenon known as mode-jumping. For a given 
set of operating conditions a valve could suddenly jump from one of the possible 
modes to another nearby one. This had obvious disadvantages, and I should 
like to describe briefly how this was overcome by Sayers at Birmingham. 

By connecting alternate segments together by short copper wires the circuit 
constants of the resonator system are altered, and the general effect can briefly 
be described as a greater separation of the modes of the multiplet. In nine- 
or ten-centimetre magnetrons this enables what is called the z-mode of operation, 
with a phase difference of 7 between adjacent segments, to be achieved. 
Examination of the wave pattern inside a cold block in these circumstances 
shows that the standing wave pattern is very regular in amplitude. The patterns 
for other modes are not suppressed but are irregular and often only weakly 
excited. ‘The CV64, a 9-10-cm. strapped magnetron designed at Birmingham 
and greatly used by the Services, is illustrated in figure 4. The possibility of 
producing magnetrons having efficiencies of about 50°% caused most of the 
British magnetron effort to be directed towards the design of strapped valves 
for particular Service applications at both 10 cm. and 3 cm. For the former 
wave-length no great difficulties were experienced. At first, considerable diffi- 
culties existed at 3 cm., but these cannot be detailed here; they were eventually 
overcome by adopting the American technique of scaling down a 10-cm. valve. 

The maximum peak power now obtainable from cavity magnetrons in the 
9-10-cm. region is of the order of a megawatt. ‘T'wo views of such a valve are 
shown in figures 5 and 6. ‘This is a great advance on the position existing in 
1940, and is an example of what can be done by cooperation between the various 
organizations (universities, firms, establishments) responsible. 

All resonator magnetrons operate above cut-off, that is to say, at a value of 
magnetic field for which there would be no anode current under steady conditions. 
There is a space-charge between anode and cathode, and the steady states without 
any high-frequency field on the anode have been analysed in detail by Hartree 


and Stoner. 
The theoretical approach to the problem of the interaction of the electron 


stream and the high-frequency wave which travels round the anode is extremely 
complicated; that many of the problems concerning the operation of cavity 
magnetrons have been solved is in large part due to the efforts of Hartree, Slater, 
Stoner and their collaborators. Perhaps the most important early theoretical 
work was carried out by the Slater group. First, calculations were made utilizing 
field theory to solve the problem of the possible number of frequencies of a 
multiple resonator system. Secondly, the significance of the space-charge in 
the magnetron was realized and investigations based on the assumption of a linear 
field distribution between anode and cathode were made. From a knowledge 
of known operating conditions it was found possible to calculate the mode of 
operation of a given valve. 

Later theoretical developments were in large part due to English workers, 
chiefly Hartree and Stoner. The team working under Hartree considered the 
problem of magnetron space charge more thoroughly and made calculations 
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based on a field distribution between anode and cathode which was a good deal 
nearer reality. This work led to the formulation of the now well-known threshold 
criterion. It was shown that exact synchronism between the outer edge of the 
space-charge cloud and the rotating wave was not required, but that operation 
of the valve, once started, would continue over a range of angular velocities of 
the cloud. The lower limit of angular velocity depended on the amplitude 
of oscillation ; it was possible to find a valve for which electrons would reach | 
the anode for vanishingly small segment amplitudes for any practical mode- 
number 7. ‘Thus an anode voltage was determined below which no current 
could reach the anode for a given magnetic field, and this voltage has been 
called the threshold voltage. The threshold relation may be expressed in 
equation form by 
V7 2H Oa ees 
10100a2 n * 21400 ~ 72? 

where Vis the threshold voltage in kilovolts, and n is the effective nm number 
of the mode. 

By the middle of 1942 the calculations, which have been very briefly outlined, 
had given considerable insight into the operation of valves in the 9-cm. to 12-cm. 
wave-length region. Large discrepancies, however, still existed between the 
calculated and observed operating conditions for 3-cm. magnetrons. In May 
1942, Sayers and Sixsmith, by phase measurements on 3-cm. unstrapped valves 
with 12 to 16 segments, showed that such valves were not operating in the 7-mode ; 
furthermore, it was shown that the direction of the rotating anode-potential wave 
was opposite to that of the rotating space-charge cloud. ‘This meant that a 
lower space-charge angular velocity was required for proper interaction for 
modes other than the 7-mode; the lower space-charge velocities corresponded 
much more closely with the observed operating conditions. The Hartree group 
worked ou the required modification to the earlier resonance criterion. Under 
the new conditions it was shown that for electrons to reach the anode at magnetic 
fields well above cut-off value, and for vanishingly small high-frequency amplitude 
a certain degree of synchronism is required between the rotating space charge 
and one of the spatial components of the rotating electric field in the anode- 
cathode space. 

The theoretical work on the magnetron which has been referred to has not 
yet been published, and it is extremely difficult to cover it adequately in such a 
short space. It is to be hoped that those directly concerned will take steps to 
give a comprehensive account in the near future. 
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ABSTRACT. This paper describes the principles and the construction of an apparatus 
for the precise measurement of the acoustic impedance, at normal incidence, of sound- 
absorbing materials by a stationary-wave method for frequencies in the range 100 to 5000 c/s. 
From measurements made on a sample of material 13 inches in diameter, backed by a 
substantially rigid wall, the magnitudes of the resistive and reactive components of the 
impedance may be calculated with an accuracy of about one per cent. 

The theory of the method is discussed with particular reference to the influence thereon 
of the attenuation of sound associated with the walls of the tube in which the standing 
wave is formed, and it is shown that in addition to the correction which must be provided 
to the elementary expression for the “‘ standing-wave ratio ’’, the finite attenuation in the 
tube leads to an additional correction to the expression for the distance from the sample 
of the first minimum of pressure. Attention is drawn to convenient methods which 
facilitate calculation of the results, and typical results of measurements of the impedance 
of a sample of porous sound-absorbing material are shown. 


SP INDRODUC TION 


N recent years interest has been steadily growing in the recognition of the 

acoustic impedance, rather than the absorption cofficient, as the chief 

characteristic of the performance of acoustical materials. ‘he advantages 
which follow from a knowledge of the impedance are twofold. Firstly, modern 
theories relating to the absorption of sound both in enclosures (Morse, 1939 a) 
and in ducts (Morse, 1939b) are formulated in terms of the impedance of the walls 
for normal incidence. Secondly, a knowledge of the separate resistive and 
reactive components of impedance of an absorbing material simplifies the task 
of relating the deficiencies in the absorbing power to the physical structure 
of the material. 

Various methods for the measurement of acoustic impedance have been 
devised, and are referred to in a recent paper by Beranek (1940). There is a 
paucity of published results by such methods, and until the publication of 
Beranek’s paper, no substantial accumulation of results of adequate accuracy 
had appeared. Beranek attributes this condition to the difficulties inherent 
n the methods of measurement previously in use, and proceeds to describe an 
apparatus of his own design which he considers to be capable of producing 
results of the required accuracy. 

In the present paper, a modification of the apparatus originally used by 
H. O. Taylor is described. This apparatus is capable of an accuracy of the same 
order as that of Beranek and possesses, in the author’s opinion, advantages of 
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simplicity in design and in theory of operation. Measurements of acoustic 
impedance can be made with this apparatus at frequencies in the range 100 
to 5000 c/s. 


§2. THE METHOD AND ITS JHEORY OF OPERATION 


The principle of the method is the same as that used by Paris, which in turn 
.s based on an extension of that devised by Tuma, and later used by H. O. Taylor 
for the measurement of the absorption’ coefficient. The sample of material 
closes one end of a uniform tube, and plane sound waves are set up in the tube 
by means of a source of sound at the other end. ‘The stationary-wave pattern 
which is formed in the tube is explored by means of a microphone, and 
measurements are made of the ratio of the maximum to the minimum pressure 
in the system and of the distance of the first minimum from the face of the sample. 
From these two quantities the two components of the acoustic impedance of the 
material are obtained. 

The derivation of the algebraic expressions which relate the resistive (R) 
and reactive (X) components of the acoustic impedance of the sample to the 
ratio of maximum (Pmax) to Minimum (P,jin) oscillatory pressure in the stationary 
wave and to the distance (dy) of the Nth minimum from the face of the sample 
is given in standard text-books on sound * for the ideal case of negligible dissipa- 
tion in the stationary wave. ‘Thus it can be shown that the acoustic impedance Z 
of the sample is given by 


Z=RYjX Spetanh(y 75), (1) 


where p is the density and c is the velocity of sound in the air. The quantities 
y and 6 are obtained from the measured values of | Pmax|/|Pmin| and of dy by the 
relations 


ee | / | Pmin | =coth {ne Le Seed (2) 
Teas oe ORES 


where A is the wave-length of sound in the tube. 

The assumption of negligible attenuation in the tube is not justifiable where 
precise measurements (i.e. to about 1% in R and X) are required. The principal 
effect of finite attenuation in the tube is to cause the maxima and minima to alter 
progressively in magnitude with distance from the sample. An example of this 
is shown in figure 1, which gives a graph of oscillatory pressure against distance 
from the sample for a specimen of rock-wool 1” deep at 1500 c/s., backed 
by a massive piston. ‘Ihe theory of this effect has been forinulated by Davis 
and Evans (1930) and by Hall (1939), who show that to a first approximation 
the maxima can be regarded as substantially unaffected, whilst the minima 
increase linearly as the distance from the sample increases. 

The finite attenuation in the tube leads also to a shift in position of the minima. 
The existence of this shift must have been clear to Davis and Evans, but its 
magnitude does not yet appear to have been discussed in detail. Neglect of the 


* Morse (1936). Note that Morse describes the time dependence by exp (—jw#), whereas 


here we use the alternative exp(jwt). The difference results in a corresponding change of sign 
in the phase angle. 
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shift can cause appreciable error in precise measurements of impedance, and 
an outline of the effect will, therefore, be given. 

Assume that the acoustic disturbance in the tube may be represented by means 
of two oppositely and axially directed trains of plane waves of frequency w/2z, 
wave-length constant 8 (=27/X) and attenuation constant «. The expression 
for the oscillatory pressure p, at a section of the tube distance x from the face 
of the sample is 


p,=Asinh {ax+y +7(Bx +65)}, 


whence the amplitude of the oscillatory pressure is given by 


P,=B [cosh 2(ax+y)—cos2(Bx+8)}?, 9 ....e, (4) 
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Figure 1. Typical graph of pressure against distance for sample. 


where B is a constant independent of «, and where y and 6 are quantities which 

are associated with the conditions of reflection at the sample, and which satisfy 

equation (1); thus e-*” is the ratio of the amplitudes of reflected to incident 

wave at the sample, and 26 is the phase change accompanying reflection. 
Stationary values of p, occur when 


, asinh 2(axn+y)+Bsin2(Bx+d)=0, = = — ...... (5) 
1.€. 
sin 2(Bx +8) = — sinh 2 (a+). 


If «/B is much less than z, the solutions of this equation can be written 


(og 


2(8x+8)=2Nr— 5 


sinh 2 (ax+y) 
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and 2(Bx +8) =(2M + 1) + 3sinh 2(ax-+y), 
whence, for p, a minimum, the distance dy to the sample is given by 
dA 
=1NA—=- —aspsinh2(ax+y), = a naan 6 
dy=3NA an — zppsinh2 (wx+y), (6a) 


and for p, a maximum, the distance Dy, to the sample is given by 


2M+1, 8A 
pee ed 
TT 


Dy a uintt 2 (ase yy eee (6) 


2p? 
In the measurement of impedance, particular interest centres on the position 
of the first minimum. Rearranging (6a), 


5 

2a BRR 
Comparison of (7) with (3) shows that the effect of finite attenuation in the tube 
is to reduce the distance between the sample and the first minimum by an amount 


dy/A=3N ——~sinh2(ax+y). «= neeee'e (7) 


oR sinh 2(ax«+y). 

The effect of attenuation on the ratio | Pyax|/|Pmin| may be found by substituting 
from (6a) and (66) in the expression for p,. Thus for the ratio of the Mth 
maximum to the Nth minimum : 


2 cosh? («Dy+y) — 
| Prax | / | Pmin| = 


2B sinh? 2(«Dy+ ilk 


oF 
The importance of the effect of attenuation in the tube depends on the magnitude 
of the attenuation coefhicient. For tubes with smooth walls and of diameters 
greater than six inches, the values of « are very small. Experiments on the 
apparatus which 1s to be described later show that for frequencies up to 3000 c/s. 
the attenuation constant for a tube of diameter 4-45 cm. is given by 1-37 x 10-54/f 
for room temperatures of 22° c. This figure corresponds to air of low humidity, 
and some modification may well be expected for damp air (Knudsen, 1933). 
The above figure is in good agreement with the theoretical attenuation constant 
deduced from the well-known Helmholtz-Kirchhoff analysis of the passage of 
sound waves along smooth-walled, circular pipes (viz. 1-33 x 10->y/f at 23° c.). 
Where y is small compared with «d,, the shift, which is a maximum at low 


frequency, is less than 0-1% of d,. As y approaches unity, however, the 
correction becomes 


2 sinh? (ady+y)+ a5 int 2(ady+y) 


NE B sinh 2y, 
and is appreciable at low frequencies. Figure 2 shows a plot of sinh 2y against 
the ratio | Pmax | /|Pmin| expressed in terms of L =201log,9|Pmax|/|Pmin |, and included 
on the graph is a table of «/26? for various frequencies in the range of the apparatus. 
For a stationary wave with |Pmax|/Pmin| =2, i.e. L=6-0, sinh 2y = 1-34 and the 
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correction to dy becomes 2-7 mm. at 100 c/s. For a tube of smaller diameter 
this correction would become more important. 


For the values of « and f appropriate to the present apparatus, the ratio of 
maximum to minimum pressure differs negligibly from 


cosh («Dy +)/sinh(ady +y). 


Successive maxima lie on a cosh curve and minima lie on a sinh curve. The 
maxima vary very little with distance, whilst the minima vary approximately 
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Figure 2. Correction to position of minimum, for effect of attenuation in the tube. 


linearly with dy, except where y is large, in which case the total variation 
of minima with distance is correspondingly small. Account may conveniently 
be taken of the effect of attenuation on the ratio of maximum to minimum 
by plotting | Pmax| and | Pmin| against distance from the sample and extrapolating 
the envelopes to the position of thesample. The ratio | Pmax|/| Pmin| extrapolated 
in this fashion then satisfies the equation (2), which is deduced by neglecting 
attenuation. 
$3) DESCRIPTION OF APPARATUS 


For most practical problems in which sound-absorbing materials are used, 
it is necessary to know the magnitudes of the components of the acoustic impedance 
to within 5%. A better accuracy than 5% is required where measurements 
of impedance are to be used in the discussion of the detailed mechanism by which 
sound is absorbed in the material. It is, therefore, desirable that the measuring 
apparatus should be capable of providing results with an accuracy of 1 °% in both 
resistive and reactive components for frequencies in the range of 100 to 5000 c/s. 
At lower frequencies than 100 c/s., available materials have such poor absorption 
that accurate measurements of impedance are at present of little practical interest. 

Where determination of the normal sound-absorption coefficient to within 
‘0-1 is concerned, measurements present no great difficulty, since the ratio 
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|Pmax|/|Pmin| alone is involved. For measurements of impedance a knowledge 
of the position of the minima is also necessary, and here the requirements are 
stringent. The impedance of a layer of almost any material used for absorption 
of sound possesses a large negative reactance term for frequencies low in the 
acoustic range. ‘The reactance is associated with the finite depth of the 
absorbing material, and is due to the high elasticity of the air trapped in the layer. 
Thus R/pc=1-01, X/pc= —16-7 for a 1” thick layer of rock-wool at 100 c/s. 
In this case d, is about 83 cm., A about 345 cm. and d, must be measured to within 
0-3 mm. in order to give R and X with an accuracy of 1%. If such precise 
measurement is to be significant, the following conditions must be satisfied : 


(i) the bore of the tube must be uniform and the walls smooth and substantially 
rigid ; 
(ii) the sound field in the tube must consist of plane, uniform, axially directed 
waves ; 
(iii) the face of the specimen of material must be substantially plane and mounted 
normal to the axis of the tube; 


From Oscillator 


Figure 3. General arrangement of measuring apparatus 


(iv) the microphone used for exploration of the sound field must notfitself 
appreciably affect the field and must be sensitive and stable; 

(v) the measurement of position of microphone orifice must be accurate to 
about 0-1 mm.; 


(vi) the pressure measured with the microphone must arise only from sound 
of single frequency. 


The apparatus was designed and constructed so as to satisfy these require- 
ments. The conditions under which plane, axially directed waves may be set up 
in a tube are well known. For a circular tube Rayleigh has shown theoretically 
that serious departures from uniform plane waves will not occur provided that 
the diameter of the tube is less than 0-58 A, where A is the wave-length of the sound 
of highest frequency used in the tube. Experiment shows that in practice 
the field is substantially uniform until the diameter of the tube reaches nearly 
twice this value. Thus, taking Rayleigh’s limit, if measurements are required 
at frequencies as high as 5000 c/s., the diameter of the tube should not exceed 
4:0 cm. A diagram of the apparatus is shown in figure 3. The tube (B) has 
an internal diameter of 4-45 cm. and a smooth brass wall 1” thick. The 
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maximum variation in cross-sectional area along its length amounts to less than 
01%. The specimen of absorbent material (C) is mounted in a detach- 
able specimen tube which is closed at the rear end by a massive (1000 gm.) 
close-fitting brass piston (D). The positions of the piston and specimen are 
adjusted so that the front face of the specimen lies accurately in the plane (H) 
of the flange of the specimen tube. The sound enters the tube through a conical 
end-piece (J) from a moving-coil loud-speaker unit (A) (Vitavox: type N), which 
in turn is supplied with electrical power from a valve oscillator of conventional 
type. The sound field in the tube is explored with a moving probe-microphone 
(E), the position of which, along the axis of the tube, is measured by means of 
a cursor sliding over an accurately engraved steel metre scale (Chesterman: 
No. 1367) (G) divided in millimetres. 

The mode of construction of the probe-microphone is important. The 
microphone consists essentially of a long tube of stainless steel of internal diameter 
0-241 cm., external diameter 0-316 cm. and length 104cm. One end of the tube 
is connected to a small cavity, of volume about 2 cm, in front of the conical 
diaphragm of an otherwise hermetically sealed rochelle-salt, piezo-electric 
microphone. ‘The passage of sound waves to the microphone by any path other 
than along the tube is prevented by the presence of the massive (}” thick) shell 
of the microphone casing, which is made from a steel bar of diameter 3’’ bored 
out so as to take the microphone unit. The other end (F) of the probe-tube 
is open and communicates with the sound field. ‘The response of the microphone 
at any one frequency is proportional to the oscillatory pressure at the mouth 
of the probe-tube, and since the impedance presented by the mouth of the probe- 
tube is high, the acoustic pressure at the mouth of the probe is effectively that 
existing at the same position inthe sound field when the microphone is withdrawn. 
The modification of the sound field caused by obstruction due to the probe-tube 
is negligible, since the area of cross-section of the probe-tube is less than $% 
of that of the main tube. The construction of the apparatus is shown in more 
detail in figures 4 and 5. Figure 4 is a side view of the trolley (with indicating 
meter removed) and of the driving end of the tube. Figure 5 is an end view of the 
trolley, looking from the direction of the tube. 

Vibration which exists in the wall of the probe is prevented from reaching 
the microphone by the interposition of a sleeve (K) of cycle valve-rubber between 
the two components. Tests were made of the response of the microphone 
in the normal condition and with the orifice of the probe-tube closed with a short 
length of brass rod. A typical test at 1000 c/s. showed that the response of the 
microphone was reduced at least by 50 db. when the orifice was closed. ‘The 
microphone has proved in practice to be stable and sensitive, and no difficulty 
is experience in determining the acoustic pressure at maxima or minima in the 
tube to within 0-1 db. for any frequency in the range from 100 c/s. to 5000 c/s. 

The probe-tube passes into the main tube through a long hole bored in the 
conical end-piece, and for the most part of its length rests on the bottom of 
the main tube. The probe-tube, rather than the microphone casing, is clamped 
rigidly to the saddle of the microphone trolley (L), whilst the microphone 
casing is resiliently supported above the trolley by means of the sponge-rubber 
(M) and the ring (N). The scale cursor (O) is mounted close to the fixing point 
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of the probe, and in consequence the scale readings give accurately the successive 
positions of the probe orifice. 

The microphone trolley is mounted on four wheels: two grooved to fit one 
rail of the angle-brass track (Q in figure 5) and the other two ungrooved. Coarse 


Figure 4. Detail of driving unit and microphone trolley. 


adjustment of position of the microphone is made by pushing the trolley bodily 
along the track, and fine adjustment by means of the knurled knobs (R), which 
form extensions to the near-side wheels. A four- rather than three-wheeled 
construction for the trolley was used to provide greater stability of the trolley. 


Figure 5. End view of microphone trolley. 


The frame of the trolley was made to be sufficiently flexible for each of the wheels 
to remain in good contact with the rails. 

The output voltage from the microphone is amplified by means of a multi- 
stage valve-amplifier, and is then passed through an electrical analyser, tuned 
to the working frequency, to an indicating meter (P in figure 5) of the moving-coil, 
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rectifier-type (2}’”) which, for convenience, is mounted on the trolley. The 
amplifier effectively excludes such harmonics of the fundamental frequency 
as may exist in the sound field. 

In the plate facing p. 264, (a) shows a general view of the apparatus, and (b) gives 
a more detailed view of the driving end of the tube and of the microphone trolley. 

Typical operation of the apparatus is as follows :—A sample of acoustical 
material is cut by means of a circular pastry-cutter so as to fit accurately into 
the specimen tube. The front face of the material is adjusted so as to lie 
accurately in the plane of the face of the flange of the specimen tube, and the 
brass piston (D) is moved so as just to touch the rear face of the specimen. The 
sound field of appropriate frequency is set up in the tube; there is, of course, 
no necessity to “‘tune”’ the tube since the form of the sound-field, looked at from 
the “‘sample”” end of the tube, depends only on the nature of the reflection 
at the sample. The reading of the indicating meter, which is scaled in db., 
now corresponds to the pressure at the end of the microphone tube. The 
scale (G) and cursor have previously been set up so as to give a reading of zero 
when the axis of the upturned end portion of the probe-tube is exactly in the 
plane of the face of the flange to which the specimen tube attaches. Correspon- 
ding readings of indicating meter and of scale give the acoustic pressure in terms 
of the distance from the face of the sample. The trolley is moved smoothly 
away from the specimen, and several successive readings are taken of pressure 
at maxima and minima. From these measurements the ratio at the sample 
of maximum and minimum pressure is obtained, if necessary by extrapolation 
of a graph of the ratio, in linear units. 

Positions of minima relative to the sample are now found, and are conveniently 
determined from the mean of two readings, equally spaced about the minimum, 
each corresponding to the meter pointer passing through an arbitrarily chosen 
value just in excess of the minimum. The curve of meter reading against 
distance shows a high degree of symmetry in the region of the minimum. When 
the number of minima in the tube allows, the positions of eight minima (four 
at the specimen end and four at the distant end of the tube) are measured. ‘These 
data allow not only a convenient check on each other, but also make it possible 
to determine directly the wave-length of the sound, which enters into the 
calculation. 

Thus for a sample of rock-wool (stillite, 5 lb./cu. ft.) 1” thick at 1500 ¢/s., 
successive maxima and minima are shown in the following table. ‘The maximum 
remains closely constant at 13-6 db. (the reference level is arbitrary, but the same 
throughout the measurements). The acoustic pressure is plotted directly against 
distance from the sample in figure 1, and extrapolation back to the sample shows 
that, with no dissipation in the tube, the minimum pressure would have a value 
of 0-99. The ratio | Pnax|/|Pmin| of maximum to minimum pressure is, therefore, 
4-83, i.e. 13-7 db. oe 

The positions of eight minima are also shownin the table. By pairing of 
readings, four values for 2A are obtained and lead to a value of A of 23-00 cm, 
and to a corresponding velocity of sound of 345-0 m./sec. ‘The position of the 
first minimum (d,) can be obtained directly from the first reading, viz. 2-895 cm., 
from the sample. This value can be checked by means of the other figures 
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Position of minima of pressure and magnitude of maxima and minima 
for a typical sample of rock-wool 


Stillite, 5-0 lb./ft®, 1” depth 


ae Intensity Distance oe a 
N | Pmax |/ PP natn | from sample colunin (4y—w- 3) 
(db.) dy(cm.) (cm.) (cm.) 

1 0-0 2°895 2°895 
13-6 

2 0:2 14-39 2°89 
13:6 

3 0:5 25°89 2°89 
3-6 

4 0:75 37°39 2°89 
13-65 Sa 

5 1-0 48-89 45-995 2°89 
13-65 

6 hows 60:395 46-005 2°895 
13-65 

7 1:45 71-895 46-005 2°895 
iLS7/ 

8 1:65 83-395 46-005 2°895 

Mean 46-00 


since the successive readings dy should differ from the first reading only by some 
known multiple (V — 1) of the measured half-wave-length. Thus the last column 
of figures in the table is found by subtracting from the measured values of dy 
the appropriate multiple of the half-wave-length, and it is apparent that the 
agreement in the resulting values of d, is excellent. 

Using the above data, 


[Dae } | Print| = 4-83 or 13-7 db., 
d, =2-895 cm, i.e. d,/A=0-1259. 


The correction to d, as calculated from (7) is 0-001 cm. The correction 
is not significant at this frequency. The corresponding values of R/pc and 
of X/pc are 

Ripc= 0399: 
X/pc= —0-927. 


The length (86°5 cm.) of the tube used in the present apparatus is such that 
two or more minima can be located only for frequencies of 300 c/s. and above. 
Direct measurement of A from the interval between minima is therefore impossible 
below 300 c/s., and A must be found by other methods. Careful measurement 
of the wave-length, and hence of the velocity of sound, in the tube in the range 
of frequency 400 to 4000 c/s. has shown that the velocity C’ follows closely 
the Helmholtz-Kirchhoff law, 
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where C is the free-space velocity, r is the radius of the tube, f the frequency 
and y’ a function of the kinematic viscosity and diffusivity of the air (Rayleigh) 
equal numerically, for r in centimetres, to 0-577 at 23°c. 

The usual procedure adopted with the apparatus is to measure C’ accurately 
at one frequency (3000 c/s., at which many minima are available) and to apply 
a small correction calculated from the Helmholtz-Kirchhoff relation, at other 
frequencies. By this means the wave-length in the tube can be determined for 
frequencies below 300 c/s. and measurements of acoustic impedance can be made 
at 100 c/s., at which frequency the tube-length is only one-quarter of the 
wave-length. 


§4. REDUCTION OF RESULTS 


Apart from the corrections for attenuation, the values of |Pmax|/|Pmin| and 
of d,/A are provided directly by the apparatus, and there remains the process 
of calculating the real and imaginary components of the impedance from them. 
This problem reduces to the solution of equations (1), (2) and (3). These 
equations can be solved algebraically, and thus it can be shown that 


on 
1+7?+(1—7?)cos4md,/A’ 


—(1—r?) sin 4zrd,/A 


Ripe = 


ites = cosara a: 
where 7—tanhy = Loni 
max 


Numerical evaluation of R/pc and X/pc is , however, tedious, especially when 
results are to be worked out at several frequencies. It is convenient, therefore, 
to have graphical methods of solution available. In this respect the position 
is fortunate. H. J. Sabine (1942) has constructed charts with rectangular axes 
R/pc and X/pc, on which are plotted families of curves of constant | Pmax!/|Pmin| 
(in decibel units) and of constant d,/A. (He plots for D,/D,, where D, has the same 
significance as our d, and D,=4A.) ‘These charts are very useful for measurements 
at medium frequencies on materials of depth greater than 1 inch. For lower 
frequencies and small depths, e.g. where R/pc may be of the order unity and 
X/pc of the order — 20, the range of curves is inadequate. 

A second chart suitable for the calculation of R/pc and X/pc has been provided 
recently for the solution of a closely analogous problem in connection with radio 
wave-guides. Such a chart has been published by P. H. Smith (1939) and 
by Willis Jackson and L. G. H. Huxley (1944). In this chart, functions of 
| Pmax| /|Pmin | and of d,/A are provided as polar co-ordinates, and curves of constant 
R/pcand X/pc are plotted thereon. The resulting families of curves corresponding 
to all values of R/pc from zero to infinity and of all values of X/pc from minus 
to plus infinity lie inside a circular boundary of finite radius. ‘This type of chart 
is easily constructed from the data provided in the above references and is very 
convenient to use in practice. 

Figure 6 shows the results of typical measurements of acoustical impedance 
obtained, in the manner described above, over the range 100 to 5000 c/s., 
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of a 1” thick specimen of rock-wool (stiltite: 5 Ib./ft?) backed by a sub- 
stantially rigid piston. 


Acoustic Impedance (in terms of pe.) 


Frequency ¢C/s. 


100 200 400 600 1600 3200 6400 


Figure 6. Typical measurements made with the apparatus of acoustic impedance of a 
soundeabsorbing material. 
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AN EXPERIMENTAL INVESTIGATION OF THE 
REFLECTION AND ABSORPTION OF RADIATION 
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National Physical Laboratory, Teddington 
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ABSTRACT. ‘The paper, which was circulated as a confidential report to the Radio 
Research Board in October 1944, describes measurements on the reflecting power of various 
surfaces at a wave-length of 9 cm. and angles of incidence ranging from 80° to 45°; these 
included level and uneven bare ground, vegetation-covered ground, tap water and a 4% 
salt solution approximating to sea water. Specular reflection was found to occur only 
from very level surfaces ; the absorptions of these surface media were measured, and 
from the combined measurements of reflection and absorption their electrical constants 
were derived. 

Rough surfaces, either of bare ground or vegetation-covered ground, gave values of 
reflection coefficient in general agreement with the optical rule that regular reflection 
is only observed from an uneven surface if the product of the depth of the surface irregu- 
larities and the cosine of the angle of incidence is a small fraction of the wave-length. 
If this fraction exceeded 4A, the values of reflection coefficient measured were about 0-1. 

The calculated values of reflection coefficient corresponding to dry soil, wet soil and 
sea water for angles of incidence varying from 0° to 85° are given in an appendix to the paper. 


SNE RO DUG asLON 


HE reflection of radio waves has been investigated experimentally at a 
wave-length of 9cm., using a wide range of reflecting surfaces; these 
included bare ground, both level and in ridges, ground covered with vege- 
tation, fresh water, and salt water approximating in composition to sea water. 
When specular reflection occurs, the electrical constants of the material forming the 
reflecting surface can generally be estimated from measurements of the reflection 
coefficient. At the high frequencies involved in this investigation, however, 
the reflection coefficient is determined almost entirely by the dielectric constant ; 
consequently reflection experiments, while providing a reasonably accurate 
indication of the dielectric constant, only enable a very rough estimate to be made 
of the conductivity of the surface material. ‘The absorption in the material, on the 
other hand, is then dependent mainly on the conductivity. Experiments were 
therefore made on the absorption, at the same wave-length, of the materials of the 
various reflecting surfaces, from which the conductivities could be determined. 
Since both reflection and absorption experiments are needed to obtain the 
electrical constants, the sequence adopted in the presentation of this paper is to 
describe each of the experimental methods before giving any of the results obtained 
‘from either, rather than to divide the paper into two distinct sections. The 
results of both reflection and absorption measurements are then given together, 
for each type of surface in turn. 
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§2. EXPERIMENTAL METHODS 


A. Measurement of reflection coefficients 


The measurement of the coefficient of reflection for radio-waves is simplified at 
centimetre wave-lengths by the possibility of using aerial systems which give beams 
of radiation which are so narrow that, except at angles within a few degrees of 
normal or grazing incidence, the signals due to the direct and the reflected rays can 
be measured separately by suitably adjusting the aerial. The magnitude of the 
reflection coefficient is thus obtained from the ratio of two measurements of field 
strength; no indication of its phase angle is obtained. 

A diagrammatic representation of the arrangement of the apparatus is shown in 
figure 1, while the plate shows a photograph of an experiment. 
were constructed of tubular scaffolding, and were 15 metres (50’) apart, with 
platform levels of 1-7 m. (5:5’), 4:3 m. (14’) and 6-7 m. (22’) above ground. The 
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Figure 1. Diagrammatic representation of ground reflection experiments. 


1:2-metre diameter aperture parabolic mirrors, each of which was used with a 
half-wave dipole aerial at the focus, were so mounted that the aerials were 1-4 metre 
above the platform level, thus giving angles between the reflected ray and ground of 
22°, 36°-5 and 46°-5 (angles of incidence of 68°, 53°-5 and 43°-5). Between the 
two towers a symmetrically disposed area of ground with sides 7 x 10 metres was 
prepared to obtain the several types of reflecting surface used. For experiments 
on water, a National Fire Service canvas dam, giving a water surface of about 
5 metres diameter, was placed symmetrically between the towers. 

The sender consisted of a CV. 67 valve oscillator mounted directly on the back 
of one of the paraboloids, and a crystal monitor was used to ensure constancy of 
output during each set of observations. Two receivers were used: one (Saxton 
and Grace, 1946) was a bolometer mounted directly behind the receiving para- 
boloid, the other was a field-strength measuring set (Grace, 1946). When using 
the bolometer, reflection coefficients greater than about 0-2 could be obtained with 
an instrumental accuracy of about + 0-01, from the square root of the ratio of the 
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values of the power received with a direct and with a reflected signal. When 
using the field-strength measuring set, the ratio of the direct and reflected 
received signals, in decibels, could be measured with an accuracy of about + 0-5 db. 
if less than 6 db., and an accuracy of about +1 db. if greater than6db. Thus the 
bolometer was considerably more sensitive than the field-strength measuring set, 
except under conditions where the reflected signal was small. 

The procedure followed at each level was to adjust the sending and receiving 
mirrors in turn until a maximum signal could be measured from the direct ray; 
next to readjust the mirrors for a maximum signal from the reflected ray, and 
finally to repeat the measurement of the direct signal to verify that the sender 
output had remained constant. 

Apart from instrumental inaccuracies in measurement, there are a number of 
possible sources of error in the results obtained. Experiment showed that the 
polar diagrams of the mirrors were such that any signal due to the direct ray, when 
measuring the reflected ray, would be more than 40 db. down, which virtually 
eliminated one potential source of error. Experiment also showed that a distance 
of 15 metres between the mirrors was somewhat less than the minimum distance at 
which the transmission law holds, that the received signal varies inversely as the 
distance from the sender. A careful series of measurements was therefore made 
to determine the variation of received signal with distance from the sender over the 
range concerned, and the values thus obtained were used in place of the geometrical 
values when correcting the measurements of reflected signal for the path difference 
between the direct and reflected rays. A slight, but perceptible, interaction 
could, moreover, be detected between the receiver and sender mirrors, due to 
reflection from receiver to sender mirror when the two were directed at each other ; 
calculation showed that this was unlikely to cause an error of more than 0-01 in the 
reflection coefiicient. ‘The agreement obtained between measurements at three 
different angles of incidence, and consequently three path differences, under 
conditions of specular reflection, makes it improbable that any of these sources of 
error was serious; consequently, the reflection coefficients measured with bolo- 
meter are believed to be accurate to + 0-02. Had the distance between the towers 
been made great enough to reduce all the sources of error mentioned to negligible 
proportions, both the height of the towers and the area of ground requiring pre- 
paration would have become undesirably large. 

To check the reliability of the method, a level platform of planks was set on the 
ground midway between the towers, and on it were placed six 1 x 2-metre iron 
sheets, making an area of 3x4 metres. The reflection coefficient of this was 
measured at angles of incidence of 53°-5 and 43°-5. ‘The reflection coefficients 
obtained varied between 0-95 and 1:03, with a mean value of 0:98. The obser- 
vations were made with the field-strength set, and the possible instrumental error 
was consequently +0-5db. The values obtained, therefore, are in sufficiently 
good agreement with the expected value of unity to indicate that none of the 
possible sources of error in the experimental method was serious. 


B. Measurements of absorption 
The absorption of radiation by a screen of the reflecting medium used in the 
previous experiments was measured directly by interposing a filled container 
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between sender and receiver. The arrangement used is shown diagrammatically 
in figure 2. The sender was exactly as used in the reflection experiments, and was 
set up about 6 metres away from a table on which stood the test sample and the 
receiver. The latter consisted of a horn in which was a bolometer, forming the 
field-strength measuring device described by Saxton and Grace. A container 
consisting of a wooden trough was made having ends of 6mm. plate glass, 46 cm. 
square, one of which was fixed and the other movable, thus forming a container 
in which a sample 46 cm. square and of any desired thickness up to 30 cm. 
could be tested. This container could be interposed between the sender and the 


receiving horn. 
The experimental procedure was to measure the field strength at the horn 


when the transmission path from the sender was unobstructed, and then to repeat 
the measurement with the container filled with a chosen thickness of the substance 
under test, placed symmetrically to obscure the horn from the sender. Successive 
measurements were made with the container at various distances from 1 to 6 cm. 


QISTANCE ADJUSTABLE, FIXEO GLASS PLATE. 


/¢tm To 6cm. 


AOTUSTABLE GLASS PLATE. 


RECEIVER 


SENDER 


Figure 2. Diagrammatic representation of absorption experiments 


from the aperture of the horn. In most cases considerable variations were found 
to be present in the received signal as the distance was varied, due, presumably, 
to reflection from the end of the container nearest the horn, and the mean value 
over a range of distance equal to half a wave-length was taken as the true trans- 
mitted signal. ‘The thickness of the obstruction was then altered by adjusting the 
movable glass end-plate and refilling the container, and the experiment repeated. 
Thus the attenuation in a known thickness of the material under test could be 
obtained without reference to the properties of the container. 

The method of measurement was subject to limitations; if the thickness of 
material under test were too small, multiple reflections between the two end plates 
would bea serious source of error, while if the thickness were too great the accuracy 
of measurement would be inadequate. With thicknesses of material such that the 
attenuation was between about 5 and 15db., however, the measurements are 
believed to be accurate to + 0-5 db. 


§3. EXPERIMENTAL RESULTS: SPECULAR REPEECTION 


A. Reflection and absorption from bare ground 


An area 7 x 10 metres midway between the towers was cleared of turf, forked 
over to a depth of about 10cm., carefully raked level, and finally smoothed by 
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dragging a plank across it. Preliminary experiments showed that the surface thus 
obtained gave specular reflection. This conclusion was based on the fact that 
a single value of dielectric constant, K, and conductivity, o, could be chosen which 
gave calculated values of reflection coefficients which, at all angles and with both 
horizontal and vertical polarizations, agreed within the accuracy of experiment 
with the measured values. Curves connecting the reflection coefficient with K 
and o, of the type described by McPetrie (1938), were drawn for each of the angles of 
incidence available to simplify the choice of the appropriate values of K and to fit 
any given set of test results. It quickly became apparent that the dielectric 
constant rather than the conductivity determined the reflection coefficient, and 
that changes in the wetness of the ground caused large variations in the value 
obtained. ‘The decision was therefore made to concentrate on the two extreme 
conditions of humidity. 

The penetration of the wave into the ground on reflection with centimetre 
wave-length radiation is only a few centimetres, which is much less than occurs at 
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Figure 3. Attenuation cf A=9 cm. radiation in various media. 


the longer wave-lengths; consequently extreme conditions, particularly of dry 
ground, will be more often encountered with centimetre waves than with the longer 
wave-lengths. 

An immediate consequence of the relative unimportance of the conductivity in 
determining the reflection coefficient is that the change in phase on reflection 1s 
about 180° for horizontally polarized waves at all angles of incidence; it is about 
180° for vertically polarized waves at angles of incidence appreciably greater than 
the Brewster angle, and at angles of incidence appreciably less than that angle there 
is no appreciable phase change on reflection. 

(a) Dry ground. A prolonged period during May 1944 in which no rain fell 
gave an opportunity to measure the ground reflection under conditions where the 
surface layer was as dry as is ever probable; no trace of moisture was apparent, and 
the top surface was almost as loose in consistency as sand. 

The decision to measure the ground absorption was not taken until some time 
after the reflection tests had been made. Consequently, to obtain ground in as 
nearly as possible the same condition, a sample of surface soil was placed under 
cover and left for about the same period as had elapsed without rain before the 
reflection test. The results of the absorption tests on dry soil are shown in 
figure 3 (a), and indicate an attenuation of 36 db. per metre. 
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The most suitable values of K and o to fit both reflection and absorption 


measurements were 
ears g=3x 10° e.s.u. 


The reflection coefficients calculated from these values, together with the 
measured values, are given in table 1. 


Table 1. Reflection coefficients for very dry ground, A=9cm. 


Angle of Vertical polarization Horizontal polarization 
incidencé, {$+ _—_ | __ 
(degrees) Calculated Observed Calculated Observed 
68 0-18 0:20 0-48 0-47 
53505 0-015 0-03 0-33 0:33 
43-5 0-08 0-09 0:26 0:27 


(b) Wet ground. ‘The preliminary experiments showed that the value of K 
corresponding to the reflection coefficients usually encountered varied between 
5 and 15, according to the degree of wetness of theground. ‘The opportunity was 
taken during a rainy spell in February to obtain a series of measurements immedi- 
ately after a twelve-hour period of steady rain, when the ground was well saturated. 
The field-strength measuring set was used, so that the results obtained were of less 
accuracy than the others given in this report for conditions of specular reflection, 
all of which were taken using the bolometer. Accordingly, a further series of 
measurements were made in which, since the weather conditions were not suitable, ° 
the ground was kept wet artificially. For this purpose a spray branch of the type 
used to deal with electrical fires was used with a fire hose to water the site, and the 
equivalent of about four inches of rain was delivered at intervals over a period of 
48 hours. A slight further watering was given about half an hour before each test. 
The resulting condition of the ground was such that it was saturated with water, 
but was not waterlogged, 1.e. no puddles were present on the surface. 

To obtain a similar sample of soil for absorption measurements the container 
was filled with soil and water, and allowed to stand until the surplus water had 
drained off. It was found that this method gave closely repeatable results, and the 
assumption that the sample in the container was then under the same conditions 
of humidity as the ground during the reflection experiments seemed justified. 
The results of the absorption tests on saturated soil are shown in figure 3 (a), and 
indicate an attenuation of 220 db. per metre. 

The most suitable values of K and o to fit both reflection and absorption 
measurements were 


K=24, o=0x 109 esin. 


The reflection coefficients calculated from these values and the measured 
values are given in table 2; calculated values assuming K=25, o=0, are also 
included to demonstrate the comparatively small contribution made by the con- 
ductivity to the reflection coefficient. 

The ground surface on which the above experiments were made is composed of 
a fine sandy loam. Other types of ground might give appreciably different results. 
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The conductivities of both dry and wet soil were measured at a frequency of 
1 Mc/s. (A=300m.) by the Electricity Division, N.P.L., using a laboratory 
method. ‘The values obtained were 3 x 10° and 8 x 107 e.s.u. for dry and wet soil 
respectively. ‘These agree with the conductivities given by Smith-Rose (1933) 
for the soil at Teddington in extreme conditions of low and high humidity, and 


Table 2. Reflection coefficients for saturated ground, A=9cm. 


Vertical polarization Horizontal polarization 
Angle of 
Pee Calculated Observed Calculated Observed 
d 
ee) tee lees Heavy | Watered] K=24 | K=25 | Heavy | Watered 
c—6><107)|Fao—0 rain | by hose | c=6x10°9| o=0 rain | by hose 

68 0:31 0:31 0:28 0:32 0:85 0:86 0-90 0.86 
53°5 0:50 0:50 0:50 0:50 0:78 0-79 — 0:78 
43°5 0:57 0:57 0:58 0:58 0-74 0:75 0:72 0-74 


thus show that the site used was not abnormal. Smith-Rose and McPetrie (1934) 
have measured the constants of the soil at Teddington by a laboratory method at a 
frequency of 200Mc/s. (A=1:5m.), and they obtained values of K=22 and 
o=3 x 10%e.s.u., and of K =4 and negligible conductivity, under extreme wet and 
dry conditions respectively. 


B. Reflection and absorption from water 


Tests on the reflection of radiation of 9cm. wave-length were made on fresh 
and salt water by placing a circular canvas dam, lent for the purpose by the National 
Fire Service, midway between the towers. ‘This was filled to a depth of about 
30cm. with tap water or with a 4% solution of coarse salt, approximating to sea 
water. ‘The diameter of the base of the dam was 5-5 metres (18 feet), and that of 
the water surface about 5 metres. The area available for reflection was conse- 
quently rather small, although experience had shown that it would probably be 
sufficient to give reasonably reliable results. "The accuracy was, however, less 
than that of the ground-reflection experiments described above. 

The absorption experiments were made using the same container as before, 
the cracks between the end plates and sides being caulked with Plasticine to render 
them watertight. ‘The thickness of liquid under test was obtained by measuring 
the volume necessary to fill the container in each case. 

(a) Tap water. The results obtained from the absorption tests on tap water 
are shown plotted in figure 3 (b), and indicate an attenuation of 380 db. per metre. 

The most suitable values of K and o to fit both reflection and absorption 


experiments were 


K=80, o=2.<10!°seis.u: 


The reflection coefficients calculated from these values, together with the 


measured values are given in table 3. . | 
(b) 4% salt solution. ‘The results from the absorption tests on a solution of 


4°/, coarse salt in tap water are shown plotted in figure 3(5), and indicate an 
attenuation of 1100 db. per metre. 
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Table 3. Reflection coefficients for tap water, A=9cm. 


Angle of Vertical ¢ olarization Horizontal polarization 
incidence 


(degrees) Calculated Observed Calculated Observed 


70 0:51 0-51 0:92 0-90 
55 0-69 0-67 0-88 0-88 
AAS o0:73 0-70 0°85 0-83 


The most suitable values of K and o to fit both reflection and absorption 


experiments were 
k=80, o=5-5x 10! e.s.u. 


The reflection coefficients calculated from these values, together with the 
measured values, are given in table 4. . 


Table 4. Reflection coefficient for 4°% salt solution, A=9cm. 


Angle of Vertical polarization Horizontal polarization 
incidence © |-ARAARAAAAA———_sF| 
(degrees) Calculated Observed Calculated Observed 


70 0-54 0:56 0-93 0-92 
55 0-71 Oo7t 0:89 0-89 
Am ss 0:74 0-73 0-86 0-84 


When making the absorption measurements on all types of material, 
observations were made using both radiation polarized parallel to and perpen- 
dicular to the edge of the container. It was found that a consistent difference of 
about 0:5 db. was present between the results with the two polarizations. A test 
on the empty container with the end plates removed showed that the container 
slightly increased the gain of the horn, and that this increase differed between the 
polarizations by the amount necessary to account for the discrepancy with the 
filled container. ‘The experimental points plotted on the curves of figure 3 are all 
mean values of the results obtained with each polarization. 


C. Reflection from turf-covered ground 


A few measurements were made on the reflection coefficient of ground covered 
with short grass. ‘The results obtained showed that, in general, specular reflection 
was not present. An exception was found when observations were made on a 
recently cut and rolled cricket wicket. The distance between the mirrors during 
the test was 21:5 metres and the height of the aerials 1-85 metres each, giving an 
angle of incidence for the reflected radiation of 80°. The results obtained are 
given in table 5, which also gives calculated values of reflection coefficient based on 
suitable assumed values of Kando. Values obtained on an adjacent cricket wicket 
on which the grass had been allowed to grow to a height of about 2 cm. are also 
included in the table. 

The value of K corresponding to the wet condition is lower than might be 
expected from the bare-ground experiments, since the site appeared to be well 
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soaked by the previous night’s rain. No opportunity occurred to take measure- 
ments under conditions such that the ground was certainly saturated. Accordingly, 


Table 5. Reflection coefficients of level grass-covered ground (Cricket 
wicket), A=9cm. Angle of incidence, 80°. 


Vertical Horizontal Assumed 
Chiditions polarization polarization values 
Calculated | Observed | Calculated| Observed | K ees 


Cut very short and 
rolled, dry. (No rain | > 0-47 0-46 0:79 0-79 8 yc ioe 
for at least 7 days) ; 

Cut very short and 
rolled, wet. (Several + 

oes 0-36 
hours’ rain in pre- 
vious night) 

Grass about 2 cm. 
high, wet. (Same day | > 0-36 0-51 0-86 0-83 © [PSs ie 
as previous test) uy 


0:37 0-86 0-86 © Sx Oe 


although the results appear to indicate that the value of K for ground covered with 
turf is less than that of bare ground of the same humidity, this conclusion cannot 
be given with certainty. 


D. Summary of specular reflection results 


‘The values of K and o appropriate to the various reflecting surfaces used in the 
investigation are brought together, for convenience, in table 6, which also gives 
the attenuation in the medium. 


(Pable-6 
Nature of medium K o (e.s8.u.) Hes 
Very dry sandy loam 2 SES OE 36 
Saturated sandy loam 24 6 x 10° 220 
Tap water 80 PSN” 380 
4°%, solution of coarse salt 80 Deo 101e 1100 
Dry turf 3 (e102) (50) 
(estimated) (calculated) 
Wet turf 6 (10°) (80) 
(estimated) (calculated) 


The values of attenuation and conductivity found for salt and tap water are 
somewhat lower than those previously obtained, using a different method, by 
Turner fora wave-length of 10cm., but agree with his results in showing that the 
ratio between the conductivities of salt and tap waters is much less than at low 
frequencies. The value of 80 for the dielectric constant of water is not precise ; 
a figure of 75 would fit the experimental results equally well. 

* 
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§4. EXPERIMENTAL RESULTS: DIFFUSE REFLECTION 


At centimetre wave-lengths, specular reflection is likely to be the exception 
rather than the rule, since the reflecting surfaces normally encountered are neither 
level nor free of vegetation. Experiments have been made to determine the 
reflection coefficient of level ground covered with vegetation, and of uneven 
ground. 

A. Effect of vegetation 

The results given in the last line of table 5 show the effect of grass about 
2 cm. high on the reflection coefficient of level ground at an angle of incidence 
of 80°. The effect of the grass was to decrease the reflection coefficient for hori- 
zontally polarized waves, and to increase it for vertically polarized waves. ‘The 
latter result was most unexpected, but the phenomenon was observed on several 
occasions (see also table 7). On longer grass the reflection coefficient for both 
polarizations was reduced; unmown grass in May, approximately 30 cm. high, 
gave coefficients of 0-19 with horizontally and 0-13 with vertically polarized 
radiation, at an angle of incidence of 80°. 


Table 7. Reflection coefficients of vegetation-covered wet ground, A=9cm. 


Appearance of site, and height Angle of incidence 
of vegetation of radiation, 0 
; 6=68" =r 6=43°.5 
Appearance Height | _ 
em ye ie Ne ree ee 

Bare = -0 0:32 | 0°86. | 0-50 | 0:78 | 0°58 | 0°74 
True leaves beginning to form, 

ground visible 3-4 0-40 | 0°50 | 0-44 | 0-55 | 0-47. | 0°56 
Dense clumps, ground showing 

in places 9-12) 0:13 8005) | 0-23 580258 0-33 O49 
Ground almost obscured 20=25 | 0°06 || 0°327 190-10) 0°39 5 ROc1 7a eO-4l 


Ground completely obscured 39=45. |) 0-04. | 0219501 50:05) 30265 NO -11e OLS 


The main experimental results were taken on the site between the towers, on 
which the bulk of the previous results on ground reflection had been made. The 
levelled area was sown fairly thickly with mustard seed (very approximately 
1000 germinations per square metre) to obtain a quick-growing carpet of vegetation. 
‘To ensure constancy of ground conditions, the site was watered thoroughly before 
each test as described in § 3, A(b). The results are given in table 7 for both 
vertical, V, and horizontal, H, polarizations. 

The site at the conclusion of the test was covered completely with a dense 
carpet of vegetation, the top surface of which was comparatively level. It is 
probable that the reflection coefficient would be even lower with a less uniform 
ground covering. ‘These results are in general agreement with those given by 
other investigators. 

B. Effect of uneven ground — 


The effect of an uneven ground surface was investigated by raking the bare 
surface of the site into a series of uniform ridges. When planning the experiments 
2, 
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the possibility of obtaining information applicable to the surface of the sea was 
borne in mind, and the system of ridges was made to simulate waves. Tests were 
made with these running along, across, and at 45° to the direction of transmission, 
using distances of 0-6 and 1-2 metres between successive crests, and various depths 
of between 5 and 16cm. between crest and trough. Table 8 shows the reflection 
coefficients observed when the ridges were at 45° to the direction of transmission, 
while table 9 gives those obtained when the ridges were either along or across the 
direction of transmission. In these tables D gives the distance (cm.) between 
the crests of successive ridges, and h the mean depth (cm.) between crest and 
trough, while the column marked “‘ Level” gives the estimated reflection coefficient 
for level ground of the same humidity as that under test. 


Table 8. Reflection coefficients of ground ridged at 45° to the direction of 
transmission, A=9cm. 


Angle of Vertical polarization Horizontal polarization 


incidence 


(degrees) | eye, | DP=O0 | D=120) D=120 D=60' | D=120'| D=120 


= laa he On ase len 140 p10. | pes 


68 0-08 0:07 0-09 0-13 0:65 0-14 0-18 0-30 


S305) 0-13 0:04 0-05 0:07 0:51 0:04 0:06 0-16 


Ae) 0-22 | 0:04 0-04 0:04 0-45 0:07 0:06 0-10 


Table9. Reflection coefficients of ground ridged along or across the direction 
of transmission, A=9 cm. 


Vertical polarization Horizontal polarization 

Angle of | | = 

incidence Along Across Along Across 

(degrees) | Level | D=60 | D=60 (| D=120| Level | D=60 | D=60 | D=120 

R14 5|) h—=16 — ANAS Sh 16 [p= PD 

78 0-23 0-2 0-1 0:3 0-86 0-4 0-2 0-4 
68 0-06 0-03 0:05 0-12 0-76 0:07 0:10 0:18 
53°5 0-28 0-03 0-02 0-06 0:64 0-10 0-12 0:16 
43-5 0:36 0-04 0:08 0-08 0:58 0-04 0:08 0-14 


The results show clearly how small an irregularity in the ground surface is 
sufficient to prevent specular reflection; the reflection coefficient, moreover, 
becomes erratic with respect to angle of incidence when it has fallen below a value 
of about 0-1. The values given for level ground, although consistent among 
themselves, only relate very approximately to the state of the ground in the ridged 
condition, since the length of time necessary to make all the measurements, 
several days, was sufficient to allow marked changes in ground humidity. The 
difference between the values for level ground given in table 8 and those in table 9 
are due to the observations recorded in the former table having been taken later 
in the year than those recorded in table 9, when the ground was appreciably 
drier. Values are given in table 9 for an angle of incidence of 78°. ‘These were 
obtained with the mirrors at ground level about 14 metres apart. ‘The accuracy 
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in these circumstances was not considered sufficient for measurements under 
conditions of specular reflection, but is adequate to give an indication of the reflec- 
tion present from uneven ground. 

The presence of a pattern in the signal reflected from ridged ground had been 
expected at the start of this investigation. Measurements were therefore made 
in which the mirrors, although always maintained at a distance of 15 metres apart, 


were moved so that the point of reflection to be expected from geometry was 


successively at different parts of the ridge system. Thus, when working on ridges 
running across the direction of transmission, the mirrors were moved together in 
the line of transmission to a number of successive positions so that the geometrical 
point of reflection varied from being on the crest of a ridge to being in the trough 
between two ridges. The somewhat unexpected result was obtained that the 
reflection coefficient remained substantially constant, at each angle of incidence, 
provided that the mirrors were set each time to obtain a signal maximum. The 
variable strength of reflected signal obtained from the constantly changing surface 
of the sea, therefore, may be partly due to the impossibility of maintaining a 


mirror, at every instant, in optimum adjustment. With ridges 120cm. apart, 


running along the direction of transmission, it was found, however, that the re- 
flected signal with the geometrical point of reflection on the crest of a ridge was 
appreciably greater than when the point of reflection lay in a trough. 

When working on ridged ground, multiple reflections were, in general, 
observed. Thus, when the ridges were across the direction of transmission, in 
addition to a maximum of reflected signal obtained with the mirrors set at approxi- 
mately the elevations indicated by geometry, further maxima could be obtained by 
adjusting the mirrors so that one had a greater and one a less angle of elevation than 
the geometrical. ‘Three maxima of received signal of approximately equal value, 
but having different settings of mirror elevations, were not unusual. With the 
ridges running along the direction of transmission, anumber of maxima of received 
signal could similarly be obtained by suitable adjustments of the azimuth settings 
of the mirrors. 

The experiments on ridged ground described above were made ona surface of 
comparatively low reflection coefficient, and some similar experiments were 
therefore made on a surface of greater reflecting power. For this purpose a mat of 
1-2cm. ($”’) mesh galvanized wire netting was laid between the towers, and its 


surface formed into a series of ridges. ‘The material was somewhat intractable, — 


but by pegging it to the ground at frequent intervals in the troughs, a system of 
ridges closely similar to those used for the ground experiments was obtained. 
The results obtained are given in table 10. 


The reflection coefficients were notably greater than those for similarly ridged — 


ground, but bore much the same ratio to the level values as held in the case of bare 
ground. Absorption tests showed that the wire netting transmitted 0-57 or 0-53 of 
the radiation incident normally on it, according as the electric vector was polarized 
parallel or perpendicular to the edge of the netting. This high proportion 
transmitted accounts for the comparatively low reflection coefficient with vertical 


polarization when the netting was level, and raises some doubts whether the result 


as a whole can be taken as typical for a surface of high reflecting power. 
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Table 10. Reflection coefficient of wire netting ridged at 45° to the direction 
of transmission, A=9cm. 


Angle of Vertical polarization Horizontal polarization 
Pouce D=60 D=120 D=60 D=120 
d =60 cm. = cm. =60 cm. = cm. 
cee evel | Rigi? Wh ges ice keel We pieraed MOR ts cen 
68 0-80 0-20 0-25 0-99 0-25 0:31 
Don 0-20 0:31 
43-5 0-23 0:26 


C. Other observations on rough ground 


When the investigation started, in February 1944, the site between the towers 
consisted of rough unmown land, having variations in ground level of about 
+7cm. on the mean over the area under test, and covered with a dense matted 
layer of grasses about 30cm. high, with dead stalks rising to a height of more than 
60cm. ‘The reflection coefficient of this surface was measured, and the grass then 
mown as short as the roughness of the ground permitted. The ground contour 
of the central 6 x 8 metres of the area between the towers was measured, and a 
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(a) Before removing turf. (6) After removing turf. 


Figure 4. Ground contour maps of site. 
Contours every 2:5 cm. above datum. Direction of transmission — 


contour map is given in figure 4(a). The reflection coefficient of this surface was 
also measured. ‘The turf was next removed, and in doing this some of the ground 
irregularities were removed. A contour map of the site in this condition is given 
in figure 4(b). Reflection coefficients were again measured on this surface. 

The results of these tests are given in table 11, in which are also given the 
estimated reflection coefficients for level ground of the same humidity. 


§5. DISCUSSION OF RESULTS AND CONCLUSIONS 


The investigation has provided a fairly comprehensive outline of the pheno- 
mena associated with the reflection of radiation of 9cm. wave-length, with two 
limitations. The first of these is the possibility of a scale effect when the results of 
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the experiments on diffuse reflection, made on an area of the order of 20 square 
metres, are applied to cases of transmission over considerable distances, when the 
area contributing to the reflected signal will be much larger. The second limita- 
tion is that no results have been obtained with angles of incidence greater than 80°, 
since the method used is clearly not applicable to cases where the angular separation 
of the direct and reflected rays is small. 

The general result of the experiments on diffuse reflection is to confirm the 
optical rule that, for regular reflection for a rough surface, the depth of the surface 
irregularities multiplied by the cosine of the angle of incidence must be a small 


Table 11. Reflection coefficients of rough ground, A=9cm. 


Angle.of Vertical polarization Horizontal polarization 
cneiel 
(ages) Level Long | Short ie Level Long | Short Bare 
(estimated) | grass | grass (estimated) | grass | grass 
68 0-20 0:08 | 0:20 | 0-23 0-82 0-12 | 0-74 | 0-81 
53°5 0-40 0:06 | 0-17 | 0:27 0-73 0-12 | 0-43 | 0-50 
43-5 0-48 0:06 | 0-16 | 0-26 0-68 0:03 | 0:36 | 0:46 


fraction of the wave-length. The results of tables 8 and 9 indicate that when this 
fraction is about 1/5, the reflected ray may be reduced to about one half its proper 
value, and that when the fraction is 1/2 a reflection coefficient greater than 0-1 is 
improbable. There is no reason to suppose that these relationships will be 
greatly altered by either an increase in the area of reflecting surface or in the angle of 
incidence. 

The following conclusions on the reflection and absorption of radiation of 
9 cm. wave-length can be drawn from the experiments described above :— 


(1) Except at grazing incidence, specular reflection will be unusual, since it can 
only occur at very level surfaces which are free of vegetation. 

(2) 'The reflection coefficient of bare ground when specular reflection occurs 
will usually lie between the values computed by assuming negligible conductivity 
and a dielectric constant of either 5 or of 15, according to the wetness of the surface 
layer. 

(3) ‘The change in phase on reflection is about 180° for vertical polarization at 
angles of incidence appreciably greater than Brewster’s angle, and negligible at 
angles of incidence appreciably less than Brewster’s angle. 


(4) The change in phase on reflection is about 180° for horizontal polarization 
at all angles of incidence. 


(5) For a very dry sand loam the attenuation is 36 db. per metre, the conduc- 
tivity 3 x 10% e.s.u., and the dielectric constant 2. 


(6) For the same ground when saturated with water the attenuation is 220 db 
per metre, the conductivity 6 x 10° and the dielectric constant 24. 


(7) For tap water the attenuation is 380db. per metre, the conductivity 
2x 101° and the dielectric constant between 75 and 80. 
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(8) For a 4% salt solution, approximating to sea water, the attenuation is 
1100 db. per metre, the conductivity 5-5 x 10! and the dielectric constant between 
75 and 80. 

(9) Specular reflection ceases to occur when the surface becomes uneven (as 
would be predicted by optics); and when the product of the depth of surface 
irregularity and the cosine of the angle of incidence is A/2 or more, the reflection 
coefficient is unlikely to exceed 0-1 for all polarizations and angles of incidence. 

(10) A layer of vegetation on level ground prevents specular reflection, even if 
it is only a few centimetres high; in this case the reflected signal with vertical 
polarization may be increased at high angles of incidence. 

(11) An appreciable thickness of vegetation reduces the reflection coefficient 
at all angles of incidence to a value considerable below that obtained by specular 
reflection from the ground surface; if the layer of vegetation is more than 30cm. 
high, the reflection coefficient is unlikely to exceed 0-2 at angles of incidence less 
than 80°. 

(12) The probable value of reflection coefficient for rough ground covered with 
vegetation is about 0-1, except in the neighbourhood of grazing incidence. 
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Calculated reflection coefficients are given in table 12 for surfaces consisting 
of very dry ground, very wet soil, and sea water. ‘The values for fresh water will 
only differ by about 0-03 from those for sea water. 


Table 12. Calculated reflection coefficients, A=9cm. 


Very dry soil Very wet soil Sea water 
Angle of © |——@-ANAO@_.s xX\H ie OMe 
incidence Vertical | Horizontal | Vertical | Horizontal | Vertical | Horizontal 
(degrees) | polariza- | polariza- | polariza- | polariza- | polariza- | polariza- 
tion tion tion tion tion tion 
0 +0-17 —0:17 +0:67 —0-67 +0-81 —0-81 
20 +0-16 —0:19 +0:65 —0-68 + 0-80 —0-82 
40 +0:10 —0:24 +0:60 —0:73 +0:-76 —0°85 
60 —0-:06 —0:38 +0-42 —(-82 + (0-66 —0-:90 
70 —0-21 —0:51 +0:27 —()-87 +0:54 —0:93 
80 —0-49 —0-71 —0-08 — 0-93 +0:26 —0-:97 
85 —0:71 —0-84 —0:39 —(0-96 —0-11 —0:98 


In the table a + sign indicates that the phase change on reflection is very small, 
and a — sign indicates a phase change of about 180°. 


19-2 
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ABSTRACT. The behaviour of luminescent solids under conditions of periodic excitation 
has been discussed theoretically and predictions made on the basis of simple models for the 
electronic processes accompanying the emission of luminescent radiation. The significant 
features of the luminescence intensity as a function of time are the shift in phase with respect 
to the exciting radiation and the ratio of maximum to minimum emitted intensities. The 
variation of these parameters with intensity and period of excitation can be used to distin- 
guish between various mechanisms for the luminescence process. Experiments are 
described, using a cathode-ray-oscillograph technique, illustrating the application of the 
method to zinc sulphide and zinc cadmium sulphide phosphors. The general conclusions 
are that these phosphors do not possess any of the characteristics associated with the 
monomolecular mechanism, but show a semi-quantitative agreement with the requirements 
of a simple ionization-recombination mechanism. Deviations from theory are ascribed to 


the compiexity of these phosphors, which are activated by more than one type of activator 
atom. 


wi IN TROD U Ci row 


OsT present-day applications of luminescent substances involve 
excitation by periodic stimuli, under which conditions phosphorescence 
is usually a desirable property ; for example, the flicker in a.c.-operated 

fluorescent discharge lamps is reduced by the use of luminescent coatings ex- 
hibiting phosphorescence. The behaviour of phosphors when subjected to 
periodic excitation has been studied very little, and no theoretical treatment 
of the effect has been attempted. It has been observed, however, that the wave- 
form of the total light output from a fluorescent discharge lamp shows a slight 
advance in phase over the excitation wave-form (Davies, Ruff and Scott, 1942; 
Uyterhoeven and Zecher, 1938). It was considered, therefore, that an experi- 
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mental study of the wave-form of the luminescent radiation alone might give some 
indication of the laws governing the effect and eventually provide some data 
concerning the mechanism of the luminescent process. This paper describes 
a limited experimental study of the effect, and compares the results with predic- 
tions based on a theoretical treatment. The results can be considered only 
as preliminary, but they do serve to indicate the possibilities of this method, 
particularly with regard to the mechanism of the luminescent process. 

Wood (1921), Gottling (1923) and Szymanowski (1935) studied certain 
luminescent materials by a method similar in principle to that described here, 
but based on the use of the electro-optical Kerr effect for the measurement of very 
short time intervals. Both Wood and Gottling interpreted their results as 
indicating a dark interval between absorption of radiation and emission of 
luminescence, and it remained for Gaviola (1925) to show that this interval 
corresponded to the shift in phase of the luminescence intensity-time function 
with respect to that of the high-frequency spark excitation. 


§2. THE BASIC LUMINESCENCE PROCESSES AND PERIODIC 
EXCITATION 

It is to be expected that the behaviour of a phosphor under periodic excitation 
will depend markedly on the nature of the mechanism of the decay of 
luminescence. ‘his mechanism is, in general, complicated (Randall and Wilkins 
(1945 a and b), Garlick and Wilkins (1945), Lord, Rees and Wise (1946)), but can 
be broadly classified into two main types: that in which eiectrons do not leave the 
luminescent centres, and that in which they do. In this preliminary investigation 
we have therefore examined the differential equations describing the luminescent 
emission of radiation under periodic excitation in the simplest case of each of 
the two types, and have noted certain significant differences in the solutions. 
The simplest case of the first type is that which for constant excitation leads to a 
simple exponential decay law. The simplest case of the second type is that which 
for constant excitation leads to a simple bimolecular decay law. More complicated 
cases of the two types are then discussed in a semi-quantitative way in § 3. 

All periodic functions J(¢) representing the time-dependence of the intensity 
of the exciting radiation can, if necessary, be expressed as Fourier series, but it is 
evident that the simplest possible function, namely, that describing a sinusoidal 
variation, should be employed in the initial theoretical considerations, and also, 
if it is proved to be practicable, in any experimental work which attempts to 
establish the general predictions of the theory. 

In the simplest case, in which the electronic transitions responsible for the 
absorption and emission of radiation in phosphors are all internal, i.e. in the 
case leading, for constant excitation, to a simple exponential decay law, the 
differential equation of the process is 

dN _ullt) 

dt hv 
where I(t) is the energy of frequency v incident per second, » the absorption 
coefficient of the phosphor crystals, N the number of electrons in the excited 
level at any time ¢, and « the probability per unit time of a spontaneous transition 
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from the excited level to the ground state accompanied by the emission of | 
radiation. For sinusoidal excitation, /(t) has the form 


{teens 
i \ eee ae | 


where J,, is the maximum value of J(¢) occuring at t=1/27, 3/27, etc., and T | | 
the period, so that the differential equation becomes | 


1 
dN _ uly nt 
B= 52 (100s 7) -aN. | 


This is a standard differential equation, the solution of which is 


ane + zd sin me) 
N=Ce-*#+ bLin ey HL in | E aT th | 5 
ahv ahvy 1442/27? J 


C being the arbitrary constant of integration. 
The “‘life” of the excited state, 7, the time taken for the number of atomsin such | 

a state to decrease to 1/e of the original value in the absence of further excitation, | 
is equal to 1/«, and, unless the excited state is accompanied by a metastable level, | 
this can hardly be greater than 10-° sec. This means that for large times e™ | 
is a negligible term in the expression as « is 10° sec>1 or greater. "The luminescence 
intensity is then described by the expression 
drt 20. Lat) 

Nat — Hn | SETS canta 
SS the ee L 1 + 407? /a?T? i) 


Cos 


This is a periodic function having a period equal to that of the exciting radiation. 
The phase shift of this function with respect to that describing the exciting 


d(«N) 
dt 


radiation can be obtained by equating to zero and solving for 


The solution is 


on 
[pen eee ) 


TT 
or tan e= pe 
The significant feature of this result is that the phase shift € is independent of 
the intensity of the exciting radiation. The ratio of intensities at maxima and 
minima of «N, as obtained by substitution of this value of « in the original 
expression, is 
de 


1 +. a 
V a2 T? + 4a? pal 
ae ie cot* se. 


= aS SS 
V o2T? + 472 
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This also is independent of the intensity [,,, and is determined only by the period 
of the excitation and the “‘life”’ of the excited electronic state of the activator 
atom in the phosphor crystal. A plot of the phase shift against the ratio of 
maximum to minimum emitted intensities is given in figure 1. 

In the simplest case of an ionization-recombination mechanism, i.e. in the case 
leading, for constant excitation, to a simple bimolecular decay law, the differential 


his describing the process, under the periodic conditions of excitation as 
efore, is 


where NV, is the number of electrons in the conduction band at time t and B is 
the probability per unit time that an electron in the conduction band will encounter 


Phase shift of minimum (degrees) 


O 10) 20 30 40 50 60 
Ratio of maximum to minimum intensity 


Figure 1. Theoretical relationships between phase shift of minimum luminescence intensity and 
ratio of maximum to minimum intensities under sinusoidal excitation. 


1. Simple ionization-recombination mechanism. 2. Simple monomolecular mechanism. 


a positive hole on an activator atom and recombine with it with the emission of 
luminescent radiation. The other symbols in the expression are defined as in the 


previous case. 
The differential equation is not of a standard type; it is a particular form of the 


Riccati equation, which may be transformed into a linear form known as Hill’s 
equation. A numerical solution has been given by Wise (in course of publication) 


in terms of the parameter 
aye eis 


hin? ? 
and has been tabulated for values of «2 between 0:25 and 64. The solutions are 
given in the form 


@ 2nt 2nt 
n= ped: =a+b'sin lear +6” sin 2g) boven 
7 ihe ih 
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where a, b’, b’,.... and ¢’, ¢”,.... are determined by the value of «x?. It is |} 
significant that the solution in this case depends on the maximum excitingintensity || 
and hence, also, will the phase and the ratio of maximum to minimum. emitted | 
intensities. In table 1, the phase of the maxima and minima of PN? in terms oft, | 
and the ratio of maximum to minimum luminescent intensity have been recorded | 
for nine values of «2 between 0-25 and 64. These were obtained by a method of || 
successive approximations from the numerical solutions given by Wise; the 


Table 1 
re Phase of Phase of Maximum luminescent intensity 
minimum maximum Minimum luminescent intensity 
+0 —> 90° —> 270° >1 
i 65° 26’ 242° 18’ 2°47 
4 59213" BRIO BKOY 3-40 
1 502 S04 PINE GY 4-85 
2 43° 34’ 2a On 6:87 
4 38° 30’ 201° 14’ 10-0 
8 Bye oy 194° 3’ 14-2 
16 DOge22, 189° 49’ 21°3 
32 Dice Og ISO OY 35-4 
64 ae Uy 186° 30’ 63-0 
8 +0 +> 180° >8 


limiting values are also given. he theoretical relationships between the phase’ | 
shifts at the maximum and minimum and the ratio of maximum and minimum 
values of BN,2 have been plotted in figure 1. 

The form of the solution here is in marked contrast to that of the mono- 
molecular case, and the predictions embodied in these solutions can be used to 
obtain certain information simply and rapidly about phosphors whose electronic 
processes during absorption and emission of radiation are uncertain or unknown. 


§3. MORE COMPLEX SYSTEMS: ELECTRON TRAPS 


There are several ways in which the electronic processes leading to the emission 
of luminescence radiation may be complicated ; the most important of these are 
(i) the presence of more than one type of activator in the phosphor, and (ii) the 
existence of metastable levels, associated with abnormalities in the crystal, which 
are capable of trapping electrons. The differential equations describing such 
systems are more complex: for the case of constant excitation intensity the 
theoretical treatment has been developed to a stage where the mechanism proposed 
can be said to describe the behaviour of certain phosphors in a qualitative way 
(Lord, Rees and Wise, 1946). It is evident that many phosphors considered 
to contain only one activator type did, in fact, contain more than one, and it is 
maintained that the failure to fit the simple exponential and bimolecular laws to 
experimental observations resulted, in part at least, from such a simplifying 
assumption. 

Neglecting the presence of electron traps for the present, we may examine the 
case of a phosphor, containing more than one activator type, whose absorption 
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and emission processes are monomolecular in nature. We note that these 
processes at each activator atom will be independent and that the solution of an 
equation of the type 


dN, 
dt 


al ae 
= 5-I(t)-4,N, 


is involved for each type of centre. The total luminous intensity is then 


— tml J, eT? 2nt | 2a . Int 
ee a BY {1 ee fe ceerene Lote 


The significant features of this function, namely, the phase and the ratio of 
maximum to minimum intensities, are, as in the case of a single type of centre, 
independent of the intensity of the exciting radiation. 

For the ionization-recombination mechanism the solution of a set of simul- 
taneous differential equations is required : 


dt Ww 
n 
where JN, is the number of positive holes of the vthkind inthe crystal, N= = N, 
r=1 

the total number of electrons in the conduction band at time ¢, 8, the probability 
per unit time that an electron will recombine with a positive hole of the rth kind, 
pv, the absorption coefficient for centres of the vth kind.” The solution of this set 
of equations has not yet been attempted, but it is certain that the phase shift and 
ratio of maximum to minimum emitted intensities will depend on the intensity of 
the excitation. 

A complete treatment of the problem should also take into account the 
variation with depth of the conditions in the phosphor (Lord, Rees and Wise, 
1946), but the solutions of the appropriate equations have not yet been obtained. 
It is not, however, anticipated that the qualitative conclusions will be modified. 

The existence of ‘electron traps”? in many crystal phosphors has been 
established with certainty (Randall and Wilkins, 1945 a and b), and it would 
be instructive to examine the way in which these traps influence the luminescence 
processes. For the case of constant excitation intensity in electronic processes 
involving the conduction band, some progress has been made towards a solution 
of the problem (Lord, Rees and Wise, 1946), but for periodic excitation conditions 
it is somewhat more complicated. However, when the period 7 is short, 10 sec. 
or less, it is expected that electron traps will play a relatively small part in the 
luminescence processes, and that they will not contribute appreciably to the 
emission, except in those systems where the electron must pass through the meta- 
stable or trapping level before reaching the level from which.a transition may be 
made to the ground state with the emission of radiation. 

Systems in which the electron is excited from atoms of the crystal matrix, 
followed by a migration of the positive hole to the activator atom, where the electron- 
positive hole recombination takes place, have not been considered here, as the case 
of constant excitation intensity for such phosphors has not yet been treated 


theoretically. 


286 Mary P. Lord and A. L, G. Rees | 

It therefore appears that it should be possible to determine the general type 
of luminescence mechanism followed by a phosphor by examination of ‘its | 
behaviour under short-period excitation. The derivation of such information i 
is illustrated in the subsequent sections of this paper. 


§4. MEASUREMENT OF TIME DEPENDENCE OF LUMINESCENCE 
INTENSITY 

The phosphors employed in these experiments were ; 

(i) ZnS, activated by 7-5 x 10-8 %, Ag, but also containing interstitial Zn and | 
a trace of Cu (<1x10-4°%). The phosphor therefore contained three types | 
of activating centre (see Lord, Rees and Wise, 1946). 

(ii) ZnCdS containing 18%, ZnS, activated by 1 x 10-? % Cu, and presumably | 
containing interstitial Zn and Cd also. 

Yt was found that the intensity of the group of lines at A3650 from the quartz | 
inner tube of an 80-watt MB/V high-pressure mercury-vapour lamp operated on 
50 cycles varied in a substantially sinusoidal manner with time. Such a lamp | 
constituted a convenient source of exciting radiation and was used for the experi- 
ments on the zinc and zinc-cadmium sulphide phosphors. 

For recording the luminescence intensity-time relation, a cathode-ray-oscillo- | 
graph method was employed. Light from the source, enclosed in a ventilated 
lamp-house, was passed through an 18A Wratten filter to isolate the A3650 | 
radiation and then on to the phosphor, coated uniformly on a thin glass plate. | 
Owing to the proximity of the lamp, the temperature of the phosphor was about | 
30°c. The transmitted radiation consisted of unabsorbed 13650 and emitted — 
luminescence, and the former must be removed by an appropriate filter before the 
radiation strikes the photosensitive surface of the detector. Anelectron-multiplier | 
photocell (Cinema ‘Television, Ltd., MA 20, No. 300) was used as detector, the 
output from which was passed through a 1:2 megohm resistor and the p.d. | 
appearing across this resistor applied directly to the Y-plates of a Mullard GM 3156 
cathode-ray oscillograph. Synchronization is automatic if the time base and light 
source are both actuated by the same a.c. supply. A biasing voltage of +48 
volts was required to bring the zero line to a convenient part of the oscillograph 
screen. ‘The H.'1. supply (ca. 1000 volts) for the multiplier was obtained from | 
a rotary D.C. transformer operating on 14v. D.c.; the generator output was 
smoothed by means of conventional filters before application to the multiplier 
potentiometer. Rigorous precautions were taken to exclude stray fields from the 
oscillograph. 

The light filter system between the phosphor and the multiplier receiver must 
allow only that emitted light which it is proposed to investigate to fall on the 
photosensitive surface. For ZnS.Ag, whose emission is predominantly in the 
blue, a Wratten No. 39 (blue) filter, to remove green phosphorescence due to Cu, 
and a glass cell containing saturated aqueous sodium nitrate, to remove transmitted 
3650, were employed; for ZnCdS .Cu, whose emission maximum lies in the 
orange-red, a Wratten No. 8 (yellow) filter sufficed to remove both transmitted 
X 3650 and visible blue or green radiation. 

Calibration of the vertical displacements in the oscillograph of the zero line 
and the trace with respect to exciting intensity proved necessary, since the zero 
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line was found to rise owing to a build-up effect in the 1-2 MQ resistor, and the 
responses of the multiplier and oscillograph were not linear. Calibrations were 
performed with the mercury lines lying in the range of the spectral emission of 
the phosphor under investigation by reducing the intensity in known ratios with 
a calibrated neutral wedge. 

‘The wave-forms of the luminescence for four different excitation intensities 
were photographed on the same plate together with the wave form of the Hg 
A3650 radiation, used to fix the origin of J(¢). For each plate, the deflection 
caused by light reaching the photocell with the luminescent layer removed, but 
with filters in place, was photographed; however, this deflection proved to be 
negligible. The true zero line appropriate to each observed wave form was 
derived, the height above it of the maxima and minima obtained and corrected for 
non-linear response by means of calibration curves. All traces were recorded 
independently on three occasions, the traces measured from optical enlargements 
and the mean values taken. A typical photographic record is shown in figure 2. 
The absolute intensities of excitation radiation were obtained from measurements 
with a Moll thermopile. 


Figure 2. Reproduction of a typical oscillograph record showing phase shift. 


§5. RESULTS AND DISCUSSION 


The results of the experiments conducted are recorded in table 2 and plotted 
together with the theoretical curve for the region of the experimental observations, 
in figure 2. 


Table 2 


I Ratio of maximum Phase shift of 


Phosphor (watts/cm?) to aun Ve a capers Wee 

ZnS . Ag : 6 ie 
48-9 
Sul o7/ 


ZnCds . Cu : : 223 
: : 24-4 

27°8 

30:7 


Owing to the presence of a certain amount of distortion of the peaks near the top 
of the oscillograph screen, the phases were measured only at the minima. 
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From the curves in figure 3 it is apparent that the qualitative behaviour of the — 
phosphors is as predicted by the theory outlined in $3, but the quantitative |) 
agreement is poor. This is not suprising in view of the thickness of the layer, 


the complex nature of the systems constituting these phosphors and the possible i 


effect of electrons trap at the relatively low frequency of excitation used (10? cycles | | 
per second). However, the results demonstrate quite conclusively that the || 
luminescence mechanism is not monomolecular in type since all points in figure 3 
for one phosphor corresponding to different intensities would then be identical, 
and, indeed, lend support to the belief that it is, ca. 30° c., of an ionization- | 
recombination type. It is not necessarily expected that the deviations from — 
theory for the two phosphors investigated should be the same; however, it can be | 
seen from figure 3 that both sets of points approximate to the same curve. As the | 


100 
80 


60 


THEORETICAL 
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40 


Phase shift of minimum (degrees) 


20 


O @ 4 6 8 lO 


Ratio of maximum to minimum intensity 


Figure 3. Comparison of theory and experiment. Theoretical curve from table 2 for ionization- | 
recombination model. Experimental observations: © ZnCdS.Cu, @ ZnS. Ag. 


deviations are in all probability largely due to interstitial lattice atoms (Zn or Cd), 
whose equilibrium concentrations at the temperature of preparation of the 
phosphor would be of the same order, this fact is not surprising. It is evident that 
all phosphors following the simple bimolecular mechanism should satisfy the 
theoretical relation of figure 1 irrespective of the value of the recombination 
coefficient, absorption coefficient, intensity or periodicity of the excitation. 

It would appear that the principles of the method here discussed should be of 
considerable value in distinguishing between various models of the luminescence 
mechanism, since the interpretation of decay curves under constant excitation 
conditions is a matter of some difficulty and uncertainty. 
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ABSTRACT. For the provision of rapid and convenient means for the examination 
of luminescent solids by oscillographic methods, the possibility of producing square-wave 
light pulses from a gaseous discharge has been investigated. Decay curves sufficiently 
accurate for the estimation of half-life periods and for the study of gross effects could be 
obtained by interrupting the D.c. supply to the discharge by means of a small vibrator. 


luminescence of a phosphor is its “ half-life’’, that is, the time taken for the 

emitted intensity to fall to half its value after the removal of excitation. 
Methods less tedious and precise than the orthodox have therefore been developed 
to derive decay curves rapidly. 

An electron-multiplier phototube together with a cathode-ray oscillograph was 
first used for the plotting of decay curves by de Groot (1939). His method 
involved the use of rotating sectors for the production of square pulses of exciting 
radiation; the apparatus described here provides for the production of these 
square pulses directly from the gaseous discharge. ‘The specific advantage of this 
method lies in the absence of motor-driven moving parts, leading to a more 
convenient and portable equipment; it also avoids the provision of an adjustable 
constant-speed motor and provides simpler synchronization for the oscillograph 
time-base. Use was made of the same multiplier-oscillograph combination as that 
employed in the study of the behaviour of sulphide phosphors under periodic 
excitation (Lord and Rees, 1946). We were concerned here, however, with the 
examination of silicate materials and, therefore, used a low-pressure mercury 
discharge as the source of A 2537 radiation. 


F« comparative and routine purposes, a sufficient criterion of the decay of 
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Square wave-forms were produced by the following methods : (i) Theuse ofa 
bright-emitter diode run under saturation conditions ; (ii) the use of a thyratron 
in series with a saturated diode; (iii) interruption of the Dio supply to the dis- 
charge with a mechanical vibrator. ‘The first method did not give a completely. 
satisfactory flat-topped wave and could not be used as such; when, however, a 


Low PFESSHE 


MIELeEUry Ascharge 2 


JS LIP, FOTSING 
ete wy pulse fo 
L, CR. Sube. 
{ 
Weshy7 
Vibrofor 
Figure 1. Circuit for production of square-wave pulses for A 2537 radiation. 


thyratron was connected in series, the wave-form was improved and variation of 
the pulse duration by control of the grid bias then became possible. 

As the current requirements of the discharge lamp were of the order of 70 ma. 
for satisfactory light output, the use of a mechanical vibrator (V 6121, 12 volt) was 
found to be simpler. ‘This vibrator was used in the circuit outlined in figure1 to 
interrupt the 260-volt D.c. mains, which were smoothed by means of a conven- 


() 


(a) 
(6) (a) 
(€) 
Figure 2. Wave-form characteristics of vibrator circuit. 


(a) Discharge current ; (6) lamp voltage ; (c) vibrator voltage ; (d) light output from discharge ; 
(e) build-up and decay of luminescence in manganese-activated zine silicate. 


tional pi filter. A synchronizing pulse was taken to the oscillograph through a 
condenser and resistor as shown. ‘The 6-megohm series resistor ensured that a 
constant small number of ions were available to facilitate starting up the lamp on 
each cycle. For currents up to 70 ma. (the lamp rating), satisfactory wave-forms 
were obtained, but with some build-up and decay and a small-amplitude high- 
frequency oscillation superimposed on the flat top. The high frequency has its 
origin in the discharge, as can be seen from the wave-forms in figure 2, (a, b and c) 
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The current wave-form of figure 2(a) was obtained by placing aseries resistance of 
150 ohms in the negative line and applying the voltage established across this 
resistor to the vertical deflector plates. Figure 2(b) was obtained by applying the 
voltage across the lamp directly to the deflector plates, using a shunt of 250,000 ohms 
across thelamp. ‘The vibrator voltage wave-form, shown in figure 2(c), was obtained 
by connecting the vertical deflector plates across 50,000 ohms of the 250,000-ohm 
shunt; however, the discharge was suppressed in this case by removing the 
filament heating supply. Bounce of the contacts on make is evident at the begin- 
ning of the trace. 

Curves showing the build-up and decay of luminescence, as in figure 2 (e), are 
obtained by interposing the phosphor layer and appropriate filter system between 
the square-wave source and the photomultiplier. Half-life periods deduced from 
such curves are in substantial agreement with those obtained by more precise 
methods (Lord, Rees and Wise, 1946), provided corrections are made from cali- 
bration for the non-linearity of the multiplier-cathode-ray oscillograph system. 

The method may be extended to a wide range of pulse durations and repetition 
rates by substituting vibrators of different frequencies and times of contact. The 
vibrators do not perform satisfactorily if maximum current is drawn for long 
periods; the wave form deteriorates if the unit becomes heated. 

Although neither this method nor that of de Groot can be used to obtain decay 
curves sufficiently precise to give useful data about the electronic processes 
accompanying luminescence, this method is quite adequate for the study of gross 
effects and for the rapid comparison of phosphors, 
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ABSTRACT. Several new band systems of N, which have been reported by us indepen- | 
dently during the war are compared. Several systems have been obtained by both of us, 
and agreement is generally satisfactory. ‘There appear to be eight progressions of bands 
in the ultra-violet which are due to transitions to the afl state. Differing vibrational 
analyses for the Fifth Positive (Van der Ziel) system are discussed, and it is shown that 
Gaydon’s interpretation is to be preferred, while for Kaplan’s Second system (Herman’s a) 
Herman’s vibrational analysis is more satisfactory. Bands in the green region are compared. 


§1. INTRODUCTION 
LD: RING the war we have independently investigated band spectra obtained 


in emission from nitrogen in discharge tubes and have each found several 
new systems (Gaydon, 1944 a, 1944b ; Gaydonand Worley, 1944; Herman, 
1945). Onthe whole our results are in fairly good agreement, but there are some 
ways in which our interpretations of the observations differ. ‘The names which 
we have applied to the various new systems are naturally quite different. Thus 
a joint publication seems desirable in order to assist others to sort out the 
numerous systems attributed to N, which have now been found and, as far as 
possible, to clear up some of the differing interpretations of the observations. 
Generally speaking, Herman has developed rather more intense sources, and 
has thus been able to obtain some of the band systems more completely, so enabling 
the more extensive and reliable vibrational analyses. Gaydon, on the other hand, 
has used spectrographs of larger dispersion and has concentrated more on the 
rotational structure of the bands; data obtained from rotational analyses of a 
number of bands have led to definite conclusions on the nature of the electronic 
levels. 
The characteristics of the discharge-tube sources used for exciting the spectra 
have been rather different, and it is not surprising to find that while some of the 


new band systems have been observed by both of us, other systems have been 
reported by only one of us. 


§2. BAND SYSTEMS INVOLVING TRANSITIONS TO THE 
Ga SeATSE 
Many of the new systems consist of single progressions of bands showing the 
same vibrational intervals as those of the a1II state, which is the upper level of the 
main Lyman system of N,. This level was formerly believed to be 1II,, but 
Herzberg (1946) has recently suggested that it may be 4II,. These band systems 
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thus consist of transitions from a number of highly excited levels to the a1II level, 
this being confirmed in several cases by the rotational analyses. In table 1 we 
have set out brief particulars of these systems, with the names which we have given 
to them in out separate papers. 


Table 1. Band systems with a'II as final level 


Type 
Approx. A Gaydon’s Herman’s of upper 
Vo (0,0) (0,1) (0,2) designation designation electronic 
state 
S57 2827 2967 3119 P C us 
S5765 <6 12795 2932 3079 — Ui () 
36394 2746 2878 3020 Q onal 17] 
37417 2672 2796 — R Dy 
41705 2397 2497 2603 5 — 15} 
A343 2359 2455 2558 — t) 
43818 2281 2371 —_ al € 1} 
46611* — 2224 2308 _ n 


* These values are derived from head measurements ; other values of ») are for band 
origins. 

If the a'I[ state is ‘lly, then the '» states will be ‘7, as given in Gaydon’s paper, and the 
upper HI state will be tly. If, however, the a'II state is 1I;, then the 1% states will all 


be 15+ and the upper 1II will be 14II,,. 


For the system with vy =35371 our analyses agree; these bands appear more 
strongly on Gaydon’s plates. 

The two progressions with vg 35761 and 36394 were tentatively grouped 
together as one system by Herman. However, since the first progression is not 
present on Gaydon’s plates, it now seems more probable that the two progressions 
belong to separate systems. The value of w,;, about 640 cm-", obtained by combin- 
ing the two progressions, also seems rather small for a state with a B value of about 
1-36. Herman obtained four bands in each progression. Gaydon obtained only 
the three strongest bands of the vp 36394 progression, but achieved rotational 
analyses of two of these. Janin (1943) has independently observed bands of both 
of these progressions, his measurements agreeing with ours. 

The 37417 progression was only reported by Gaydon, who made a rotational 
analysis of the (0,0) band and also found the (0,1) band, which, it now appears, 
very nearly coincides with the (0,0) band of the 35761 progression. ‘The (0,0) 
band at 42672 is, however, present on Herman’s plates. 

The v)=41705 progression again was only reported by Gaydon, but some of 
these bands are also present on Herman’s plates. 

The 42373 progression was reported only by Herman, but there is some weak 
structure in the region of these bands on Gaydon’s plates as well. 

The 43818 progression has been obtained by both of us. Gaydon discussed 
the possibility that this might be the v’=1 progression of the system with 
¥y= 41705, but decided the evidence was not sufficiently conclusive. 

The progression with vg at 46611 was reported only by Herman, but the not 
(0,1) and (0,2) bands at A 2224 and 2308 are present on Gaydon’s plates, although 
higher members of the progression are absent or masked. 
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Thus there is evidence for the existence of eight independent systems of N, in 
the ultra-violet all having the common final level aI. Most of these systems have 
been observed by both of us, although our original publications only gave three 
systems 1n common. 

The close numerical coincidence between the energies of the upper levels of 
Herman’s « and ¢ systems with those of the levels # and c obtained from absorption 
measurements by Watson and Koontz and by Birge and Hopfield led her to 
think that these levels, and c, were the upper levels of the « and ¢ systems, thus 
indicating a violation of the selection rule that the symmetry u or g must change. 
Herzberg (1946) has now suggested that the a 1II state is of g symmetry and that the 
forbidden transition is the Lyman system itself, a4IIy=x1Xj. This accounts for 
the great path length of gas required to produce the Lyman bands in absorption. 
It brings the upper level of the 43818 progression into coincidence with the level 
h, as suggested by Herman. ‘The upper level of the 36394 progression would then 
also fit with Worley’s level m, but the upper level of the 35371 progression could 
only be identified with the level c by assuming that the electronic type of this state 
c is 1X} and not 4II,, as formerly supposed; the B values agree satisfactorily. 
There is, therefore, much to be said for this suggestion by Herzberg, but it should 
be pointed out that it does not explain why strong absorption bands to the upper 
levels of our other progressions have not been reported. 


§3. THE FIFTH POSITIVE OR VAN DER ZIEL’S SYSTEM AND 
KAPLAN’S SYSTEMS 

Three other band systems of N, with related energy levels have also been 
studied by both of us. Our designations for these systems are given in table 2. 

Van der Ziel originally reported two progressions of bands separated by an 
interval of 1827 cm-1, and he assumed that these were the v’ = 0 and 1 progressions 
of a system. ‘The intensity distribution in this system on this assumption was 
obviously unsatisfactory. Gaydon obtained these two progressions on_ his 
plates and also found an additional progression separated by 1869 cm-! from the 


‘Table 2 
Gaydon’s Herman’s 
designation designation 
Fifth Positive (van der Ziel) Van der Ziel 
Kaplan’s First System B 
Kaplan’s Second System a 


first of the older progressions. By treating the new progression as the v’ =0 and 
the original progressions as v'=1 and 2, a satisfactory intensity distribution, in 
conformity with the Franck-Condon principle, was obtained. Rotational 
analysis of one band of the new progression proved that it belonged to the same 
system as the older progressions. 

Van der Ziel’s rotational analysis had shown that the bands were due to a 
transition between two & states, and Herman made an attempt to identify the upper 
state with that of a system of bands observed by Collins and Price in absorption, 
it being assumed that van der Ziel’s two progressions arose from levels with 
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v=1 and 3 of this state. By assuming preferential excitation to these levels a 
reasonable intensity distribution could be obtained. However, it is now found 
that Gaydon’s additional progression is also present on Herman’s plates, and as 
this cannot be fitted into Herman’s scheme it seems that Gaydon’s interpretation 
is to be preferred. Moreover, Gaydon points out that Herman’s scheme would be 
inconsistent with the strong degradation of the bands to shorter wave-lengths, and 
that it is now known that the bands are due to a transition between two 1» levels, 
whereas the upper level of Collins and Price’s system is believed to be a triplet. 

Kaplan originally reported a system with (0,0) band at 21534. and gave 
measures of six bands to the nearest angstrom, and suggested that they might 
involve a transition to the same final level as that involved in van der Ziel’s system. 
Gaydon observed only three bands clearly, with a suggestion of a fourth, but 
was able to make an analysis of the (0,2) band, showing that the transition was 
1[[—X and obtaining the value of B’, and confirming that B’’ had the same value 
as for the Fifth Positive system. Herman obtained a large number of bands 
(seven in all) and was able to show that w;’=1708cm7=". , 

Kaplan’s second system of bands, which he also believed to have the same final 
level as van der Ziel’s, has been observed by us. Gaydon only obtained the bands 
weakly and was unable to make a rotational analysis, and said that the vibrational 
intervals agreed as well as could be expected with those of the lower state of the 
Fifth Positive system. Herman, however, obtained this system more strongly, and 
recorded eleven bands, including new bands at AA 2354 and 2263, which are also 
subsequently found to be present quite clearly on Gaydon’s plates and which prove 
that the transition does not involve the lower state of van der Ziel’s system, as had 
been thought by Kaplan and Gaydon. ‘The vibrational energy interval for the 
upper state is about 1709cm-', which suggests that this system has the same 
initial II level as Kaplan’s first system. 


§4. OTHER SYSTEMS 


For two other systems our respective designations are given in table 3. 

Herman’s y system consists of two fairly long progressions, and the analysis 
shows that the transition is to the A®& level of N,. Bands of this system are not 
present on Gaydon’s plates. Kaplan had previously observed a number of these 


Table 3 
‘ : Herman’s 
Gaydon’s designation decade 
(Kaplan’s third system) y 
Green system x 


bands and had suggested that some of them formed a progression involving the 
lower level of van der Ziel’s system. This is now shown to be incorrect. 

A system of rather complex bands in the green region was observed by Gaydon, 
who arranged seven bands into a tentative vibrational scheme. ‘Herman also 
reported three bands at AA5594, 5330 and 5092 under small dispersion and 
referred to them as the y system. More recently Herman has obtained this system 
with greater intensity, and revised measurements of five bands show clearly that 


this y system is identical with Gaydon’s Green system, 
202 
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Mutual comparison of our results has led to the identification of a number of 
the unassigned bands on our plates, but there are still a number of weak bands, 
probably due to N,, especially at the far end of the quartz ultra-violet region, which 
have not been analysed. 
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DISCUSSION 


Dr. S. P. Stnua. I would like to know the experimental conditions under which 
the new systems of nitrogen bands discovered by Dr. A. G. Gaydon were developed and 
if it was possible to explain their appearance in the light of the conditions maintained during 


the experiment. 

AuTHoR’s reply. The experimental conditions for exciting the new band systems 
were fully described in earlier papers (Gaydon, 1944 a, 1944 b) and the conditions of excitation 
were discussed. Briefly, the ultra-violet systems were obtained by a condensed induction 
coil discharge through N, in a simple Pyrex tube, while the Green system required a 
Tesla discharge at relatively high gas pressure. 
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ABSTRACT. Using those differential transfer coefficients of a lens system which specify 
the rate of change of the tangential aberrations with the refractive indices of the component 
singlet lenses of the system, a simple method is available for the rapid calculation of the 
chromatic variation of these aberrations. In a similar manner, other transfer coefficients 
(the chromatic coefficients of the system) provide a rapid means of obtaining a complete 
analysis of the secondary spectrum of the system at any zone in the case of the axial 
aberration, and at any obliquity of the transverse aberration. Examples of an actual 
system are given. 


Sie IN TRODUCGRLON 


“Fue aberrations of an optical system are determined normally for incident 

| light of some mean wave-length in the spectral range over which the system 

is expected to give satisfactory images. The full significance of the chro- 

matic aberration of the system is really apparent only when, in addition to the usual 
analysis, the variation of the aberrations with the wave-length of the incident light 
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has been examined. Hitherto the analysis of the dependence of the aberrations on 
the wave-length of the incident light has required repeated ray-tracings in different 
colours. It is of interest, therefore, to examine this matter in the light of the pro- 
perties of the transfer coefficients which have been developed in the earlier papers 
of this series. 

Suppose that the aberrations of the system for incident light of some mean 
wave-length, d, have been obtained by a ray-trace and the general transfer coeffi- 
cients have been calculated from the trace. If the mean refractive indices of the 
glasses of the components a, 6, ...h, ... of the system are changed by amounts 
5N,, 5N;y,...5N,,..., the change in the value of the aberration A’, which may be 
any of the tangential aberrations of the system, resulting from the glass changes is 
given by 

(aytee 
0AF= 2 aN, GIN 

If, however, instead of considering glass changes we consider a change in the 
wave-length of the light which is incident on the system from the object, then the 
value of the aberration, A’, for light of wave-length, 7, will be given by 

Ay= Ay +E St Na) 
Hence, using the computed values of the transfer coefficients, 0A’/ON,, for the 
various aberrations, we can investigate the chromatic variations of the aberrations 
quite quickly. 

§2. A TYPICAL EXAMPLE 


Asan example we consider a rough design of a wide-angle photographic objective 
of the Ross Xpres type. In table 1 the computed values of the transfer coeffi- 
cients of the aberrations with respect to the refractive indices of the components 
are shown, the oblique aberrations being calculated from a pencil of about 22° 
obliquity. 

Table 1 
LAp LA, Xin Coma Dist’ 
0/0Na 150-13 52302 —100:48 0-891 7°780 
a/@Np —140°51 —49-202 1 15°59 0-167 —7 029 


0/ONe ail 1:960 — 23:54 —1:887 1-382 
0/0Na Paes) D220 8:48 1195 Say 
a/OoNe —115°64 —42°750 47-16 —0-651 Sy 


0/ON¢ 117-16 A300" =—="35°34 —0:240 . = 6-180 


The glass selected for components a and f was DBC 615553, for components 
b and e LF 549467, and for components c and d HC 519604. ‘The dispersions of 
these glasses are given intable2. The ray-trace of the system having been carried 
out for sodium light, we calculate next the aberration differences, A’— A}, for 
wave-lengths from 7065 a. to 4341 a. Since the glasses in the front and rear halves 
of the lens are similar, we may combine the coefficients for the components which 
have similar glasses, as, in the evaluation of the differences by the equation given 
above, these coefficients will each be multiplied by the same dispersion values, 
In table 3 the full calculation is set out for the case of the marginal spherical 
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aberration, and finally, in table 4 the chromatic differences for all the aberrations 
are collected, the detailed computations being omitted. | 

‘The accompanying graphs in figures 1 to 3 show the variation of the aberrations | 
plotted against the wave-length, and thus completely summarize the analysis of the 


Wave lengths in AU 
Wave lengths in A.U. 


LA’—LA‘y 
Figure 1 


Comayp—Comafpp (22° pencil) 


Figure 2. 
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method. ‘The results are given in table 5, the discrepancies being little greater 
than the uncertainty of the last figure of the traces. 


Table 2 
DBC DE eke 
Np-Np —0:00517 —0:00539 —0-00406 
Nco-Np —0:00328 —0:-00342 —0-:00256 
N.-Np 0-00274 0-00291 0:00211 


Ny-Np 0:00785 0:00835 0-00603 
Ne-Np 0:01422 0-01529 0:01089 


Table 3 
Contrbn. Contrbn. Contrbn. Contrbn. Contrbn. 
Comp. Coeff. b_D C_D eD F_D G_D 
a, f 267-29 —1-3819 —0:8767 0:7324 2:0982 3-8008 
b,e —256:15 1:3806 0:8760 —0°7454 —2°1388 —3:9164 
c, d 19-64 —0-:0798 —0-0503 0:0415 0:1185 0:2140 
LA’—LA>y = —0-0811 —0-0510 0:0285 0:0779 0-0984 
Table 4 
Range LAy EA Xin Coma‘, Dist ’ 
b-D —0-0811 —0-0263 —0:0450 0-00205 —0-00043 
C-D —0-0510 —0-0165 —0-:0291 0:00129 —0-:00027 
e-D 0:0285 0:0087 0:0339 —0-:00109 0-00015 
F-D 0:0779 0:0236 0:0996 —0-00310 0:00040 
‘-D 0:0984 0:0275 0:2082 —0-00568 0:00047 
Table 5 
Marginal Zonal 
Trace Coeffs. Trace Coefts. 
LAj-LAp —0-066 —0:051 —0:018 —0-0165 
LA;-LA} 0:071 0:0779 0:022 0:0236 


§3. THE ANALYSIS OF THE SECONDARY SPECTRUM 


To complete the discussion of the chromatic variation of the aberrations, the 
matter of the secondary spectrum should be considered. It has been shown 
previously that the final intersection length, L/., for an axial ray of some colour 7 
is related to the corresponding quantity, Lj, for the traced ray for some mean 
wave-length, d, by 

l= 3 be wy — Nj) 
iP d ON, r d/h 


and for a principal ray of an oblique pencil 
0H; 


eae 


CNT ING) a 
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where H’ measures the height at which the principal ray intersects a fixed 
plane, conveniently the paraxial image plane. These transfer coefficients, 
aL'/AN,, and @H'/AN,,, which may be termed appropriately the chromatic coeffi- 
cients of the system, provide a means of analysing the importance of the secondary 
spectrum as it affects the axial and transverse chromatic aberrations of the system. 
It will also be obvious that they furnish a simple means of determining the effect 
on the colour correction of changes in the glasses of the system provided these are 
not too drastic. 

The continuation of our example will afford a clear insight into the use of these 
chromatic coefficients. ‘In table 6 are collected the computed values of the chro- 
matic coefficients for the marginal, zonal, and paraxial rays of the axial pencil, and 
the principal ray of an oblique pencil at 22°. 


Table 6 
Comp. OL;,,/ ON 0L;/0N al’/ oN 0H’/ aN 
a —737:99 — 640-16 — 587-86 — 39-917 
b 768:53 677-23 628-03 30-309 
c —162:67 —157°51 —155-56 — 6:300 
d —127-07 —119-76 —114-53 8-881 
e 574:99 502°11 459-53 — 30-148 
ie — 521-13 —446-98 —403-98 38°215 
Table 7 
C Chr: Contrbn. Contrbn. Contrbn. Contrbn. Contrbn. 
ORRR-* i Cocté: b-D ep e-D -FD GoD 
Ey ie — 1087-14 56205 3-5658 —2-9788 —8-5340 —15-4591 
b,e 1179-33 —6-3566 —4-0333 3-4319 9-8478 18-0320 
c,d — 277°27 1:1257 0:7098 —0-5850 —1-:6719 — 3:0195 


Let 0:3896 0:2423 —0-1319) =—0°3585 — 0:4466 


In table 7 the contributions made by the pairs of components of similar glass 
to the zonal chromatic aberration, L;— L/y, are computed for the various spectral 
ranges, the last line of the table giving the final value of the aberration. Omitting 
the corresponding detailed computation for the other three rays under consider- 


Table 8 
Aberration b-D C-D e-D F-D G’-D 
Liy—LinD 0:4443 0:2767 —0-1518 —0:4128 —0:5176 
Li-Lin 03896 0:2423 —0-1319 —0-°3585 —0-4460 
l’-lp 0:3634 0:2258 —0-1232 —0-3349 —0-4192 
’-Hy —0:00255 —0-00158 0:00126 0:00354 0:00638 


ation, we collect in table 8 the values of the aberrations calculated for each ray. 
This provides a sampling of the axial colour correction in the paraxial, zonal, and 
marginal regions, and a single sampling of the transverse colour at one obliquity. 

The curves marked aa in figures 4 to 7 show these aberration values plotted 
against the wave-length, and hence summarize completely the condition of the 
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system as regards its secondary spectrum. ‘The axial chromatic aberration is very 
undesirable and calls for attention. A stu dy of the chromatic coefficients and the 
contributions to the aberrations shows that an improvement should be effected by 
increasing the V-number of the DBC components and also the HC components. 


Wave length in A.U. 


Bs 


aa 
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Wave length in A.U, 
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Figure 6. 


Figure 7. 


By way of example we consider the result of changing to a DBC of the type 613568, 
and also substitute a BSC of the type 517641 for the HC of components c and d. 
Ten minutes’ computation with the new dispersion values in place of the former 
ones gives a new set of aberration values. Omitting the details, these are collected 
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in table 9, from which the curves bb in figures 4 to 7 have been drawn. These ‘| 
curves show a very great improvement in the axial colour correction and a slight 
improvement in the already good transverse aberration. This example affords || 
clear evidence of the general usefulness of these chromatic coefficients of the | 


system. 
Table 9 
Aberration b-D C-D e-D F-D ‘-D 
Lig Li 0:1888 0:1073 —0-0250 —0-0160 0:2198 
LL) 0:1627 0:0920 —0:0199 —0-0059 0:2097 
Ut, 0:1525 0-0861 —0:0196 —0-0069 0:1920 
HH, —0-:00221 —0:00136 0:00116 0:00298 0:00523 


As a final check on the method we investigate the reliability of the curves | 


aa as an analysis of the secondary spectrum. Ray-tracings were made for the three 
axial rays in colours C and F and the final intersection values for each compared 
with the values deduced by using the chromatic coefficients. 


Table 10 
Marginal Zonal Paraxial 
Wave: Trace Coefts. Trace Coeffts. ‘Trace Coefts. 
length 
D 108-445 — 106°522 — 106°525 —- 
C 108-734 108-721 106:766 106-764 106:748 106°751 
F 108-030 108-032 106:159 106:163 106°181 106-190 


The first line in table 10 gives the intersection lengths for D light for the three } 


rays in the original trace. From these values and the aberration quantities in 
table 8, the predicted values of the intersection lengths of the rays for C and F light 
are obtained and entered in the table. The values obtained from the complete 
traces are entered alongside the predicted values, from which it is seen that the 
agreement is within the limits of uncertanty of the traced values. 
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ABSTRACT. Bose-Einstein statistics are applied to Mendelssohn’s 2-particles and shown 


to give a satisfactory qualitative explanation of the phenomena of liquid helium. 


§1. INTRODUCTION 


which 2-particles (either atoms or electrons) occupy a set of quantum states | 
lower than the ordinary thermal states and separated from them by a finite | 
energy gap. ‘The z-particles form an aggregate in momentum space with zero |} 


D R. MENDELSSOHN has recently propounded the concept of a z-state in | 


—— 


| 
| 
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thermal energy, and transport phenomena are caused by a diffusion of g-particles 
under their own zero-point momentum. 

As Mendelssohn remarks, this concept should aid the development of a 
satisfactory theory both of liquid helium and of superconductivity. However, he 
appears to doubt whether the Bose-Einstein statistics are capable of handling the 
situation for helium atoms, and favours the possibility that some kind of com- 
promise statistics may be found that will work alike for electrons and helium atoms 
in the z-state. The present writer believes this to be essentially erroneous, if only 
because of the fact that the two types of statistics arise from mutually inconsistent 
symmetry conditions on the wave functions representing the particles ; and it is 
difficult to see what better compromise can be found than the classical statistics, 
which obviously do not apply to either case. 

The object of the present note is to show that, in fact, Mendelssohn’s concept 
can be harmonized with the Bose-Einstein statistics when the latter is suitably 
applied to the model proposed by Mendelssohn. Previous attempts to apply the 
Bose-Einstein statistics to the liquid helium problem have all neglected the 
possibility that two competing sets of states exist, and that a statistical equilibrium 
has to be set up between them. When this equilibrium is set up, the Bose- 
Einstein statistics give a very satisfactory explanation of the phenomena in 
liquid helium. 


See AMi om CAT ACHAR AChE Ris DIGS VOR hr b7-s i) Audick 


The z-particles with zero, or at any rate with less than normal, thermal 
energy must be treated statistically as having fewer degrees of freedom than the 
three normally characterizing the thermal states. If it is provisionally assumed 
that they have but two degrees of freedom, we at once arrive at a suggestion of a 
mechanism whereby this reduction may in fact occur: because the z-state is 
certainly associated very intimately with the surface layers of liquid helium, 
especially in the frictionless transport of helium atoms in films. 

Without, however, attempting to specify the mechanism, let us accept the 
hypothesis that the z-particles have two degrees of freedom only, and that as a crude 
approximation we may regard the energy states of the z-particles as determined 
by Schrédinger waves in resonance with the surface of the assembly. This 
concept agrees qualitatively with Mendelssohn’s requirement that the z-particles 
shall have a higher degree of order than the thermal states, and form a separate 
phase in momentum space alone. It is further necessary to assume, so as to 
coincide with the Mendelssohn model, that the energy zero of the surface states 
is at a finite level below the lowest thermal states, so that there will be no 
momentum transfer from the surface or g-states at sufficiently low temperature. 


The energies of the thermal states are 
e, = h?(s? + 2? +u?)/8m,V78, s,t,u=1,2,... sees, (1) 


and of the zg-states 
é,= —yx + h2(s? + t?)/8m,A, aaa ees rye ene ete} 


where A isa crude approximation for the effective area of the assembly, and — x is 
the Mendelssohn energy gap: by hypothesis x is positive. 
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§3. STATISTICAL EQUILIBRIUM BETWEEN Z-STATES AND 
THERMAL STATES 


If there are in all N particles distributed among the two sets of states, then | 
the Bose-Einstein statistics will lead to the equilibrium given by | 


N=N,+N,, 7 a ee (3) 
ae jh?(s? + ¢? + u?) 
N,== 24 id SE 4 
BCS ° Gnas exp | 8m,kTV22 (4) 
2 yy : jh?(s? + t?) 
=> =D Te ee 
ee perio) E Xie lt neTA ) 
are the mean numbers in the thermal and g-states respectively. 
Write Ae = AeXPT SA, ee (6) 
ea) " jh?(s? ah t?) 
N, == DA*s A ee 7 
pete ae exp 8m,RTA (7) 


The expression (4) was studied by Fowler and Jones. It was summed by means | 
of a theta-function transformation and shown to approach the following limit for _ 
large N: 


N, = (22mkTP2Vh die" eS (8) 
Ka 
provided LST, SA ale 
where 2amk T,>1 G2 =(N/VPPR aaa (9) 


The corresponding summation of (5) leads to 
N, = 2amkTAh? > Aj ee (10) 
, fom 


So long as A¥ <1 the number N, < N; and can be neglected. The parameter A 
is determined then by the equation 


N = (2amkT)2V h-3 & di j3!2, 
j=1 
As T falls, A increases, and as soon as 
A> exlRT <1, A* +1, 


the series for N, will begin to diverge, so that NV, can no longer be neglected. There 
is thus a critical temperature, where 


N = (2rmkT,)2VhS e-FXRPo/ jz, (11) 
j=1 
Below this temperature A remains almost at the value e-X/"7, so (3) becomes 
Tal ? 
N= (2a mkT)2 Vi eAXR [782 4 207 mRTAR-2 & A*Ij. (12) 
geil j=1 


The population of the z-states is thus negligible for T>T » and below T ¢ it is 
given by ; 


Se ixikr [pe | 
feat 

Soe i 2 ihn, eae (13) 
=1 


j= 


Ne, ay 
Neth 


c. 
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As Mendelssohn remarks, the variation of N, is more rapid than would be given by 
Neel ATT yes 


and (13) will in fact give good agreement with experimental data provided y is 
suitably chosen. 


§4.. CONCLUDING REMARKS 


The present writer, with Dr. Hsu, formerly (1941) worked out a theory of 
liquid helium based upon the idea of surface resonance, and obtained remarkable 
agreement with observational data on specific heat and thermomechanical 
phenomena. ‘The chief weakness of that theory was that no definite explanation 
was found for the frictionless nature of the film transport, although indeed the 
concept of diffusion by zero-point momentum (there called intrinsic pressure) 
was suggested. 

Mendelssohn has now pointed out that a finite energy gap between the z-states 
and the thermal states will prevent frictional transfer, and by incorporating this 
idea, as in the present paper, the Bose-Einstein statistics can certainly be made to 
give a satisfactory qualitative explanation of the phenomena of liquid helium. 

There are two important details which the theory cannot handle in its present 
form: First, the critical temperature is in fact hardly more than half that predicted 
by (11), and second, the critical temperature follows the A-line in its dependence on 
external pressure. 

But it is almost obvious that the discrepancy here is due to molecular forces, and 
the fact that the energy states (1) and (2) are derived on the assumption that the 
container is free from potential irregularities. It is clearly quite hopeless to 
improve matters in this respect until a better theory of liquids is available, and the 
Bose-Einstein statistics applied to the helium system as a liquid rather than a gas. 
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DISCUSSION 


on paper by S. Baxter entitled “‘ The thermal conductivity of metals” 
(Proc. Phys. Soc., 58, 105). 


Mr. J. B. Carne. I am prompted to draw attention to the account of the early work 
of Sir Benjamin Thomson, which was presented to the Royal Society in 1786 e¢ seq. 

From a systematic series of beautifully simple experiments, Sir Benjamin deduced not 
only the now generally accepted explanation of the mechanism of heat transfer by con- 
vection, but also the existence of a layer of stagnant air which acts as a barrier to the 
transmission of heat from a hot surface, the postulate invariably attributed to Langmuir. 

Reference to the papers will also show that Thomson found the same variation of 
conductivity with density for very low bulk densities of air-fibre mixtures as recorded 
by Baxter. It is interesting to note also that Sir Benjamin Thomson also enquired into 
the relation between water content of air-fibre mixtures and their conductivity. 
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ABSTRACT. The “ fly’s eye”’ is the name given to a multiple camera for producing 


two-dimensional diffraction gratings consisting of repeated patterns representing projections | 


of crystal structures. This paper describes the construction of a fly’s eye using lenses 


embossed on Perspex ; also the necessary adjustments, the method of photographing | 


a structure, and the apparatus for viewing the spectra produced by the two-dimensional 
gratings so obtained. 


§1. INTRODUCTION 


HE purpose of the “‘fly’s eye”’ is to produce on a photographic plate a regularly 
repeated pattern representing the projection along some crystallographic 
direction of acrysta] structure. This behaves as a two-dimensional diffrac- 
tion grating and gives spectra whose intensities are proportional to the intensities 
of x-ray reflections in a zone corresponding to the direction of projection. The use 
of the fly’s eye thus obviates the calculation of structure amplitudes of a proposed 
structure in the trial-and-error, stages of structure determination. ‘The purpose 
of the present paper is to describe the construction and setting-up of a fly’s eye; 
the theory of its use is to be given in subsequent papers. 


§2. DEVELOPMENT OF THE INSTRUMENT 


The first method by which a repeated pattern was produced (Bragg, 1944) was 
by means of a pin-hole camera with an array of pin-holes lying on a square lattice. 
This was made by photographing an array of black dots on a white background, 
giving an opaque plate with an array of transparent holes. ‘These were spaced at 
0:025cm., each hole having a diameter of 0-:0012cm., and were arranged on a 
square with 26 holes along each side. ‘To produce a photograph representing a 
structure a photographic plate was supported 1 mm. above the pin-holes, and a light, 
which could be moved about in a horizontal plane 80cm. below the plate, was 
switched on for a short time in each of the several positions representing atomic 
positions in the unit cell. The pin-hole camera gave a reduction in linear dimen- 
sions of 800: 1, and the size of the “unit cell” on which the lamp was moved was, 
therefore, 0-025 x 800 =20cm. 

This method gave good results but suffered from two disadvantages: the 
pin-hole images are somewhat diffuse, and specks of dust may block the pin-holes 
and give rise to images of unequal intensities. 

It is therefore advantageous to use lenses rather than pin-holes. Small lenses 
can be made in the following way: if a steel ball-bearing is pressed into a copper 


| 
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plate, a spherical depression of nearly optical perfection is produced; if, now, a 
piece of Perspex is pressed at its softening temperature into the copper plate, a 
spherical lens will be embossed on the surface of the Perspex. 

To make a fly’s eye we require a regular array of such lenses, and we want, 
therefore, to make a regular array of depressions in a copper plate. Two ways of 
obtaining such an array are (1) by packing a large number of balls into a tray and 
pressing a copper plate on to them, (2) by making each depression individually 
and moving the copper plate through the repeat distance of the pattern between 
successive depressions. 

The first method suffers from one inherent disadvantage, namely, that it fixes 
the spacing between the depressions relative to their radius of curvature. The 
distance between adjacent balls in the tray must be twice their radius, r, and since 
the lenses will have a focal length of r/(u—1), this distance is necessarily 2(—1) 


Giza? 


Photographic 
Plate we 
p=l 


Direction of 
cell diagonal 


Figure 1. Diagram showing maximum inclination of rays to axes of lenses. 


times the focal length of the lenses, where y is the refractive index of Perspex, 
equal to 1:49, Therefore rays which would come to a focus at the corner of a 
unit cell (see figure 1) would have to be inclined to the axis at an angle of 


tan! (uw —1)V2=35°. 


This is too large an angle between rays and axis for a plano-convex lens of reason- 
able aperture, and it cannot be reduced except by making (4 —1) smaller; this 
could be done by placing another transparent medium over the lenses, but this 
would introduce complications in construction. 

The second method has the disadvantage that the spacing and diameter of the 
lenses are not so easily kept constant, but with suitable precautions it was found 
possible to make the variations in these unimportant. 


Construction * 


A cylindrical piece of copper 1-5 in. in diameter and 0-5 in. thick was polished 
on one of its plane faces, first with emery paper and then with metal polish. In 
order to provide an accurately measurable translatory motion, the copper was 
screwed firmly, with its polished side upwards, to a lathe tool-rest with two gradu- 


* All dimensions are given in inches. since machine tools and ball-bearing balls are calibrated 


in inches. 
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ated screw-motions at right angles, which in turn was screwed to the base plate of a 
drilling machine. A steel ball, 3/16in. in diameter, was fixed with wax into a || 
depression in the end of a steel rod, which was supported in the drill holder of the 
machine. A weight was suspended from the arm of the drilling machine and, | 
by lowering this weight slowly, the ball was brought down with a constant thrust 
(as in a hardness-testing machine) on the copper plate. (The usual rotary motion | 
of the drilling machine was not, of course, used.) The weight was adjusted by | 
trial so that a depression of 0-025in. diameter was made. The plate was moved | 
0-05 in. after each depression was made. A slight “burr” was thrown up around 
the depressions, and it was therefore important not to make them too large in | 
relation to their distance apart. | 

The depressions were made to lie on a square lattice with a spacing of 0-05in. | 
and to cover a circular area 1-05in. in diameter; there were 341 depressions 
altogether. Their diameters were constant to about 10%, which was sufficiently 
good for this purpose. 

The lenses were then made from this negative by pressing a piece of Perspex 
between the copper plate and a polished brass plate in a moulding press. The | 
Perspex was originally 0-15in. thick and was cut to a circle 1-5in. in diameter. 
The moulding press was heated to about 140° c., at which temperature the Perspex 
softened, and was screwed down slightly; about 0-02 in. was sufficient to emboss 
the lenses. The pressure was not released until the temperature had fallen 
below 80°c., so that the Perspex was hard again. 

The lenses so produced had a numerical aperture of about f/7:5, and were 
found by photographic tests to give good definition over the whole field of view 
employed. ; 

In order that light should pass only through the lenses, it was necessary to make 
the plane parts of the Perspex opaque. Methods of inking or painting over the 
plane parts were tried but were found impracticable, as no ink or paint could be 
found which was opaque in layers as thin as the embossed lenses. Instead, a 
mask was made by punching holes of the same diameter as the lenses in a piece of 
black paper. A needle with its end ground flat was used as a punch and the paper 
was supported on a flat piece of lead. ‘The drilling machine and lathe tool-rest 
were used, as described above, to place the holes correctly. The mask was stuck 
to the Perspex with “ Durofix”, which was applied only round the edges of the 
fly’s eye, to avoid contaminating any of the lenses. 


Setting up 


In order that a photographic plate may be supported in the focal plane of the 
lenses, three holes were drilled and tapped near the circumference of the Perspex 
disc, and three 8-B.A. screws, pointed at their tips, were inserted in these holes 
and provided with nuts for locking. The fly’s eye was fixed by these screws on toa 
circular brass plate with a hole in its centre large enough to include all the lenses. 
This was mounted on a horizontal shelf of adjustable height with the lenses 
downwards and the tips of the screws upwards (figure 2). Below this was a base 
plate on which a lamp could be moved to various known positions by reference to 
two scales at right angles. ‘This could be replaced when desired (see below) by an 
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illuminated ground-glass screen. The whole apparatus must, of course, be used in 
a darkened room. 


Photographic Plate 
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Figure 2. Set-up of fly’s eye. 


Adjustments 


(1) Focus. ‘The three screws were adjusted, before the fly’s eye was finally 
set up, by placing the fly’s eve on a microscope stage with a piece of ground glass 
resting on the screws, and turning the screws so that all the images of a distant 
object were seen in the same plane as the ground glass. A blue filter was used so as 
to obtain the focus for blue light. 

A slight adjustment of the focus was now made to allow for the fact that the 
photographic plate responded mainly to ultra-violet light. This was done by 
lowering one of the screws by about half a turn and taking a fly’s eye photograph 
of the lamp in a number of positions. Since the photographic plate is tilted during 


(0,20) (20,20) 
(20,12) 
0,10 
10) (20,8) 
(0,0) (20,0) 


Figure 3. Unit cell and positions of lamp. 


this exposure, the images are not in focus all over the plate; by finding those images 
for which the definition is best, the height of the focal plane could be calculated, 
and hence the screws were adjusted to the correct height. 

(ii) Scale and orientation. ‘Two adjustments are now necessary to ensure that 
the contents of the unit cell shall be reproduced in the fly’s-eye photograph at the 
correct scale and in the correct orientation. ‘This was done as follows: a photo- 
graph was taken for which the lamp was switched on in positions whose coordinates 
were (0,10), (20,8), (20,12), measured in cm., the “unit cell” being 20 x 20 
cm. (figure 3), The fly’s eye should, if correctly adjusted, give a pattern as 
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shown in figure 4 (a), where the three black circles represent the spots belonging to | 
one unit cell. A typical out-of-adjustment photograph is represented by figure 4(b), 
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Figure 4. Photographs obtained with (a) correct adjustment, (6) incorrect adjustment. 


which shows that a rotation of the fly’s eye and a reduction of scale are necessary. 
To reduce the scale, the shelf on which the fly’s eye is mounted must be raised. 


§3. METHOD OF PHOTOGRAPHING A PROPOSED STRUCTURE 


The method used with the original pin-hole fly’s eye was to switch on the lamp 
for a suitable length of time in each of the positions at which an atom was to be 
represented, so that each atom was reproduced as a black spot in each unit cell of the 
fly’s-eye photograph. A negative of this photograph, in which each atom would be 
represented by a white spot, was then taken, and diffraction spectra from it were 
observed. ; 

. This method may also be used for the lens fly’s eye, but a more convenient 
method is available. Ifthe lamp is replaced by a uniformly illuminated screen, 
20 cm. square, and a number of opaque discs are placed on the screen, the fly’s-eye 
photograph will have a black background with light circles on it representing 
atoms. This method has the advantages that no negative need be taken, and that 
only one exposure is needed, once the discs have been arranged on the screen, 
instead of one for each atom. 

In practice, a difficulty arises here because the lenses of the fly’s eye may not all 
be accurately spaced, and, therefore, the black squares which are the images of the 
light square 20 x 20 cm. may overlap in some places and fail to meet in others. 
Overlapping does not matter provided the plate is sufficiently over-exposed for all 
the black part to be practically opaque: the failure of some of the squares to meet 
may be overcome by making the illuminated squares slightly larger than 20 x 
20 cm., so that all the squares overlap a little. Regions near the edge of the unit cell 
in the fly’s-eye image of the screen are now illuminated twice over, by light from 
opposite edges of the screen, and no parts are left unilluminated. Any atom that 
lies on or near the edge of the unit cell must now be represented twice over, as 
shown at AA in figure 5 (a), to ensure that the image of A shall bea complete shadow. 
Similarly an atom near the corner of the cell must be represented four times over, 
as shown in figure 5 (4). 
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Figure 5. Repetition of atom near edge or corner of unit cell. 


S44), THE DIFFRACTION APPARATUS 


The apparatus for observing the diffraction spectra was similar to that already 
described by Bragg (1944). It consisted essentially of a mercury-vapour arc lamp 
with a pin-hole stop, a lens of about a metre focal length to produce a collimated 
beam, in which the two-dimensional grating was placed, and a similar lens to focus 
the spectra. Brass tubes were arranged to eliminate extraneous light. The 
separation of the spectra thus produced was only about 0-3 mm.; an eyepiece or 
low-power microscope was therefore used for observing the spectra. A green or 
deep blue filter was found useful for separating a single line of the mercury 
spectrum; this was almost essential when the spectra were to be photographed. 

As the emulsion and the glass of the grating are not of very uniform thickness, 
the diffraction pattern may be improved considerably by placing the plate between 
two optical flats and filling the spaces with cedarwood oil. 


SBS Db OF HY Ssh VE ON DURE NE Ss PROC TURE 


The fly’s eye was tried out on the known structure of durene (Robertson, 1933). 
The unit cell has not a square cross-section, but if all the atomic-position coordin- 
ates are referred to equal axes at right angles the intensities of diffraction are not 
altered. The projection of the structure along the b-axis was taken, and since this 
gives two identical parts to each unit cell, the cell was halved before referring it to 
equal axes at right angles. Each carbon atom was represented by a disc laid on an 
illuminated screen as described above, while the hydrogen atoms were neglected. 
An enlarged reproduction of the fly’s-eye photograph is shown in figure 6, while in 
figures 7(a) and 7(b) are shown the diffraction photographs obtained from this: 
the first was taken with the fly’s-eye photograph enclosed between optical flats 
and cedarwood oil, while the second was taken without the use of optical flats. 
The figures under the spots in figure 7 (a) are relative intensities calculated from 
Robertson’s data. The agreement is seen to be good when allowance is made for the 
‘act that the discs do not properly represent the carbon atoms, giving a diffracted 
ntensity which does not fall off rapidly enough with angle; also the effect of 


emperature on intensity is neglected. 
21-2 
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Figure 6. Fly’s-eye photograph of durene. 


(a) (0) 


Figure 7. Diffraction pattern from durene grating: (a) optical flats and cedarwood oil used ; 
(b) no optical flats used. 
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DISCUSS LON 


Mr. D. E, THomas. Attempts to follow Dr. Stokes and Prof. Bragg’s technique have 
been made in my laboratory at the Research Department, Woolwich, and one or two 
minor differences have arisen. The original mould was made with a Vickers hardness- 
testing machine, the brass block being held in a clamp which could be traversed in two 
directions at right angles by means of micrometer screws. The technique of making the 
mask differed considerably. An image of the Perspex fly’s eye was made by focusing a 
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distant small source of light on to a lantern slide. The image was intensified by multiple 
printing until it gave clear holes on a black background. It was then fitted immediately 
over the fly’s eye. 

After forming the diffraction grating it was reduced to one-third (linear dimensions) 
ona fine resolution plate, and this small grating used to give the diffraction patterns of a 
point source of monochromatic light. These patterns can be seen with the naked eye and 
photographed, and are of a convenient size for records. . 


Dr. R. S. Cray. I suggest that the author could have made the impressions more 
easily if he had used the group of balls in the frame, and obtained closer grouping of the 
impressions by moving the frame, after making one impression, an exact fraction of the 
distance apart of the balls (e.g. a third of the distance apart), then making a second 


impression, and so on. With nine impressions (in the case assumed) he would get a 
complete mould. 


Prof. A. O. Rankine. Dr. Stokes’s paper emphasizes a certain ambiguity in the use 
of the word “spectra”? as between x-ray and ordinary light phenomena. It sounds a 
little odd to hear of a blue filter being inserted to render the spectra more definite, although 
the meaning, on contemplation, is evident enough. Must diffraction patterns originating 
from single wave-length radiation on the one hand, and arrays of images produced by 
composite wave-length radiation on the other, both be designated by the same term— 
Spectra ? 


AuTHoR’s reply. I am interested in Mr. Thomas’s method of making and using his 
fly’s eye. With regard to Dr. Clay’s suggestion, I doubt whether the advantage of this 
method would outweigh the trouble of both fitting balls in a frame and providing the 
n-cessary translatory motion ; moreover, a much greater force would be needed in 
order to produce many indentations at once. It is, incidentally, somewhat difficult to 
persuade balls to form a square array in a tray ; they always tend to go into a hexagonal 
pattern. 

In answer to Prof. Rankine’s point about spectra, my view is that, in x-ray or optical 
diffraction, one spectrum appears on a photographic plate generally as a series of spots, 
one for each wave-length of the incident waves ; if these be monochromatic, then a 
spectrum consists of one spot only. The several spots which arise from the different 
orders of diffraction may thus be correctly called spectra. It is, however, more precise, 
though less convenient, to speak of “ orders of diffraction’ rather than “ spectra’. 
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FLUCTUATION NOISE IN A RECEIVING AERIAL 
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ABSTRACT. The factors which determine the signal/noise ratio at the terminals of 
a receiving aerial are discussed. ‘The aerial noise considered is the random fluctuation 
type, consisting of 


(i) thermal noise associated with the loss resistance of the aerial, and 
(ii) noise associated with the radiation resistance which is induced by the surroundings. 


‘The effective temperature of the radiation resistance is expressed in terms of the tem- 
perature distribution of the surroundings and the distribution of power dissipation when 
the aerial is transmitting. Radiation from the sun and from the Milky Way are briefly 
discussed, and it is shown that the detection of solar radiation at radio frequencies requires 
the use of highly directional aerials. 
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The limitations imposed by the receiver noise on 
(a) the sensitivity of the reception of signals, and | 
(b) the accuracy of measurement of aerial noise, are discussed and the results presented | 

graphically. | 

The conclusions are of most practical interest at the higher radio frequencies (above 


about 20 Mc./s.) where atmospheric noise is negligible. | 


§1. INTRODUCTION 


. 
an enclosure at temperature 7, the mean square noise e.m.f. has the value | 


I: has been shown (Burgess, 1941) that for an aerial at thermal equilibrium in i 
given by Nyquist’s theorem, 


€,2=4RTRB =4kT(Ry+R,)B, veers (1) 
where R is the total resistance and B is the band-width of the receiving system. 
This noise is composed partly of the intrinsic thermal agitation e.m.f. generated in 
the loss resistance R, of the aerial and partly of the noise received from the sur- 
roundings which can be attributed to the radiation resistance R,. 

In practical aerial systems these ideal conditions do not obtain, and the effective | 
temperatures 7, and 7", of the ohmic and radiation resistances will not in general be 
equal. The mean square noise e.m.f. is thus 

€,° =4kR(T Ry T,R,)B. PB 0) | 

It is assumed throughout that all the noise is of the random fluctuation type, | 
in which the mean square e.m.f. is strictly proportional to band-width. 

The present note is supplementary to the earlier paper, and discusses the | 
influence of the aerial characteristics and the receiver noise on the limiting sensi- | 


tivity of signal reception and on the detection and accuracy of measurement of 
noise fields. 


§2. THE EQUIVALENT CIRCUIT OF A RECEIVING AERIAL-FOR 
SIGNAL AND NOISE 
The effective height of an aerial for reception is a function of the direction of 
arrival and polarization of the received wave. If the wave has electric intensity E 


at some reference point in the aerial system, the induced signal e.m.f. appearing at 
the terminals is given by 


€,= ERO; d:c).- - | oe a err (3) 
where A is a function of the direction of arrival (6, 4) and polarization («) as well as of 
frequency and of the choice of terminals. ; 


The effective height for transmission has the same value and the radiation 
resistance of the aerial can be found readily from Poynting’s theorem: 


SOs 
R= da. [42 dQ, Se GS 
Jo 
where A is the wave-length, while dQ =sin 6. d@ .d¢ is the element of solid angle. 
The equation gives R, in ohms if X and A are in metres; M.K.S. units are 
assumed throughout this paper. Itisnot customary to take account of the depen- 
dence of effective height on polarization, but this is essential in the present analysis 
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since, in the transmitting condition, an aerial in general produces elliptical polariza- 
tion ata distant point. Furthermore, the random polarization of noise fields requires 
integration over all angles of polarization in order to calculate the induced noise 
e.m.f. in the case of reception. 

The directivity or gain of an aerial is a measure of the extent to which its polar 
diagram is concentrated in one direction and for one polarization. Thus if h,, is 
the maximum effective height which occurs for this direction and polarization, the 
directivity can be defined by 

4n*h 2 
= 2 ee a 
| "dx [42 do 
0 fw 
The directivity is unity for a hypothetical omnidirectional aerial giving linear 
polarization (h=h,,sin«). For linearly polarized aerials the above equation 
agrees with the usually accepted definition of directivity, e.g. for a Hertzian dipole 
D=1-5 and for a half-wave dipole D= 1-64. 
Combining equations (4) and (5) gives 


_ 1207*h 

MD 
The maximum available power for such an aerial in a plane wave of electric intensity 
E volts/metre is 


ie. (ohms);? <0 FRYE: = ots (6) 


Eh?  E22D 


= aR, 480n 
St A Didarwatts, 9 aes (7) 
FE? 
where Ee =T00 watts/m? 


is the energy flux in the wave. 


In practice the resistance R of an aerial is the sum of the radiation resistance R, 
and the resistance R,, which represents the losses in the conductors, the dielectrics 
and the earth (and in the terminating resistance in the case of a progressive wave 
system), all at the ambient temperature 7). It is convenient to write for the ratio 
of the radiation resistance to the total resistance 

R, 
Vani? 
wherey may for obvious reasons betermed the efficiency. For aresonant dipole y 
may closely approach unity, but for a short aerial or a progressive wave-system it 
may be small. 

Now the effective noise temperature of a linear two-terminal system whose 
elements are at different temperatures 7’, is readily shown by an application of the 
Reciprocity Theorem to be given by 


IO, ONE 


SE, a Otte (8) 


where P,, is the power dissipated in the kth element when a generator is applied to 
the terminals of the system and P is the total power dissipated. 


/ 
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This equation is directly applicable to an aerial. The loss resistance Ro, 
which corresponds to the power dissipated in the aerial itself and in its immediate 
surroundings, is at the ambient temperature 7). The radiation resistance k, 
corresponds to the power which will eventually be dissipated at a relatively 
great distance from the aerial, and thus its effective temperature cannot be so readily 
determined. Forexample, onlow radio frequencies (say below 1 Mc./s.), radiation 
can neither escape from nor enter the space bounded by the earth and the iono- 
sphere, and in this case the effective temperature 7’, of the radiation resistance will 
lie somewhere between the temperature of these boundaries. At higher radio 
frequencies, where penetration of the ionosphere occurs, the temperature of bodies 
lying outside it must be taken into account. 

The equivalent circuit of the receiving aerial for signal and noise may now be 
readily constructed. The signal e.m.f. is hE, where h is the effective height for the 
particular conditions of the incident wave, and will be 4,,, if adjustment is made for 
maximum reception. The noise e.m.f.s are V(4R7T)R,B) in series with the 


ohmic resistance Ry and V(4Rk7..R,B) in series with the radiation resistance R,. 
The noise e.m.f.s may be combined by ascribing an effective temperature 7” to 
the whole aerial resistance R, where 

ly Teena 
Se aa 
The overall signal/noise ratio at the aerial terminals in the condition of maximum 
reception is thus 


see iL EVE) Tae ae (9) 


Cy hl 
a za 
e. a Cee T,) 
. _ AE eee (10) 
V 48072RBT, Ei 1. lero 
+ —= 
ye dT, 


The first factor represents the signal/noise ratio of a perfectly efficient aerial (y= 1) 
of directivity unity. To render practical interpretation easy, figure 1 shows the 
field strength for which the signal/noise ratio is unity in a loss-free Hertzian dipole 
(D=1-5, y=1) with a band-width B of 10kc./s. and an ambient temperature of 
290°k. ‘The corresponding equation is 


V(320m?kBT,) 0:35 |T 
1 — ns Tr. eeeees (i) 


(uv./m. when the wave-length is in metres). 

It is seen from equation (10) that if the noise field is relatively large (T,,>> Ty) the 
efficiency of the aerial is of secondary importance. It is on account of the pre- 
valence of high noise levels that aerial systems, which theoretically appear to be 
poor noise collectors of signal energy, are in practice found to be little inferior to 
systems designed to be of high efficiency. 

If the noise radiation is not perfectly diffuse, 7, will itself be a function of the 
directivity and, in fact, may either increase or decrease as D increases, depending 
upon whether the polar diagram is concentrated on a region of relatively higher or 
lower temperature. . 
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§3. THE EFFECTIVE TEMPERATURE OF THE RADIATION 
RESISTANCE 


The density of black-body radiant energy at absolute temperature T in a 
frequency band B is given by the Rayleigh-Jeans law at radio frequencies : 


87kRT 
CSTE 18}. 
The energy flux received from a solid angle dQ is given by 
dieu. 2kTB 
df= Wore oe dl ere dQ. watt/m?, 


where dE,” is the element of mean square noise field intensity.* Thus the total 
mean square noise e.m.f. induced in an aerial from its surroundings is given by 


2407RB (?2"™dax 
tae | LS el etal 
Cea he - | Th dQ. 


‘Therefore the effective temperature of the radiation resistance is given by 


| ae i Th .dQ 
pee See ee Pa 1 ed or) Sa 12 
iN DEAS [ "dx. [22.dO a 


by virtue of equations (5) and (6). This formula for 7, is seen to be formally 
identical with equation (8) for a linear network, and, in fact, h?. dQ. dx is propor- 
tional to the element of power dissipated in a direction dQ for angles of polarization 
between « and (« + da). 

If the source of radiation is a black-body at temperature T which subtends only 
a small solid angle Q at the aerial, and furthermore, if the aerial is linearly polarized 
and directed towards the body, equation (12) assumes the simple form 


QD 
Does Dae setae (13) 
Ne 34: 
provided that the aerial directivity D is small compared with = To take a 
specific instance, consider the case of black-body radiation from the sun for which 
approximately 
== O000- Ki = 7 x0, 
F=1-12x10-/\? watts/m? per Mc./s., 
I,=0-03D°x., 


and thus it will be appreciated that very high directivity of the order of 10 is 
required if the solar radiation is to contribute appreciably to the aerial noise. 
Furthermore, it is not certain that the sun will radiate as a black-body or that 
absorption in the intervening space will be negligible, and it is, therefore, possible 
that the noise field at the aerial will differ from the ideal case calculated above. 
Measurements of noise at radio frequencies of about 20 Mc./s. by Jansky (1935 

* dEy* is twice the mean square field measured in any one direction in the wave-front, the 


factor of 2 taking into account the two independent orthogonal polarizations which contribute 
equally to the energy flux of a randomly polarized wave. 
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and 1937) did not reveal any evidence of solar radiation, nor did those of Reber | 
(1940 and 1942) at higher frequencies (160, 900 and 3300 Mc./s.).* 
Jansky (1935 and 1937), Reber (1940 and 1942), and Franz (1942) have established | 
that an intense source of fluctuation noise exists in a fixed region in space in the | 
direction of the Milky Way. The exact nature and origin of this radiation has not 
yet been conclusively established, but if it is thermal, the temperatures at the 
source must be very high (probably greater than 100,000°xK.). ‘The discussion by | 
Henyey and Keenan (1940) is of interest in this connection. | 
It is suggested that refined measurements of noise radiation at relatively high 
frequencies for which the ionosphere and the earth’s atmosphere are sensibly 
transparent, and at which directive aerials are readily constructed, would provide 
much useful information regarding the temperature of, and processes in, regions of | 
outer space | 


| 


§4. THE SIGNIFICANCE OF RECEIVER NOISE 


If an aerial is attached to an ideal receiver which produces no noise, the output | 
signal/noise ratio (measured at a point prior to detection) would be the same as that | 
at the aerial. Similarly any source of noise radiation whose ne. acts on the aerial | 
could be detected provided that sufficient amplification were available. | 

In practice a receiver contributes to the output noise level, so resulting in a ' 
reduction of signal/noise ratio in the case of signal reception or in masking of the _ 
aerial noise where measurement of the latter is desired. 

A convenient measure of receiver noise is the “‘ noise factor’ usually denoted by 
the symbol N. A receiver of noise factor N connected to a source of signal, 
whose resistive component R of impedance is at temperature 79, will give a signal/ 
noise ratio which is 1/N times that given by an ideal noise-free receiver of the same 
band-width connected to the same source. The noise e.m.f. of the receiver itself 
referred to the aerial terminals is, therefore, 


é2=4kT(N—1)RB. |) 2 eee (14) 


(In practice it is more convenient to express the noise factor as 10 log,) N decibels, 
but in this equation and those that follow the numerical power ratio is used.) It is 
seen from this equation that the noise of the receiver can be expressed as an 
equivalent field intensity by taking (N —1) for T,/7, in figure 1. 
The loss in signal/noise ratio when the receiver is joined to an aerial at effective 
temperature 7” will not in general be N, but will be given by 
: 6.2 
Ni=l+gtten =14(N-1)38, RA (15) 
which is greater or less than N according as 7” is less or greater than Ty. In 
figure 2, the effective noise factor N’ is shown as a function of T’/T, for various 
values of N. It is readily seen that the values of T’/T, which occur in practice 
result in a considerable improvement in the effective noise factor, e.g. if for a 
given receiver N= 25 (14 db.), then for T’/T,=10 it is found that N’ =5-3 db., 


* Since the present paper was written, Reber (1944) and Southworth (1945) report measure- 


ments of solar noise, using highly directional aerials which give intensities of the same order as 
the black-body value calculated above. 


Fluctuation noise in a receiving aerial 220 


while for T’/T,=100, N’=0-9 db. This serves to explain why in practice an 
improved design of receiver or aerial frequently fails to give the expected improve- 
ment in sensitivity. 
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Figure 1. Equivalent noise field intensity for a dipole receiving radiation of temperature Ty. 
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Figure 2. The influence of the aerial noise on the effective noise factor of a receiver. 


The effect of receiver noise on the observation of aerial noise can readily be 
estimated. A convenient method of determining the aerial noise (i.e. 7’) is to 
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measure the ratio m of the output noise powers when the receiver is connected in 
turn to the aerial (at 7’) and to the dummy (at 7\): 


ees ee 
\4RT RB Ce ta ee 


nN 


Figure 3 shows the noise ratio m as a function of T’/T, for various values of N. 
The accuracy of the measurement of 7’ can be expressed by means of the differ- 
ential ratio: 


Ty , 1 
97 | a LA (N= la > cuore eege) 
Thus the fractional uncertainity of 7” is always greater than that in n, and if 


AT’/T” becomes the order of unity it may be said that the aerial noise is not 
measurable. 
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Figure 3. The influence of receiver noise on the measurement of aerial noise. 


As the frequency cf measurement increases, it is possible to obtain higher © 
aerial directivity for a given aperture and thus the effective aerial temperature TJ’ 
due to noise from a single source (such as the sun) is raised. This improvement 
is, however, slightly offset by the increase in the receiver noise factor at the higher 
frequencies. 
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ABSTRACT. Equations have been derived for the steady-state distributions of 
temperature or concentration across membranes in the form of plates, tubes or spherical 
shells, when the thermal conductivity, or diffusion, “‘ constants ”’ are functions respectively 
of temperature, or solute concentration, in the diffusion medium, and also where these 
“constants ’”? are functions of the positional co-ordinates. Equations are similarly given 
for the quantity of heat or solute transported through the membranes. 

Procedures based on these equations are given for measuring diffusion and thermal 
conductivity ‘“ constants’? and their dependence upon concentration and temperature 


respectively, or on positional co-ordinates. All observational data are based only on the 
a) 


\ 2 0c 2c 
steady state of flow, thus avoiding the need to deal with equations such as PP ral Pax) 


which have so far not been integrated. "The methods are suitable for measurements 
of flow of liquids, vapours or dyestuffs through natural and artificial plastics, of gases 
and liquids through zeolitic media, in surface migration of adatoms on metals, and similar 


problems. 


§1. INTRODUCTION 
HE focus of attention in quantitative studies of the flow of matter or of heat 
has long been Fick’s law : 


Oc 
i) ieee) 01 ist bee Sideabea oe 1 
a NEC (1) 
or 
Ss) anit Yon, ie | dae beeen la 
at V (1a) 


where c and 7 are concentration and temperature respectively, and D is the dif- 
fusion or thermal diffusivity constant. It is now clear that D is often not a 
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constant at all, so that this law is a limiting law only, just as in other spheres are the 
Langmuir isotherm, the Debye-Hiickel equation or Raoult’s law. Equation (1) 
must then be replaced by the following equation, which allows for the non- 


constancy of D: 
Glee Gl dc 0 dc 0 0c 
pale a = ae = =} ye Re oes 2 
a1 5 (P x)+ 5(P 5) 5 (2 a ) 


No exact solution of equation (2) has been found for any boundary conditions, 
although in several instances approximate solutions have been derived (e.g. 
Hopkins, 1938). Variants of a graphical method have been evolved (Matano, 1933; 
Eversole, Peterson and Kindsvater, 1938) for determining D as a function of c from 


the equation 
ae 0c 
ae 5 (P a) apes (2a) 


The method has been used by Wells and Mehl (1941), Matano (1933), Mehl (1937) 
and Rhines and Mehl (1937) for diffusion in metallic media and for diffusion of 
salts in solutions by Eversole et al., 1941, 1942). 

In general D may be a function of concentration of the diffusing species, or a 
function of the positional co-ordinates, x, y and x. D depends on c when metals 
and non-metals diffuse into metals (see above references; also Wells and Mehl, 
1940 ; Smoluchowski, 1942, 1943 ; Barrer, 1941a); salts diffuse in solution 
(Eversole et al., 1941, 1942); gases and liquids diffuse in zeolitic media (Barrer, 
1941b; Barrer and Ibbitson, 1944); adatoms diffuse on metallic surfaces (Lang- 
muir, 1932, 1934; Langmuir and Kingdon, 1925; Bosworth, 1935, 1937); and 
liquids or vapours are imbibed by polymers which thereby swell (Daynes, 1937; 
Barrer, 1941c; Garvie and Neal, 1938; King, 1945). D depends on positional 
co-ordinates when heat or matter flows in physically or chemically inhomogeneous 
media. ‘he unequal thermal gradients which arise at different places on the 
earth’s crust (Bullard, 1945), or the permeability of inhomogeneous rock strata to 
fluids such as petroleum, afford examples of interest. Diffusion of thorium along 
a tungsten filament where a temperature gradient exists along the filament is a 
more complex example, for here D depends both on the thorium concentration 
(Langmuir, 1932, 1934) and also on the positional co-ordinate, x. Dependence 
on x in such problems arises independently of dependence on ¢ because D isa 


E 
function of T ( D=Dye°2f and T depends on x. 

The rather laborious character of the graphical method of measuring D as a 
function of ¢, and the absence of a procedure when D = f(x, y, z),render it important 
to evolve new simple procedures. ‘The theory of diffusion in condensed phases 
often awaits methods which permit data to be readily obtained. In this paper 
procedures are developed for measuring D when this “ constant” depends either 
on concentration or on the positional co-ordinates. 'The new methods require 
observation of the stationary state of flow only, with a corresponding increase in 
directness and practicality. They apply equally to the diffusion of heat or of 
matter. 


| 


| 
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i 
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§2. THE FUNDAMENTAL EQUATIONS 


In Cartesian co-ordinates, equation (2) is the fundamental expression for 
material diffusion; its counterpart for flow of heat is 


ii hae at) oT 7) oT ry) oT 
Gi Sie CL ne es Ole 2k 
Caaret a(x =) 5(K 5, )+ 5(K ff ) nae (25) 


In these equations D has the dimensions /?¢, and K, the thermal conductivity 
constant, the dimensions m'J-1t-!._ Also C,, is the heat capacity at the point 
(x, y, 2), and p the density at this point. Equation (2) is thus a particular case of 
(2b), where C,,p is independent of the dependent or independent variables. How- 
ever, when D and K are true constants, equations (2) and (26) reduce to equations 
(1) and (1a) respectively, the thermal diffusivity constant D=K/C,p in equation 
(1 a) now having the same dimensions as the diffusion constant D in equation (1). 

Transformation of equations (2) and (26) to cylindrical co-ordinates gave for 
flow of heat and matter, respectively, 


oY AG OLN EA oT 0 oT 
lous Sorts Jaa ® 7) * i oe (x a rete (3) 
dc_ 1d CCG. dc 7) dc 
Ot oT A" Ds\+ D (PD x) =-(P =) Ke ores (3a) 
Transformation to polar co-ordinates gave similarly 


am eyes hing or ee eae 
Cp. = Sr KG, )+ faa an heer se 78) 


Oca lee dc tea, dc 1G dc 
ot 5dr De) + sind sol sin pss) + 7? sin? 0 5(P 55) 


Many important problems concern diffusion in cylinders and spheres with con- 
centration or temperature independent of @ and z (for cylinders) or of 6 and ¢ 
(for spheres), so that the corresponding equations follow from equations (3) and 
(4) on omitting terms in these variables. In the steady state in which we are 


interested s ss = = 0, and for heat flow, for example, one has 
A p ("x =) Cpe 1 nee (5) 
v= 5 (rk =) (cylinder). .§ 7 > “esease (6) 
To these may be added the equation for heat flow in the x-direction only: 
= = K =) ities Boe. A eS (7) 


The corresponding equations for diffusion of matter follow at once by 
substituting D for K and c for T. 


| 
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The equations (5), (6) and (7) may be integrated readily. Thus for equation (5) 
de aT | 


ayes PR ad) re 8) 
PDT or PK (8) 


where A is aconstant. Now D (or K) may bea function of ¢ or of 7(x). In the) 
former case one has quite generally 

D=Dil+foly eee (9) >} 
and in the latter 

D=D 5 [LtfO” yo ee ae (10) || 
In equation (9), f(c) is conventionally chosen so that Dy is the value of D when | 
c->0; in equation (10), similarly, Dy is the value of D at a chosen facer=7,. In| 
either case the variables may be nacre and one obtains for the steady-state | 
distribution of solute across a spherical shell: 


A 
D,{[1+folde= -—< +B cae i 
when D is a function of c, and 
peel sole 
= lance * eee | 
when Disafunctionofr. Similarly for the distribution across a hollow cylindrical | 
tube 


Dy{(1+f(c))de= A ing bo Se ee (13) 
= dr 
| Dil ee Tale. oleae (14) 
while for a plate 
D of (1 +f(\\de=Axt+ Bo ~~ ores (15) 
Dee 4 ene Bi be 8 ieee (16) 


§3. DETERMINATION OF D, AND f(c) WHEN D=D,(1+f(0)] 
Equations (11), (13) and (15) all provide means of measuring both D, and f(c) 
when D=D,[1+/(c)]. 
When the boundary conditions for diffusion across the plate are. 
c=c, for all Hat v0; 
e=c, for alll? at x=, 
equation (15) for the stationary state becomes 
¢, + F(c,)—(¢+ F(c)) _ x (17) : 
Cr+ P(e) —(eg+ Flea) gee 
where* F(c) = I f(c)de. Similarly for the tube or spherical shell with boundary 
conditions 
c=¢, for all ¢ at r=n,, 
C=C, for all ¢ at r=r,, 


* Integration is between the limits 0 and c. F(c) thus tends to zero as c>0, the laws of flow 
in the limit then reducing to the simple Fick laws in accord with experiment. 
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one obtains from equation (13) (the tube) 


y+ F(c,)-(c+ F(c)) — Inr,—Inr 


€,+ F(c,)—(e,+ Fles)) nr, —In7, ee) 
and from equation (11) (the spherical shell) 
e+ P(e) — (c+ F(c)) £ ica (19) 


SiG SCPeric)) n+, Pets 
Thus the steady state concentration distributions depend on the way in ee D 


depends on concentration. From this distribution F(c) and so f(c)= <( F(c)) 


may be derived. Alternatively, one may derive © , which, when c, = 0, is (from 
x 


de y+ F(q) 
dx l(i 1(1 +f(c)) 


function of c, one sees the intercept on the axis of x to be 


Then on plotting — sae as a 


equation (17) for example as 
C+ F(e) 
l 


, and from 


values of a between c¢, and 0, f(c) easily follows. 
; x 


Dy in equation (9) can now be obtained by measuring the steady-state flow of 
solute through the medium. This flow per unit area and per unit time is 


P= -(Dg) . eee (20) 


By differentiation of equation (17) with respect to x one finds ae and then 


by integration of (p =) from 0 to ¢ the quantity O of solute passing in time t¢ 
X} c=l 


through a membrane of area A is found: 


t 
Q=4| Pdt= [+ Fe)—(a+ Fey) |e See (21) 
0 
For a tube of length Z a similar procedure gives 
= y+ F(G) — (C+ =e) t 22 
Q=2nrel| ‘Pat =2n LDg) eeeoie aera (22) 


and for a spherical shell 


aft r. 
O=4nr2 | Pat= SPEC) = (CaF (Gy) ty nde (23) 
0 
Since F(c) has previously been determined using equations (17), (18), or (19), 
D, can be evaluated in the appropriate equation (21), (22), or (23). 
An alternative procedure for determining Dy and f(c) makes use only of 
equations (21), (22) or (23). Whence, <c,, equation (21), for instance, becomes 
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where all quantities on the left-hand side of equation (24 a) are measured. If, now, | 
c, is made small, departures from the simple Fick laws tend to zero, so that 


(inh oes Ep: Doe (215) 


1>0 tAc, 


Next F(c) is given by the relation 
l 
221) oh SRG) ae Speen eee (24) 
0 
the left-hand side of which is determined for a range of values of c. Finally | 


fle) = (FO). 


§4. DETERMINATION OF D, AND f(x) WHEN D=D,(1+f(#)) 


We now return to consideration of equations (12), (14) and (16) for the steady- 
state distribution of cwhen D=D,[1+/(x)]. Ifthe integral on the right-hand side | 
of each of these equations be denoted by /, and by J, and J, at each boundary 
surface, the constants A and B may be eliminated to give for the boundary con- 
ditions in the previous section 


Moreover, the quantity Q of solute diffusing through the outgoing surface of a 
plate of area A is 


t dc CLG 
=A| er pe ota tml Pein Ai 
O : ( a) a ADot ia (26) 


while for a hollow tube of length L 


t dc ce 
= D, L| _ D — = : 2 hd sa eat 
O=2nr, : ( 5) aed 2rL. eae (27) 


and for the spherical shell 


dc 
dr 


C1 — Ce 


Le 


t 
Q=4nr,*| ~ (> ) dt = —4nDpt. 
0 r=1o 


Then to determine D, and /(x) for the plate, for example, the steady-state con- 
centration distribution can be measured, and from it the value of 


dc _ 4% 1 
dx 1,—I,(1+f(x)) 


| 
Ceaaic . 
1? and from its values for 
Jee 


Cee dc : ; ; 
n plotting Ja B8Ainst x, its value at x=0 is 


0<x <1, I(x) may be found. Equation (26) then serves to determine Dy. 
. Details given in this and the previous section apply equally well to the deter- 
mination of Ky and f(T) or f(x) in the analogous problem of heat flow. 
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§5. CHARACTER OF STEADY-STATE CONCENTRATION DISTRIBUTION 
WHEN D DEPENDS ON c OR ON x 


In equations (9) or (10) the functions f(c) and f(x) (or f(r)) must be such that D 
is never negative or infinite. For a suitable choice of constants appropriate 
functions include 


f(e)= + ac" + be™ + dc? ....(exponents positive or zero), ...... (29) 
ac 

= nes stOU 

WOeS Saas (30) 

TC) =F gees, pee 2031) 

fc) = + acre#¥e (n>0). Eos (32) 


CC —_——_ 


i, 
Ld 


be 


D x 10° (cm* sec 


Atom % Sn, Be, Si or Zn 


Figure 1 a. Experimental (Rhines and Mehl, 1938) (full curve) and calculated (dotted curve) relations 
between D and c. Calculated curves employ the exponential functions D=D,(1+ae*) and 
D=D,(1+ac*e-4/e), 

Calculated Curve (1): Sn—Cu. D=7:7 Xx 10719(1+-4:46 x 10-2 €0°862c) cm?sec>}, 

Curve (2): Si—Cu. D=2-31 x 1071°(1+0-306 €0°402c) cm?sec=1. 

Curve (3): Be—Cu. D=4-23 x 10-1%(1+7:08 x 10-2 €0°352c) cm?secr}. 

2°304 
Curve (4): Zn—Cu. D=3-46x 10-10 14+1:51x 10-33 = cm? sec:1. 


>} 


29 
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Figures 1 a and 16 give actual and calculated D-c curves according to several 
of these equations, and illustrate the degree of success obtained with them. 

The steady-state concentration distribution* across a plate of the diffusion 
medium is given by equation (20), which for c, ~0 becomes 


c+ Me)= (1-7) let Fa) Pee 7a) 


* In the absence of rate-controlling phase-boundary processes (For a discussion of these 
processes see Barrer, 1941 a). 
22-2 
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It is convenient to choose units for c such that c,=1; then in equations (29) to 
(32), 0<c<1. Also, if the solute has a limited solubility in the diffusion medium, 
c in these equations may conveniently be replaced by c/c, where c, denotes the 
saturation concentration. Using these conventions, figures 2. a and 26 give typical 
concentration distributions for some of the functions of equations (29) to (32). 
When D increases as c increases, equilibrium distributions are convex away from 
the axis of x; but when D decreases with increasing c, the distributions are concave 
towards that axis. When D at first increases and then decreases, or vice versa, 


10° x D for methylene blue (cm? sec>*) 
10° x D for copper sulphate (cm? sec>*) 


Concentration (millemoles/litre for CuSO, ; mgrm./litre for methylene blue.) 
Figure 1b. Experimental (Eversole et al., 1941, 1942), (full curve) and calculated (dotted curve) 
relations between D and c. Calculated curves employ the functions D=D,(1+ac”+bce™+- ....). 
Calculated curve 1: Methylene blue—H,O. 
D=1:-635 x 10-81 +0-0824c!/4—0-00128c +0-0000738c3!2) cm?secz?. 


Calculated curve 2: Copper sulphate—H,O. 
D=11-65 x 10-°(1 —0-00311c!/4+0-002645c1/3) cm? sec=}. 


with increasing c, the behaviour is less predictable (e.g. figure 26, curve 3). 
Curves 4 and 6 in figure 2a give the extreme range of steady-state distributions 
of concentration for a pure zeolitic diffusion process where (cf. equation (29)) 
D=D,(—eje,): 

Equilibrium distributions when D=D,[1+/(x)] follow from equation (25). 
Two typical values for f(«) will be considered : 


f(x) = ax, Aen) 
Ke) = bie GR?) a Re ere (34) 
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When f(x) =ax, I in equation (25) has the values 
1 
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Figure 2 a. 
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Typical steady-state concentration distributions across a membrane 


when D=D,(1-+f(c)). 


i(@)—ae7a—100: 

ie) —acena—105 

ii) ach a2 

f(c)=0 (simple Fick law obeyed). 
f(c)=—ac ; a=0°5. 

f(c)=—ac ; a=1°0. 


Figure 2b. Typical steady-state concentration distributions across a membrane 


when D=D,(1+/(c)). 


Curve (1): f(c)=aebe ; a=1, b=3. 
ac 
Curve (2): f(c)= jecn I= NOY), == 
Curve (3): f(c)=—act?+bc? ; a=1, b=2. 
ac 
Curve (4): f(c)= Tbe? O—Nao ae 
Curve (5): f(c)=0 (simple Fick law obeyed). 
5) - a ac . =()- = 
Curve (6): f(c)= Tbs a—0:95b—1), 
Curve (7): f(c)\=—acl? ; a=1. 
Curve (8): f(c)=—acl4 ; a=1. 


and values for J, and J, follow on setting x=0 and x=1 for the plate, or r=7, 
and r, for the tube and spherical shell. When f(x)=bx + ax*, several values of 


r={ dx 
—J14+bxe4+ ax? 
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result for different values of a and b: 


zZ 2ax+b 
———- t = on eS 4 = )7) 
Wy ream omy remy 
1 2ax+b—Vb?—4a 


[= =. |» —_—— (b?> 4a), (for plate) ...... (38) 
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Figure 3. Typical steady-state concentration distributions when D=D,(1+/(«)). The straight 
line corresponds to thes simple Fick law (f(x)=0). 
Curve (1): f(x)=—ax ; a=0-99. 
Curve (2): f(x)=—ax ; a=0-90. 
Curve (3): f(#)=ax ; a=1-0. 
Curve (4): f(x)=ax ; a=2:0. 
Curve (5): f(x)=bx+ax? ; a=1, b=2. 
Curvel(6)i3fi) bx -- ax): .a—2°250—se 
Curve) io) —cxnna—o- 
Curve (8): f(x)=ax; a=99, 


and J, and J, follow as before. Similarly, when f(r) =6r + ar’?, the value of 


I=f dr 
~ Jr(1+br+ar®) 
for the tube is 


Ta Te b 
Splat Pepe: — 5! plates ASN ER) 


where T pate denotes the appropriate integral for the plate but with 7 replacing x. 
Finally, for the spherical shell, 


1={ dr 
~ JP(1+br + ar?) 
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: 6b, 1+br+ar? 1 6%-2a 
1S I=, In a TE = e a ae Tinea doratelexs (40) 


Figure 3 shows some steady state concentration distributions across a plate of 
unit thickness when c;=1 and c,=0; and D=D,(1+<ax) (curves 8, and 1 to 5) 
and D = D,(1 + bx + ax*) (curves 6 and 7). When D is an increasing function of x 
the curves are concave towards the axis of x, and when D is a decreasing function 
of x the ¢ : x curves are convex away from this axis.* 

There is little material available for quantitative interpretation according to 
methods developed in this paper. In rubber, steady-state concentration distri- 
butions of water concave towards the x-axis (Taylor, Herman and Kemp, 1936), 
suggest that D decreases as the moisture content grows.t There seems little 
regularity in available data, however, for in other polymers, such as viscose silk and 
wool, D for certain diffusing species increases with c (Garvie and Neal, 1938; 
King, 1945). 
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OBITUARY NOTICES 


JOHN EDMUND CALTHROP 


By the sudden death in June 1945 of J. E. Calthrop, Senior Lecturer in Physics, Queen 
Mary College has lost one of the most loyal members of its staff, and his fellow-teachers 
have lost a singularly gentle and unexacting friend. 

Calthrop was born in Lincoln in July 1891. He took his London B.Sc. degree in 1912 
as a student at the Borough Road Training College, and after two years spent in teaching 
science at the Lincoln Technical School he enlisted in the Honourable Artillery Company. 
He was severely wounded at Beaumont Hamel in November 1916, while serving as an 
infantryman ; after he had recovered from his wounds he was given a commission in the 
Royal Regiment of Artillery, and he saw active service again in Flanders with a Heavy 
Battery until the armistice. Calthrop’s war service was of a pattern with his whole life : 
he was the last man one would expect to see sitting in a job at the Base, or collecting 
decorations on lines of communication. 

After the war he entered at Pembroke College, Cambridge, taking a Natural Science 
Tripos, and in 1921 he joined the Physics staff of Queen Mary (then East London) College, 
under Professor C. H. Lees. For a period after his appointment he did research work 
in general physics and heat, publishing a number of papers in the Philosophical Magazine 
and in our Proceedings, but he soon decided that teaching was his real vocation, and for 
the last twenty years of his life he devoted himself almost entirely to the training of students 
and to doing all he could for their general welfare. He had a real affection for the young, 
and he never spared himself in their service. His interest in teaching is shown in a number 
of notes he published in the School Science Review, and in a highly original course of 
Advanced Experiments in Practical Physics which appeared in 1938 and was very favourably 
reviewed both here and abroad. His last publication, at the end of 1943, was a German 
Physics Reader. 

His chief handicap as a teacher was an unconquerable diffidence of manner. He built 
up a fine library of modern physics books in English, French and German, which he studied 
assiduously (and, to mention his only vice, under-scored freely), he took pains to find the 
best way of treating a subject and he prepared admirable lectures—which he then, in class, 
delivered in a low tone and an almost apologetic manner, as though he were ashamed of 
their content. As a result, although the best students got a lot from him, the less discerning 
tended to under-value his teaching, and I fear that he was often conscious of this. His 
most effective teaching was done in the laboratory and in informal talks with the students, 
to whom he willingly sacrificed much of his spare time. 

Calthrop was one of the most transparently honest and least self-seeking of men. 
Although of a naturally retiring disposition, and a bachelor, he was—as I have known since 
his earliest contacts with my own family—entirely at his ease with very young children, 
and outstandingly happy in his relations with them. 

From 1939 to 1945, while Queen Mary College was evacuated to Cambridge, Calthrop 
resided in King’s College. He enjoyed immensely the privilege of sharing in the corporate 
life of a great Cambridge College, and his obvious integrity, his eagerness to help and his 
evident delight in his surroundings brought him many new friends. There is much 


comfort in the thought that his last years, spent in such congenial surroundings, were among 
the happiest of his life. H. R. ROBINSON. 


JOHN JOB MANLEY 


MANLEyY was born on 13 July 1863 near Wellington, Somerset, and the family migrated 
later to Devon. At the age of 25 he became Daubeny Curator of the Laboratory of 
Magdalen College, Oxford, and after the first three years he taught practical physics and 
analytical chemistry to undergraduates reading for the Honours School of Natural Science. 
His first research work at this laboratory was carried out in collaboration with V. H. Veley 
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and resulted in a series of papers on the preparation of anhydrous nitric acid, and on the 
physical properties of aqueous solutions of this acid and of sulphuric acid. This work, 
which was completed in 1905, was recognized by the award of an M.A. degree, honoris 
causa, in 1903. 

Manley was typical of the best school of experimentalists—those who took a pride in carry- 
ing out operations with their own hands. His opportunities for research were increased by 
the award of a Research Fellowship in 1917, the year in which he took his M.A. He gave 
up his teaching work in 1923 and his curatorship in 1929, in which year, at the age of 66, 
he left Oxford and equipped a laboratory in Bournemouth. He took an Oxford D.Sc. in 
1931. He continued work at Bournemouth until 1938 as the holder of a Research 
Fellowship and then privately until his death on 24 March 1946, which occurred after 
only a fortnight’s illness. 

After his partnership with Veley had ceased, Manley devised a differential densimeter 
for the study of sea-water ; a duplicate of this apparatus, mounted on gimbals, was 
successfully used by N. P. Campbell on a voyage from London to Colombo. Later, 
whilst standardizing a set of weights, Manley encountered certain difficulties which led 
him to devise his enclosed or ‘‘ protected’ balance beam (1910), and ultimately to the 
revision of Landolt’s work on the apparent change in weight during chemical reaction 
(1912). For this revision he chose the straightforward reaction between barium chloride 
and sodium sulphate, instead of the somewhat complicated ones selected by Landolt. 
He found that the change in weight did not exceed one part in 100 million. His interest 
in the subject of the precision balance and the study of analytical weights continued up to 
his last illness, and several papers were published in the Proceedings of the Physical Society, 
the last of them, on the advantages of rhodium plating, only last year (Proc. Phys Soc. 
57, 136 (1945)). 

His experimental skill, and in particular his mastery of the art of glass blowing, was 
successfully applied to the study of extremely thin wires for platinum thermometers, 
gases at low pressures, and the interaction between mercury and helium. His method 
of purifying phosphoric acid by ozonized air, worked out in 1922, was used to good effect 
in 1929 to 1934 in an extension of Baker’s work on the intensive drying of benzene. 

J. J. Manley’s chief recreations were those of a quiet, modest man, who liked the 
open air—cycling, climbing, gardening and photography. He was a devout churchman, 
with an assured faith in things spiritual. Many will remember him as an unassuming 
worker, an interesting speaker, with that self-confidence which comes from certainty of 
knowing his subject. He was twice married, and is survived by his second wife and 
by his son, C. H. Manley, the City Analyst of Leeds and Wakefield. One of his two 


grandsons is a civil engineer. 


RAOUL FREDERIC SCHMID 


R. F. ScuMip, Fellow of the Physical Society, died in Budapest on 2 October 1943. 
He was carried away by a disease of tragic suddenness at the age of 39. A man of science, 
he worked and taught all his life. As a true scientist he professed the freedom of thought 
and human rights in an age when this was dangerous and difficult, trusting a better period 
would come. This period, however, he did not live to see. With a loving heart he helped 
everybody, teaching also through the example of his own life. 

At the age of 22 he had graduated at three Universities and possessed the diplomas 
of chemical and economic engineer and the degree of Ph.D. He was the winner of several 
scholarships and, as a Rockefeller Foundation Fellow, worked in Chicago in 1932. 
Originally a chemist, he saw in molecular spectra the firm experimental ground for 
studying the fundamental problems of chemistry. Working in this field, he added some- 
thing in nearly every domain of molecular spectroscopy with an admirable fertility. He 
united wide knowledge and exceptional memory with experimental skill and invention. 

Through the inscrutable ways of providence, his collaborator, L. Ger6, was destined 
to follow him in the short space of one year and a half. He died, at the age of 35, a victim 
of the war, at a moment when he could already expect to have got to the beginning of 
the eagerly expected new period. Their names are now closely connected by the scientific 
work they have done together. 
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To name only some milestones of this work: a minute investigation of the CO | 
spectrum in order to get as complete an analysis as possible of the term system of this | 
molecule, of such importance from the point of view of organic chemistry ; pioneer |) 
work in discovery and interpretation of perturbations and predissociations ; establish- 
ment of the dissociation schemes for CO, CN, Ny, NO and diatomic hydride and deuteride 
molecules ; an interpretation of the thermochemical data of organic chemistry and the 
emphasizing of the importance of excited atomic states from the point of view of chemistry. 

Those who part from this world continue to live in their works and in the hearts of 
those who knew and loved them. J. G. VALATIN. 


REVIEWS OF BOOKS 


Mathematical Theory of Plasticity, by W. Pracer. (Brown University, Provi- 
dence, Rhode Island, 1944.) 


A considerable literature now exists dealing with various aspects of the mechanics of 
solids beyond the elastic range. Only in one or two instances, however, as far as the 
writer is aware, has any attempt been made to co-ordinate the available information in 
the form of treatises dealing with the subject as a whole, and none of these treatises appears 
to merit the description of ‘‘ mathematical theory of plasticity ’’. 

The present manuscript therefore meets a considerable need, giving as it does a 
summary of the present state of knowledge of mechanics of solids beyond the elastic range, 
together with a study of the mathematical theory of plasticity as applied to certain types 
of ideal materials. 

In the opening two chapters the analysis of stress and strain is developed in terms 
of both Lagrangian and Eulerian co-ordinates, and some properties of the stress and strain 
tensors are discussed, particularly in relation to certain invariants of the stress and strain 
systems. Expressions for “ work” are obtained in terms of the stress and strain tensors. 
These analyses are the basis of the consideration of stress-strain relations in subsequent 
chapters. , 

The third chapter considers the relations between the stress and strain tensors for an 
isotropic elastic body. For this purpose the stress and strain tensors are decomposed 
into functions of the mean stress or strain and the stress and strain deviators. The 
equations of equilibrium and compatibility are derived, and the relations between the 
invariants of stress and strain are considered. An expression for strain energy is obtained. 

The nature of elastic breakdown is discussed, with consideration of the various 
theoretical criteria, and experimental evidence supporting them. 

In Chapter IV the stress-strain relations for a visco-plastic material (defined as a 
viscous fluid having a definite yield point) are derived from consideration of the corre- 
sponding relations for a viscous fluid, by the use of the tensor of ‘‘ over stress ”’ expressing 
the extent by which the yield limit has been exceeded. The relations corresponding to 
several arbitrary conditions of yielding are stated. 

Experimental methods of determining the coefficient of viscosity of a visco-plastic 
material are described. As a limiting case of the visco-plastic material, the stress-strain 
relations of a perfectly plastic material are derived, and the experimental evidence of Taylor 
and Quinney in this connection is discussed. 

Chapter V deals with the effects of strain hardening on the stress-strain relations of a 
visco-plastic or plastic body, suitable modification being made of invariants of stress and 
strain concerned in the relations to account for this effect. Both isotropic and anisotropic 
strain hardening are considered. The latter form is assumed to be such that the three 
principal planes of strain constitute three equivalent planes of asymmetry, and that rotations 
interchanging principal axes of strain do not change the value of invariants. 

In Chapter VI, the stress-strain relations of composite material are derived from the 
relationship obtained in previous chapters for the relatively simple materials. Materials 
composed of two or more of the following elements are considered—elastic solids, viscous 
fluids, visco-plastic materials and perfectly plastic materials. These materials include 
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those of Voigt and Maxwell. The constituent materials are either placed in series, with 
Superposition of deformations and the same stress; or in parallel, with superposition of 
stresses, and the same deformation. The modification of the stress-strain relations for a 
composite visco-plastic and elastic material used by F. Odquvist for the analysis of creep 
is stated. 

Chapter VII deals with the subject of materials displaying a gradual transition from 
the elastic to the plastic state. Geometrical methods of representation of stress and strain 
deviations are first evolved. Using these methods of analysis, and making the assumption 
that materials of this type obey a general relation of the type 


ae ip{ tao) 


where gq is stress and is strain, it is shown that the equation relating the stress (Sk) 
and strain (¢,k) tensors is of the type 
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C, » and n are constants, ¢ is time and S,, a second-order invariant. Modification of the 
stress and strain relations for materials of this type undergoing isotropic or non-isotropic 
strain hardening is discussed. 

In Chapters VIII to XI the author shows the application of the theory of plasticity 
to the problems of pure shear, torsion of a cylindrical body, and plane strain with and 
without axial symmetry. In each case general formulae are first developed, and are 
subsequently applied to the cases of the elastic, visco-plastic, St. Wenant and Prandtl- 
Reuss materials, and also to materials with gradual transition from the elastic to the plastic 


%} state. Where precise mathamatical solutions of the problems presented are not found 


possible, reference is made to various experimental methods of obtaining the values of 
stress or strain functions, including the flow, electric field, and membrane analogies. 
Numerous references are given to the publications of other workers who have dealt with 
these particular problems. 

» In the concluding chapter two boundary-value problems involving the variation of 
the stresses and displacements in certain areas of incompressible bodies, including elastic 
material, Prandtl-Reuss material and material with gradual transition to the plastic state, 
are considered. Solutions are obtained by the employment of variational principles 
postulating minimum values for certain combined functions of stresses and displacements. 

Two minor criticisms have occurred to the reviewer. Regarding the content ot the 
manuscript, it seems a matter of some regret that in view of the practical importance of 
the subject of creep or time-dependent strain, the author has given only passing notice to it. 

Regarding the actual arrangement of the text, the provision of a glossary of symbols 
used would have improved the readability of the manuscript, since these symbols are 
quite numerous and constitute a considerable tax on the reader’s memory. 

It is to be hoped that in a future manuscript the author may apply the theory developed 
in this present work to further problems of plasticity, many of which, e.g. those involved 
in the engineering design of gas turbines, pipe flanges, and aerofoils, are now acquiring 
very considerable practical importance. A. E. JOHNSON. 


Mathematical Theory of Optics, by R. K. LUNEBERG. Pp. xvii-401, (Brown 
University, Providence, Rhode Island, 1944.) $4. 


These mimeographed notes are a record of a course of lectures given at Brown Uni- 
versity in the summer of 1944. The subject-matter is dealt with in a comprehensive and 
unbiased manner, and the book should prove of value to those looking for a thorough 
treatment of the theory underlying image formation in, and the design of, optical systems. 
It is likely to be of interest to those designers who make use of trigonometric ray tracing 
and are seeking an alternative method. No mention is made of such tracing, the theory 
being developed from a point of view which differs fundamentally from, say, that adopted 
in the well known book by Conrady. Considerable use is made of the characteristic 
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function method and the power and breadth of application of this method is amply 
demonstrated. "The book may be regarded as a companion to the Cambridge Tracts by || 
Synge and by Steward and the classic, Optik, by Professor Max Born. i} 

In the first chapter the wave theory (based on Maxwell’s equations) and the ray theory || 
of optics are developed simultaneously. Transparent isotropic media are considered, || 
The kernel of this theory is the equation of the wave fronts (the eiconal equation), 
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the light rays being the orthogonal trajectories of a set of wave fronts, 4=ct. Differential | 
equations satisfied by the light rays are found and Fermat’s principle is introduced. It is | 
shown that wave fronts may be constructed with the aid of light rays ; and, conversely, 
by means of Jacobi’s theorem, a method is obtained for constructing the light rays when a | 
solution to the above equation is known. Applications are made to (a) spherical waves | 
in a homogeneous medium, (5) wave fronts in a medium of discontinuous optical properties, | 
(c) periodic waves of small wave-length. 

Hamilton’s theory is presented in Chapter II. The three characteristic functions, 
the point-, mixed-, and angle-functions, are defined and their properties and physical 
significance are studied. Integral invariants and their application to ray congruences 
are discussed and an interesting geometrical interpretation is given of Poincaré’s integral 
invariant. There are a few examples of the characteristic functions, but the important | 
applications are reserved for Chapter III, in which some of the topics studied are the | 
correction of optical instruments by aspheric surfaces, axially symmetric systems, spherical | 
aberration and coma, and the problem of designing condenser systems. 

A comprehensive treatment of first-order optics is given in Chapter IV. Considerable | 
use is made of a matrix notation in the manipulation of the mapping equations. The 
problem of first-order optics is to investigate the properties of the four-parameter system, S, | 
of rays which lie in the neighbourhood of a given ray, for example the principal ray. | 
This is shown to be equivalent to a study of the linear transformations 


(ee B Ge . 
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where X;= i) P;= ie (i=0, 1) and the matrices A, B, C and D, each of order two, are 
v LY 
such that 

A*C—C*A=0, B*D-—D*B=0, D*A—B*C=1. 
For this purpose Hamilton’s method of characteristic functions is used, and the results 
are interpreted optically in regard to the formation of images by the narrow bundle, S. 


When | C|=40, it is possible, using the above matrix equation, to express Xj and Xj, in 
terms of Py and P,, in the form 
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where yp and y, are constants and 4 and C are symmetric. In the case of an axially 


symmetric system we obtain the Gaussian first-order optics. Then 
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and by means of these matrices the properties of the Gaussian system are deduced. A 


generalization is made to the case of a bundle of rays about an oblique meridional ray in a 
system of revolution. Then 
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and we have what is called an orthogonal ray system, 
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The theory of the third-order Seidel aberrations in a system of axial symmetry is given 
in Chapter V. The image functions are developed one stage beyond the first order in 
| the neighbourhood of the principal ray. The usual expressions are found for the five 
‘|, image errors, and, in the case of a medium with a continuously variable refractive index, 
the analogous integral expressions are given. Petzval’s theorem on field curvature is 
discussed and a section deals with chromatic aberrations. 

The last chapter deals with the Diffraction Theory of Optical Instruments, that 1S5 
the propagation of light waves in a non-homogeneous optical medium. Here we return 
to the electromagnetic theory. The diffraction of unpolarized light and the diffraction 
patterns for various types of aberrations are discussed. Problems of resolution also receive 
attention, including the resolution of objects of periodic structure. Some valuable 
supplementary notes are given. About forty pages are devoted to the theory of electron 
optics, and there are some notes by Hertzberger, including information on the optical 

- qualities of glass. J 

After this survey of the contents of the book little need be said further. It is a pleasure 

to find that the characteristic function method is emphasized and used extensively. More 

publicity on behalf of this method might be given if the notes were used as the basis for a 

printed book more readily available than the mimeographed notes from Brown University. 
L. S. GODDARD. 


Table of Arc sinx. Pp. xx+124. $3.50; Tables of Associated Legendre 
Functions. Pp. xlvit+ 306. $5.00. (Columbia University Press, New York: 
British agents, Scientific Computing Service, Ltd., 23 Bedford Square, 
London -W.C.1, 1945.) 


These volumes are numbered (17) and (18) in the, by now, well-known series of tables 
produced by the New York Mathematical Tables Project, directed by Dr. Arnold N. 
Lowan, and now sponsored by the National Bureau of Standards. They are typical and 
worthy members of the series. 

The purpose of the first volume is to enable 12-decimal values of arc sin x in radian 
measure to be obtained for any value of x from zero to unity. ‘The function is, of course, 


{ fat : 

a tov) 

and since previous volumes have tabulated the transcendental functions (natural logarithm 
and inverse tangent) arising in the integration of rational algebraic functions, the present 
table may be regarded as the first instalment of tables of integrals of functions involving 
the square root of a quadratic. It is believed that work has also been done on the 
companion functions, arc sinh x and arc cosh x, and it is to be hoped that the publication 
of similar tables of these functions may be not long delayed. 

The main Table gives 12-decimal values of arc sin x at interval 0:0001 from x=0 to 
x=0-9890. In this range the greatest attainable accuracy in interpolation is obtained 
by the use of the second differences, or ‘‘ modified ”? second differences, provided. The 
difficult range, for greater values of x, has been treated as follows. First, the interval has 
been reduced to 0:00001. Second differences suffice to x=0-:99750, and modified second 
differences thence to x=0:99900. Fourth differences are given from 0:99900 to 0-99950, 
but no differences in the remaining range. Here interpolation is to be effected by use of 
the formula 


are sin (1—v) = +/(2v) . f(v), 


and f(v) is tabulated at interval 0-00001 from v=0 to v=0:00050. Tables of interpolation 
coefficients for use with the differences provided are included and worked examples of 
their use are given in the introduction. 
The introduction gives also a full account of the methods of computation and checking 
adopted and a bibliography. It is clear from the latter that the present Table has more 
decimals and/or a smaller interval than any of its predecessors, and that it constitutes a 


fundamental table of great vaiue. 
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Considering the importance of the Legendre and Associated Legendre functions 
(e.g., in spherical, ellipsoidal, and toroidal harmonic analysis), it is somewhat surprising 
that tables of these functions have not been more readily available. The present volume 
goes far to fill the gap. It tabulates the functions Pio) and Q'"(x) and their first 
derivatives for integer values of n from 0 to 10 and of m from 0 to 4 (or to n if n<A4), and 
for half-odd-integer values of n from —1/2 to 9/2 and integer values of m from 0 to 4; 
and for real and pure imaginary values of x. 'The argument ranges from 1 (in some cases 0) 
to 10 at interval 0-1. Values of P’” (cos @) and d{P%’ (cos 6)}/d@ for 8 from 0° to 90° 
at interval 1° are also given. Normally six significant figures are tabulated. No pro- 
vision for interpolation is made, but a schedule showing the accuracy attainable by various 
Lagrangian interpolation formulae precedes each main Table. Some supplementary 
tables give exact values and normalization factors. 

It is clear that this volume will be of very great use to physicists and others whose 
work requires numerical values of the functions tabulated in it, and the accuracy should be 
ample for almost all practical applications. W. G. BICKLEY. 


Thermodynamics, by PuILipp Frank. Pp. iv+123+24. (Brown University, 
Providence, Rhode Island, 1945.) 


These notes, in mimecgraphed form, were prepared in connection with one of the 
very valuable summer schools run at Brown University. The scope of the course is, 
roughly speaking, that covered by the ordinary textbook, restricted to the classical (and 
perfectly general) thermodynamics ; the statistical approach is not employed. The 
course, however, is far from being conventional. It is original in being one of the first 
to adopt Caratheodory’s methods, so that, for example, the reader sees clearly that entropy 
is a concept suggested by, and associated with, the first law of thermodynamics, and not, 
as students commonly suppose, only introduced because of the second law. This fact 
was, as it happens, known to Kelvin, but it has tended to slip out of sight. In any case, 
Caratheodory’s treatment is of an altogether higher order of rigour than that of any of 
his predecessors. ‘This is not to say that it uses more difficult mathematics, but its outlook 
is more critical. 

After the bases of the science have been laid, the notes proceed to deal with the 
applications to pure substances, engines, magnetic and electric phenomena ; the third law 
is then introduced and its application to problems of chemical equilibrium described. 
A valuable note on the thermodynamics of gas mixtures is contributed by John A. Goff. 
In this, the methods of Gibbs are supplemented by recent developments, including the 
determination of free energies (or equivalent quantities) from spectroscopic data. 

Whilst we may envy those who were privileged to attend the summer school, and 
may find these notes useful and interesting, it would not be true to say that they would 
make a good book if published as they stand. For that, they would need either to be 
pruned to make more definitely a set of notes, or else expanded further to give those 
full explanations which were doubtless added orally when the lectures were delivered. - 


Hos 


Atomic Artillery and the Atomic Bomb. By J. K. Rogertson. Pp. viii+173. 
(New York: Van Nostrand Co., 1945.) $2.50 


The title of this book appears more sensational than the author really intended, for his 
artillery is, in fact, alpha particles, protons, electrons and the like. The book is intendes 
for the intelligent layman, and covers roughly the work on atomic physics subsequent to 
the discovery of the electron by Sir J. J. Thomson. After an explanation of the nature 
of elements and the significance of atomic weights, the author discusses the concept of 
the atom, proceeds to the electron and discusses methods of accelerating charged particles 
He then mentions the early work on positive rays and the mass spectrograph, gives a shoo 
account of radio-activity, and then comes to gamma rays. . After this, he deals with the 


work on cosmic rays which, in the period 1930-35, gave us the new particles of positron 
and meson (or, as he calls it, the mesotron). 
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The work with which Rutherford’s last period of activity was associated was the trans- 
mutation of elements by bombardment with various ,Particles, and the author gives a clear 
account of the nature of this nuclear chemistry, explaining how both mass number and 
atomic weight must balance in each reaction. From this point his story becomes quite 
modern, and we are introduced to fission and chain reactions, and it is explainzd how 
these could be used to keep a reaction going, and thus provide the world with a new source 
of power. 

As is now well known, this work led not only to a possible source of power, but 
incidentally to a means of production of elements not previously known and, in particular, 
of plutonium. Some short account of this work is given and its application to the atomic 
bomb is mentioned. 

The book, on the whole, seems accurate. There is one slip where positron occurs 
for proton, and the reader is led to suppose that Aston is still alive and working, but these 
seem to be the major slips. 

It is satisfactory to read a book on this subject which does not spin words at great 
length on the sociological and political implications of the work being done, and it is still 
more satisfactory to know that we can now recommend to our non-scientific friends a book 
which they can read for themselves and understand in the knowledge that they will not 
be misled as to the facts. J. H. A. 


Circular and Hyperbolic Functions, Exponential and Sine and Cosine integrals, 
Factorial Function and Allied Functions, Hermitian Probability Functions. 
(British Association Mathematical Tables, Volume I). Second edition. 
Pp. xi+72. (London: Cambridge University Press, 1946.) 10s. 


The collection under review is a welcome second edition of the British Association 
Mathematical Tables, Volume I, originally published in 1931 and now for some years 
out of print. The functions tabulated (together with even differences, sometimes up to 
the sixth) are as follows: $7x[1(1)100 ; 15], sin x, cos x [0:0(0-1)50-0 ; 15 and 0-000 
(0-001)1-600 ; 11], sinh zx, cosh 7x [0-0000(0-0001)0-0100 ;° 15 and 0-00(0-01)4-00 ; 15], 
Sinhx, cosh x [0-0(0-1)10-0 ; 15], Ei(x), —Ei(—x) [0-0(0-1)15-0; 11], Sit), Ci) 


x 
[0-0(0-1)20-0(0-2)40-0 ; 11, 10], x ! [0-00(0-01)1-00 ; 12], { 1+ ie logy) t! dt ae [0-00(0-01) 


sy 


1:00; 10], loge =! [n=1, 2, 3, 4, 0-00(0-01)1-00 and 10-0(0-1)60-0 ; 10, 12], Hh,(x) 
xX 
@ 
alt Hhy,—1 (x) dx where Hhg(a)= | e—** dt [n=0(1)21, x= —7-0(0:1)+6-6 ; 10], Hho(x) 
x 


Hh,(x)/{Hh,(x)}? [—7-0(0-1)+50-0 ; 9 to 3]. In tabular content the new edition thus does 
not differ from the old, but the opportunity has been taken to provide an auxiliary table 
extending the values of the Factorial Function, x !, from 12 to 18 places of decimals and to 
correct eleven errors (of which nine merely affected end-figures by one unit) that have been 
detected. Possessors of the first edition are reminded that they may obtain a list of these 
corrections by applying to the Office of the British Association at Burlington House, 
London W. 1. 

A drastic pruning of the new edition of the tables has removed the historical note, 
and the 30-page introduction which formed such an interesting and informative a feature 
of the first edition. This loss may not be felt acutely by those who, as it were, already 
possess a function and merely require its numerical values, but, from the standpoint of 
general interest and suggestiveness for fruitful application, one experiences a feeling of 
slight regret at the omission of the material illustrating the analytical properties and uses 
of the ‘“ Polygamma ”’ Functions and the repeated integrals of the Error Function (now 
called Hermitian Probability Functions). 

Mathematical tables, although (when accurate) essentially perdurable, are yet susceptible 
to loss of pre-eminence by comparison with subsequent more extensive tables. Few users 
will find the number of significant figures inadequate in the present volume of British 
Association tables, but, as is obvious from the list of contents given above, for many 
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purposes subtabulation of several of the functions would be welcome, although inevitably 
bulky. In this connection it is relevant to remark that since 1931 a new standard for the 
close tabulation of a large number of functions has been established by the publications. 
of the New York Mathematical Tables Project. 

Criticism of a book because it is not a whole library is, however, unfair. The volume 
under review contains a diverse, interesting and most useful collection of tables which yield 
their full value to those prepared to interpolate accurately. The price is very moderate, 
despite the fact that the present binding is not sufficiently stout for its hard future. Itisa 
pleasure to consult tables so elegantly set and printed. M. S. JONES. 


The Music Review, Vol. VII, No.1. (Heffer & Sons, Cambridge.) 5s. 


This number contains an article by B. van der Pol on “‘ Music and Elementary ‘Theory 
of Numbers”, which deals with time and rhythm, absolute and relative pitch, Euler’s 
theory of dissonance, musical scales, relaxation oscillations and variations of musical 
pitch. The author’s investigations on the variation in pitch between orchestras which 
broadcast was described in Vol. VII of the Progress Reports (pp. 33-34). The remaining 
topics are considered either in Helmholtz’s Sensations of Tone or in modern textbooks 
of sound, but the musical physicist will find in this article a useful summary of present 
knowledge in this field, and also some useful observations on the sense of absolute pitch, 
studied by the author, himself a musician with this gift and with the scientist’s sense of 
impartiality. E. G.-R: 


Statistical Thermodynanucs, by ERWIN SCHRODINGER. Pp. 88. (Cambridge: 
The University Press, 1946). 6s. 


This very small book contains a great deal of material, which was originally presented 
as a course of lectures at the Dublin Institute for Advanced Studies, in the spring of 1944. 
It is hardly a reference work to which a man would turn when presented with a specific 
problem which he wished to solve, but rather a background exposition which he should 
read and ponder when he is not busy applying statistical methods to particular problems. 

Schrédinger adopts the attitude that there is only one problem in statistical thermo- 
dynamics, that of determining the distribution of an amount EF of energy over N similar 
systems, which may be only mental copies of one real system. From this starting point 
he considers how the distributions may be enumerated, and in doing so he is led to discuss 
the logical bases of the procedure, and also to sketch the mathematical apparatus normally 
used. He is more interested in the difficulties of the theory than in those parts where 
all is straightforward, and consequently has written a most stimulating book, which is 
well worth its price. Incidentally, though produced under war-time conditions, the 
book is pleasant to handle and to look at—a tribute to the skill of the printers. Among 
matters which are discussed are the setting equal to zero of the constant entropy which 
remains in a crystalline system at absolute zero, and the relationship between the classical, 
the Bose-Einstein and the Fermi-Dirac statistics. All these are subjects on which a great 
deal has been written before, but the author’s peculiar quality of illuminating those matters 
on which he writes makes his treatment very well worthy of study. J. Haare 


CORRIGENDUM 


“ Crystal structure of double oxides of the perovskite type’’, by HELEN D. Mecaw 
(Proc. Phys. Soc., 58, p. 133). 


The address given under the author’s name should read “ Philips’ Lamps, Limited, 
Mitcham, Surrey ; now at Birkbeck College, London”’. 


THE PROCEEDINGS OF THE PHYSICAL SOCIETY 


THE PHYSICAL SOCIETY 


MEMBERSHIP 
Membership of the Society is open to all who are interested in Physics; 


Fettowsnip. A candidate for election to Fellowship must as a rule be recommended 
by three Fellows, to two of whom he is known personally. Fellows may attend all meetings 
of the Society, are entitled to receive Publications 1, 3, 4 and 5 below, and may obtain the 
other publications at much reduced rates. 


STuDENT MemsBersHip. A candidate for election to Student Membership must be 
between 18 and 26 years of age and must be recommended from personal knowledge by 
a Fellow. Student Members may attend all meetings of the Society, are entitled to receive 
Publications 1, 3 and 4, and may obtain the other publications at much reduced rates. 


Books and periodicals may be read in the Society’s Library, and a limited number may 
be borrowed by Fellows and Student Members on application to the Honorary Librarian. 


Fellows and Student Members may become members of the Colour Group, the Optical 
Group and the Low-Temperature Group (Discussion Groups formed in the Society) with- 
out payment of additional annual subscriptions. 


PUBLICATIONS 


1. The Proceedings of the Physical Society, published six times annually, contains original 
papers, lectures by specialists, reports of discussions and of demonstrations, and reviews. 
Universal Decimal Classification Index Slips for the Proceedings are issued gratis to those 
members who ask for them in advance. 


2. Reports on Progress in Physics, published annually, is a comprehensive review by 
qualified physicists. 

3. The Catalogue of the Physical Society’s Annual Exhibition of Scientific Instruments 
and Apparatus. ‘This exhibition, which is held normally in January of each year, is recog- 
nized as the most important function of its kind. 


4. The Agenda Paper, issued at frequent intervals during the session, informs members of 
the programmes of future meetings and business of the Society generally. 


5. Physics Abstracts (Science Abstracts A), published monthly in association with the 
Institution of Electrical Engineers, covers practically the whole field of contemporary 
physical research, 

6. Electrical Engineering Abstracts (Science Abstracts B), published monthly in association 
with the Institution of Electrical Engineers, covers the whole field of contemporary research 
in electrical engineering. 

7. Special Publications, critical monographs and reports on special subjects prepared 
by experts or committees, are issued from time to time. Particulars of these publica- 
tions are given in the advertisement pages of the Proceedings. 


MEETINGS 


At approximately fortnightly intervals throughout each normal session, meetings are 
held for the reading and discussion of papers, for lectures, and for experimental demonstra- 
tions. Special lectures include: the Guthrie Lecture, in memory of the founder of the 
Society, given annually by a physicist of international reputation ; the Thomas Young 
Oration, given biennially on an optical subject ; the Charles Chree Address, given biennially 
on Geomagnetism, Atmospheric Electricity, or a cognate subject ; and the biennial Ruther- 
ford Memorial Lecture. A Summer Meeting is generally held each year at a provincial centre, 
and from time to time meetings are arranged jointly with other Societies for the discussion 
of subjects of common interest. 


SUBSCRIPTIONS 


Fellows pay an Entrance Fee of £1. 1s. and an Annual Subscription of £2. 2s. Student 
Members pay only an Annual Subscription of 10s. 6d. No entrance fee is payable by 
a Student Member on transfer to Fellowship. 


Further information may be obtained from the Honorary Secretary (Business) 
at the office of the Society, 


1 Lowrner GarpDENs, Prince Consort RoapD, Lonpon S.W.7 


THE PROCEEDINGS OF THE pevereal cece | 
ELECTRICAL MEASURING 
INSTRUMENTS © 
OF THE HIGHER GRADES 


ERNEST TURNER 
ELECTRICAL INSTRUMENTS 


LIMITED 
CHILTERN WORKS 
HIGH WYCOMBE 
BUCKS 


Telephone : 


Telegrams : 
High Wycombe 1301/2 


Gorgeous, High Wycombe 4 


